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1.1  lack  around  In.f?nitl*1 9fl 

Tha  International  Gravity  Standardization  Net  1971  (XGSN  71) 

produced  a  world-wide  gravity  reference  system  of  gravity  base  stations 

having  a  standard  error  for  each  station  less  than  0.1  ogal, 

-5  2 

1  mgal  r  10  ne/ a  ,  using  absolute,  pendulum,  and  gravity  meter 
measurements  (Morelli,  et  al.,  1974].  See  Table  1  for  a  summary  of  the 

i 

measurements  used  in  the  XGSN  71  which  indicates  that  the  most  widely 
used  gravity  meter  in  the  XGSN  71  was  the  LaCoste  1  Romberg  gravity 
meter.  However,  since  the  XGSN  71  results  were  adopted  by  the 
International  Union  of  Geodesy  and  Geophysics  (XUGG)  in  Moscow,  1971, 
and  subsequently  published  (Morelli,  et  al.,  1974],  there  has  developed 
an  ever  Increasing  desire  to  obtain  better  values  of  gravity  for 
gravity  base  stations. 

Given  an  existing  gravity  base  station  network,  the  only  way  the 
standard  error  for  the  stations  can  be  improved  is  by  using  new 
information.  This  new  information  can  be  in  the  form  of  either  new 
measurements  or  Improved  modelling  of  the  relationship  between  the 


a 


Table  1  —  Suaaary  of  data  uaad  In  ISSN  71. 
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Nuabar  of 
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4  trips 

North  Amor 1  can 

a 

S  trips 

Hast  am 

s 
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Tha  quantltlas  In  paranthosos  rapraaant  tha  approxlaata  nuabar  of 
aaasuroaonts  aada  with  aach  class  of  Instruaanta. 
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If  new  MtsurMOTti  are  made,  they  can  be  classified  as  either 

absolute  or  rolativa  in  nature.  Absolute  aieasureatents  eade  with 

absolute  Measuring  apparatus  are  used  to  determine  the  value  of  gravity 

which  is  the  magnitude  of  the  vertical  gradient  of  the  goopotential 

[Mueller  and  Reekie,  lHi,  pp  43—44],  at  a  station  by  measuring  the 

time  it  takes  an  object  to  fall  a  specific  distance.  The  value  of 

gravity  obtained  is  the  resultant  of  gravitation  and  the  centrifugal 

force  caused  by  the  rotation  of  the  earth  [Mueller  and  Reekie,  1944, 

pp  43—44].  The  physical  dimension  of  gravity  Is  an  acceleration  Mith 

its  magnitude,  in  the  geodetic  community,  given  in  either  gal,  mgal  or 

2 

ygal  where  1  gal  =  0.01  m/s  ,  1  gal  s  1000  mgal,  and 

1  mgal  =  1000  ygal.  The  value  of  gravity  varies  on  the  earth's  surface 
from  about  970  gal  at  the  equator  to  about  903  gal  at  the  poles 
[Heiskanen  and  Moritz.  1947,  p  40]. 

Relative  measurements  are  wade  using  gravity  meters.  Brevity  meters 
do  not  have  the  direct  capability  of  measuring  the  absolute  value  of 
gravity  but  they  are  used  to  measure  the  difference  in  gravity  between 
two  stations.  For  gravity  base  station  networks,  the  most  commonly 
used  gravity  meter  is  the  LaCoste  0  Romberg  *B*  gravity  meter 
manufactured  by  LaCoste  0  Romberg,  Inc.,  Austin,  Texas.  This  gravity 
meter  was  designed  to  be  able  to  be  used  anywhere  on  the  surface  of  the 
earth,  tines  the  approximate  gravity  difference  between  the  equator 
and  the  poles  is  about  0  gal,  in  order  to  insure  world-wide  measuring 
eapabilites,  the  gravity  meter  has  a  measuring  range  of  approximately 
7  gal  and  is  adjusted  to  work  within  the  range  of  absolute  gravity  from 
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•pprexlMttly  977  gal  to  994  gal. 

It  Is  obvious  that  If  tho  gravity  valuos  of  a  few  stations  were 
desired.  and  If  tlaa  and  funds  prosantad  no  problaa.  tho  bast  aothod  to 
obtain  tha  stations'  gravity  valuas  would  ba  to  aake  absoluta  gravity 
dat anal  nations  at  aach  station  using  althar  a  transportabla  absoluta 
gravity  aaasurlng  apparatus  or  by  Installing  a  paroanant  absoluta 
gravity  aaasurlng  apparatus.  however.  dua  to  rastrletlons  on  tlaa 
and/or  funds,  this  Is  not  always  faaslbla  whan  tha  valuas  of  gravity  at 
a  largo  nuabor  of  stations  ara  doslrod.  Zt  takas  approxlaataly  2  to 
4  days  to  sat  up  and  to  aaka  aaeh  absoluta  gravity  dotaral nation. 

Thus.  If  only  ona  absoluta  gravity  aaasurlng  apparatus  wars  usad.  It 


would  taka  about  S  yoars  to  astabllsh  a  gravity  basa  station  notwork  o 
300  stations,  assualng  tha  absoluta  apparatus  workad  continuously.  To 
raduca  tha  tlaa  required.  gravity  basa  station  notworks  ara  astabllshoi 


In  which  tho  gravity  valuos  of  stations  ara  rasolvod  froa  ralatlva 
gravity  aatar  obsarvatlons  and  tha  control  Is  provided  by  a  fow 
stations  whoso  gravity  value  has  bean  deteralned  by  absolute  gravity 
aaasurlng  apparatus. 

Tha  standard  errors  that  ean  bo  associated  with  tho  gravity  valuas 
of  tha  stations  In  such  a  notwork  depend,  not  only  on  tha  distribution 
and  quality  of  tho  gravity  aatar  obsarvatlons  and  tha  absoluta  gravity 
value  deteral not Ions  aada,  but  also  on  how  wall  tho  relationship 
between  tha  obsorvad  quantities  and  tha  derived  quantltes  ean  bo 
represented.  I.e.  aathaaatleally  nodal led. 


1.2  Description  of  Prtunt  Atudv 

The  nature  of  this  study  Is  to  Investigate  tho  behavior  of  the 
LaCoste  4  Roaberg  *0*  gravity  aeter  and  develop  a  aatheaatieal  aodol 
which  approxlaates  this  behavior.  In  ordor  to  test  the  aatheaatieal 
aodels  developed*  data  froa  the  United  Stated  Gravity  tase  Station 
Network  will  be  used.  The  aajorlty  of  this  data  was  obtained  along  tho 
Ml d-Cont 1 nent  Calibration  Line  which  Is  located  along  the  eastern  side 
of  the  Rocky  Mountains  froa  Now  Mexico  to  Montana.  Along  this 
Mid-Continent  Calibration  Line*  eleven  absolute  gravity  stations  wero 
established  at  Intervals  of  approxlaately  200  agal  to  provide  control. 
The  vast  aajorlty  of  the  gravity  aeter  observations  were  aade  with  a 
nuaber  of  LaCoste  S  Rea berg  *0*  gravity  eaters  with  the  reaalnder  being 
aade  with  two  LaCoste  A  R saber g  *0*  gravity  aeters*  VD17*  and  *D4S*. 
LaCoste  A  Roaberg  *D*  gravity  aeters  work  over  a  very  Halted  gravity 
difference  of  approxlaately  200  agal  but  can  be  adjusted  to  work 
anywhere  In  the  range  of  the  LaCoste  A  Roaberg  *G*  gravity  aeter. 

In  order  to  understand  what  causes  the  LaCoste  A  Roaberg  *0*  gravity 
aeter  to  behave  as  It  does*  a  review  of  the  asseably  and  testing 
procedures  used  In  Its  construction  will  be  done.  Most  of  tho 
Inforaatlen  about  the  gravity  aeter  was  based  on  first  hand  accounts 
obtained  during  a  visit  In  Deeeaber*  1*00,  to  the  LaCoste  A  Roaberg* 
Inc.  facilities  In  Austin*  Texas. 

Free  the  inforaatlen  obtained  about  what  goes  into  the  construction 
and  production  of  a  *0*  gravity  aeter*  atteapts  will  be  aade  to  develop 
a  aore  represented ve  aatheaatieal  aodel  of  the  Instruaent's  behavior. 


« 

Another  aroa  to  bo  explored  Is  the  eethod  that  could  be  uaed  to 
Indicate  where  new  aeasureaonts*  abaolute  or  relative*  ahould  be  aado 
in  order  to  laprove  any  existing  gravity  base  station  network. 

1.3  Review  of  Previous  »tud1es 

Various  aatheaatleal  aodels  have  been  proposed  and  used  to  aodel  the 
behavior  of  the  LaCoste  A  Roaberg  *0*  gravity  aeter  [herein  et  al. 
1171;  Torge  and  Kanngleser*  1979].  There  are  basically  two  different 
types  of  aodels.  One  type  of  aodel  uses  as  Its  observable  the  value  In 
allllgals  of  the  observed  gravity  aeter  counter  reading  Interpolated 
froa  the  Calibration  Table  1  supplied  with  each  gravity  aeter  [llotlla. 
1974] ;  the  other  type  ef  aodel  uses  the  difference  In  value  In  allllgal 
between  two  consecutive  gravity  aeter  observations  as  It  observable 
[McConnell  and  Gantar*  1974;  Whalen,  1974;  Torge  and  Kanngleser,  1979]. 
Each  aodel  differs  further  In  what  paraaeters  are  Included  and  how  they 
are  related,  teae  aodels  Include  possible  relationships  for  linear 
gravity  aeter  drift  [llotlla*  1974).  One  aodel  even  postulates  e 
relationship  Involving  the  square  of  stations'  gravity  values  [Torge 
end  Kanngleser*  1979]. 

One  thing  that  Is  esaasn  for  aodels  proposed  to  date  Is  that  the 
Calibration  Table  1  which  Is  used  to  convert  gravity  aeter  counter 
readings  to  their  values  In  allllgal  Is  assuaed  to  be  eorrect  except 
for  a  linear  scale  factor  which  needs  to  be  applied  to  all  the  values 
In  allllgal.  Models  have  been  proposed  that  Include  additional  higher 
order  scale  factor  teres  [llotlla*  1974|  Torge  and  Kanngleser*  1979]. 


Du*  to  th*  construction  of  th*  LaCost*  A  Romberg  *  O’  gravity  motor, 
thoro  Is  •  possibility  of  •  periodic  variation  In  the  counter  readings 
being  Introduced  due  to  Imperfections  In  th*  gear  reduction  system 
[Klvlnieml,  1974],  Th*  amplitude  of  this  effect  has  been  estimated  to 
be  as  large  as  0.04  mgal  [Tor go  end  Kanngieser,  1979].  Th*  attempts  to 
solve  for  this  periodic  effect  have  resulted  In  no  clear  conclusion  of 
Its  existence.  On*  reason  might  be  that  th*  Calibration  Table  1  Is 
being  used  as  a  standard  and  assumed  to  be  error  free,  thus  masking  th* 
existence  of  the  periodic  effect.  Another  possibility  is  that  the 
mathematical  models  proposed  by  Torg*  and  Kanngieser  [1979]  are  not  be 
appropriate.  Further,  the  control  provide  by  th*  absolute  gravity 
si  tea  might  not  be  distributed  appropriately  and  be  of  high  enough 
accuracy  to  solve  for  th*  periodic  effects. 

To  better  understand  th*  requirements  of  th*  mathematical  model 
needed  to  represent  the  behavior  of  th*  LaCoste  A  Romberg  *6*  gravity 
meter.  It  Is  necessary  to  know  how  th*  gravity  meter  Is  constructed  and 


how  it  works. 


CHAPTIR  TWO 


kAfifisu  *  BoniHg  2sz  cmvity  am 

2.1  a>»«l«Bwtnt  snd  History 

The  LaCoste  A  leubtrg  "0"  gravity  attar  i»  a  rtlativa  gravity 
aiaaauring  davlca  daalgnad  to  ba  usad  anywhere  on  tha  land  surface  of 
tha  aarth.  To  aehiava  world— wide  measuring  capability,  tba  "0**  gravity 
eater  Must  ba  abla  to  measure  diffarancaa  in  gravity  as  largo  as  5  gal. 
To  satisfy  this  requirement,  tha  "0"  gravity  aatar  was  dasignad  to  ba 
abla  to  aaasura  diffarancaa  of  approximately  7  gal  and  adjusted  to  work 
within  tha  range  of  absolute  gravity  froa  approxiaately  977  gal  to 
9<4  gal.  Tha  saallast  difference  that  can  bo  datarainad  directly  froa 
tha  "0"  gravity  aatar  observations  corresponds  to  approxiaately 
10  ugal. 

Tha  first  LaCoste  0  loabarg  land  gravity  aatar  was  Manufactured  in 
1919  and  was  tha  pradoeassar  of  tha  present  "0"  gravity  aatar. 

Coaparad  to  tha  present  "0"  aadals  being  aenufactured,  the  first  land 
gravity  aatar  was  auch  bigger  and  heavier.  Tha  currant  "0"  gravity 
eaters  being  aanufaeturod  and  used  are  approxiaately  20  ea  in  length, 

10  ea  in  width  and  20  ce  in  height  with  each  weighing  approxiaately 
0.2  kilograas  excluding  batteries  tLaCoste  0  Romberg,  1901].  To  data, 
in  tha  neighborhood  of  000  LaCoste  0  Romberg  "0"  gravity  eaters  have 
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been  built,  with  the  currant  production  loyal  running  about  30  to  40 
■otars  par  yaar.  Each  ”0"  gravity  natar  raquiras  approximately  3  to  12 
weeks  to  construct  and  tost  from  the  start  of  motor  building  until  tha 
gravity  motor  is  dalivarad  to  tha  customer.  However,  due  to  tha 
currant  largo  demand  for  tha  **0"  gravity  motor  and  tha  shortage  of 
qualified  gravity  mater  builders,  tha  load  time  is  running  around  two 
years  for  tha  delivery  of  a  "0"  gravity  mater  [Perry,  1930,  private 
eoaimunl  cation]. 

2.2  Changes  Midi  la  ZSH  OraviAv  Mater 

Tha  actual  internal  workings  of  the  "0"  gravity  mater  have  basically 
remained  the  same  since  the  first  one  was  manufactured  in  1939.  Tha 
major  design  changes  that  have  occurred  have  been  cosmetic  'to  the  meter 
case.  The  changes  were  made  to  reduce  the  sice  and  weight  of  the 
gravity  meter  and  to  allow  for  electronic  Improvements  and  options  such 
as  electronic  readout  capability  which  is  one  of  the  options  available 
on  the  gravity  meters.  However,  starting  with  meter  ,6-4SSv,  a  new 
gear  reduction  system  <aaar  box)  was  installed  into  the  ”0”  gravity 
meter.  This  gear  box  changes  the  gear  ratios  in  addition  to  changing 
the  type  of  gear  bex  [Perry,  1910,  private  communication].  A  detailed 
description  of  what  this  change  affected  is  given  in  Section  2.4.2. 

2.S  What  Does  a  Brevity  Meter  £gf 

An  observation,  0.  made  with  a  "0"  gravity  meter  at  a  station  is 
related  to  the  absolute  gravity  value  at  the  station  by  the  following 
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relation: 

f(0)  ♦  C  +  S  -  G  <2-x> 

Mhoro 

f  (0)  —  a o«o  functional  relationship  of  the  gravity  meter's 
observation, 

C  —  corrections  that  make  the  function  f(0)  independent  of  the 
epoch  of  the  observation, 

t  —  an  unknown  offset  value  that  must  be  added  to  obtain  the 
correct  absolute  gravity  value  for  the  station, 

C  —  absolute  gravity  value  for  the  station. 

Assuming  that  the  unknown,  S,  is  constant  for  an  instrument  during  a 
particular  time  period,  and  the  functional  relationship,  fCO),  does  not 
change  during  that  time  period,  the  gravity  difference  between  two 
consecutive  observations  at  station  1  and  station  j,  can  expressed  by 
f(0,)  -  f(0j)  +  C,  -  Cj  -  G{  -  Gj  (2.2) 

where  the  subscript  i  refers  to  events  at  station  i  at  one  epoch  and 
subscript  j  refers  to  events  at  station  )  at  another  epoch. 

eliminating  the  unknown  offset,  >,  from  equation  (2.2)  does  not  mean 
that  the  value  of  t  is  net  needed.  By  this  method,  the  value  of  the 
unknown  offset,  t,  remains  unknown  but  dees  not  enter  directly  into 
equation  (2.2).  However,  the  value  ef  t  is  still  needed,  as  can  be 
seen  in  equation  (2.1),  to  determine  the  gravity  value  of  a  station. 

In  order  to  determine  e  value  for  B,  it  requires  that  at  least  one 
station's  gravity  value  in  the  network  be  known.  This  implies  that  no 
matter  how  many  equations  similar  to  equation  (2.2)  are  formed,  the 
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system  of  equations  will  have  a  rank  deficiency  of  at  least  one  as  the 
result  of  the  unknown  value  of  offset  S.  The  exact  rank  deficiency  of 
the  system  will  depend  on  the  functional  relati onshi p,  f(0>,  used.  The 
rank  deficiency  of  the  system  determines  the  minimum  number  of  stations 
whose  value  of  gravity  is  required  to  be  known. 

2.4  Internal  Components  gf  !&£  SC  gravity  Meter 

The  LaCoste  4  Romberg  "0"  gravity  meter  can  be  thought  of  as 
consisting  of  two  major  components,  the  meter  case  and  the  meter  box. 
The  meter  case  which  houses  the  meter  box  consists  of  the  external 
casing  Including  l££.  insulating  material,  electrical  components. 
llV<H1.ng  screws,  temperature  probe  and  the  heater  ggx.  The  levelling 
bubbles,  however,  are  part  of  the  meter  box  and  not  part  of  the  meter 
case. 

The  meter  box  is  a  mechanical— optical  device  which  can  be  thought  of 
consisting  of  four  major  inter-connected  assemblies:  oear  train. 
Mflgyri.flq  jflcm.  lever  linkage  and  gpt.lffll  *  Yalta  - 

2.4.1  Meter  Case 

The  meter  case  houses  the  meter  box.  The  meter  case  is  Insulated  to 
protect  the  meter  from  changes  in  the  ambient  temperature.  Changes  in 
the  operating  temperature  of  the  meter  box  have  a  marked  effect  on  the 
gravity  meter  readings  [Klvlnlemi,  1974].  To  insure  a  constant 
operating  temperature,  the  meter  box  is  Installed  In  a  heater  box  which 
requires  a  small  amount  of  electrical  power  to  keep  the  Instrument  at 
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operating  temperature.  The  optimum  operating  temperature  to  be 
ntained  is  determined  from  teat  procedurea  after  which  the  heater 
ia  adjuated  to  maintain  that  operating  temperature  for  the 
trument.  A  temperature  probe  la  1  natal led  in  the  meter  case  to 
mit  the  uaer  to  verify  that  the  inatrument  ia  at  ita  operating 
perature  and  that  the  heater  box  la  working  properly.  The 
ulation  alao  acta  as  a  shock  absorbing  material  in  case  the  gravity 
er  were  to  be  accidentally  jarred  or  dropped. 

The  gravity  mater  can  be  used  to  make  consistent  observations  only 
er  the  instrument's  operating  temperature  has  bean  attained  and 
ntained  for  a  length  of  time.  When  this  occurs  the  Instrument  is 
d  to  be  on-heat .  The  recommended  length  of  time  of  being  on— heat 
ore  observations  should  be  made  is  about  four  hours  (Perry,  1980, 
vate  communication],  if  the  power  is  Interrupted  for  any  length  of 
e  and  the  temperature  of  the  Instrument  falls  below  its  operating 
serature,  the  Instrument  is  said  to  be  off-heat .  If  an  instrument 
i  off-beat,  it  must  be  put  back  on-heat  before  it  can  be  used  to 
i  additional  observations.  Any  gravity  meter  differences  determined 
i  the  meter  was  off—heat  must  not  be  considered  as  part  of 
trvation  set.  This  means  that  the  gravity  meter  must  be  on— heat 
ng  its  transportation  between  stations  when  observations  are  being 
>. 

he  top  of  the  meter  ease  is  removable  to  permit  the  installation  of 
meter  box.  In  addition,  on  the  top  will  be  a  name  plate  which 
t  if  lea  the  Instrument  and,  generally,  below  the  thermometer  opening 
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will  be  the  value  of  tha  currant  null  or  reading  Una  'for  the 
Instrument. 

Tha  meter  case  contains  the  levelling  screw  by  which  the  gravity 
meter  Is  levelled.  For  a  period  of  time,  the  LaCoste  3  Romberg  "0" 
gravity  meters  were  manufactured  with  the  thumb  screws  used  for 
levelling  the  gravity  meter  located  on  the  bottom  of  the  meter  case, 
lut  now,  the  Instrument  Is  being  manufactured  so  the  levelling  can  bo 
adjusted  via  knobs  that  extend  above  the  top  of  the  meter  case.  This 
modi fl elation  does  not  affect  the  behavior  of  the  Instrument  but  makes 
the  levelling  of  the  Instrument  more  convenient  and  reduces  the 
possibility  of  jarring  or  moving  the  Instrument  when  It  Is  being 
levelled  [Ferry,  1930,  private  communication]. 

In  addition,  all  the  electrical  connections  for  the  Instrument  are 
housed  In  the  meter  case.  These  Include  the  connections  for  the  power 
supply  to  operate  the  heater  box  and  lamps  and  any  optional  electronic 
devices  such  as  the  electronic  readout.  The  electronic  readout  Is 
really  part  of  the  meter  box  since  It  basically  consists  of  a  sot  of 
capacitor  plates  Installed  above  and  below  the  beam  and  a  nulling 
meter.  As  the  beam  moves  between  these  plates,  the  change  In  the 
capacitance  Is  recorded  on  the  nulling  meter  Installed  In  the  meter 
ease  top.  This  nulling  meter  Is  then  adjusted  so  when  the  Instrument 
Is  In  the  null  position,  the  nulling  meter  will  be  In  Its  center 
position  (Hemlngson,  1930,  private  communication). 
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2.4.2  Hater  Box  Saar  Train  Assembly 

Tha  gaar  train  assembly  shewn  in  Figure  1  consists  of  the  dial,  dial 
shaft .  counter  and  tha  gear  box.  Tha  dial*  which  is  turned  to  null  tha 
instrument,  has  100  equal  divisions  Marked  on  it  with  each  division 
correspondi ng  to  0.01  counter  units.  The  counter  is  attached  to  tha 
dial  shaft  and  is  used  to  keep  track  of  tha  number  of  rotations  of  the 
dial  shaft  with  one  counter  unit  equivalent  to  one  cowplete  rotation  of 
the  dial  shaft.  Tha  counter  records  in  0.1  counter  units  fro*  0.0  to 
4999.9.  Physical  stops  within  tha  counter  prevent  tha  dial  from  being 
rotated  outside  this  range.  On  the  end  of  the  dial  shaft  opposite  the 
dial  is  a  toothed  gear  which  drives  the  gears  in  the  gear  box. 

There  are  two  different  types  of  gear  boxes  that  have  been  installed 
in  the  "G"  gravity  meters.  The  original  gear  box  installed  in 
instruments  prior  to  *B-4SS*  is  referred  to  as  the  old  gear  box.  The 
old  gear  box  used  a  floating  pivot  gear  system.  In  this  system,  the 
dial  shaft  with  its  17  tooth  gear  drove  a  134  tooth  floating  pivot 
gear.  The  floating  gear,  in  turn,  had  a  20  tooth  smaller  gear  which 
drove  a  ISO  tooth  gaar  on  a  shaft  to  which  the  measuring  screw  was 
attached  as  shown  in  Figure  2.  The  floating  pivot  gear  was  held  in 
contact  with  the  dial  gear  and  the  measuring  screw  gear  by  spring 
tension. 

The  gear  box  installed  in  meters  built  since  meter  *8-430*  uses  a 
fixed  pivot  gear  system  and  is  referred  to  as  the  new  gear  box.  The 
new  (jear  box  uses  a  30  tooth  gear  at  the  end  of  the  dial  shaft  which 
drives  a  220  tooth  fixed  pivot  gear.  The  fixed  pivot  gear  in  turn  Has 
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Figure  2  —  Sehewatie  illustration  of  tho  eld  gear  box  used  in  the 
LaCosto  A  Roaberg  "8"  gravity  Motor. 
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•  SO  teeth  gear  which  drives  a  S00  teeth  gear  on  the  shaft  to  which  the 
■assuring  screw  Is  attached  as  shewn  In  Figure  S. 

Mlth  the  eld  gear  ben,  1200  rotations  of  the  dial  are  required  to 
■aka  the  Measuring  scraw  rotate  17  tines,  a  ratio  of  approximately 
70.041:1;  with  the  new  gear  box,  220  rotations  of  the  dial  results  In  2 
rotations  of  the  Measuring  screw,  a  ratio  e f  approximately  72.222:1. 

Just  because  a  gravity  mter  originally  had  an  eld  gear  box 
Installed  In  It  dees  net  Imply  that  It  will  always  have  an  old  gear 
box.  If  the  gravity  mter  la  returned  to  the  factory  fer  repairs,  and 
the  gear  box  needs  to  be  replaced,  a  new  gear  box  Might  be  used  as  a 
replacement.  Therefore,  It  Is  very  Important  to  knew  what  components 
are  presently  Installed  In  the  gravity  meters  being  used  because  they 
effect  the  modelling  of  the  gravity  meter's  behavior.  If  there  Is  any 
doubt,  the  manufacture's  leg  on  the  construction  of  each  gravity  meter, 
which  Is  maintained  by  LaCoste  4  Romberg,  Inc.,  should  be  consulted. 

*.*.2  Hater  ftfiS  Measuring  term  fttamhta 

The  measuring  screw  moves  within  a  hollow  shaft  which  has  threaded 
fingers  at  the  end  furthest  from  the  gear  box  as  shown  In  Figure  4. 

The  rotation  of  the  dial  causea  an  angular  motion  of  the  measuring 
screw.  This  angular  motion  of  the  measuring  screw  Is  converted  to  a 
linear  motion  by  the  measuring  screw  threaded  being  In  contact  with  the 
stationary  threaded  fingers.  The  maximum  linear  motion  of  the 
measuring  scraw  In  a  "0"  gravity  meter  Is  on  tho  order  of  20  mm  and 
this  motion  Is  accomplished  In  less  than  100  turns  of  the  measuring 


DIAL  SHAFT 


Schematic  Illustration  of  tho  now  gaar  bon  usod  In  tho 
LaCosto  t  Sombor g  "8"  gravity  motor. 
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serew  [Ferry,  1990,  private  coaaunleatlon] . 

At  the  and  of  the  ■••■urlng  screw  furthest  froa  the  gear  box,  a 
hardened  eetal  jewel  In  the  shape  of  a  denut  Is  set  by  a  press  fit  Into 
the  Measuring  serew  as  shewn  In  Figure  S.  The  dlaaeter  ef  the  hele  In 
the  jewel  Is  less  than  1  aa  [Ferry,  19*0,  private  coaaunleatlon). 

There  are  feur  set  screws  In  the  Measuring  serew  which  can  be  used  fer 
Miner  centering  adjustaents  ef  the  jewel.  Since  there  Is  a  spherical 
aetal  ball  en  the  end  ef  the  lever  linkage  which  aakes  contact  with  the 
jewel.  It  Is  laportant  that  the  jewel  la  shaped  and  positioned  such 
that  the  spherical  aetal  ball  Is  always  In  contact  with  the  edge  ef  the 
jewel's  hele.  Zf  this  Is  net  the  ease,  then  the  uni  fora  rotation 
aetlen  ef  the  screw  eeuld  be  translated  Into  a  non-uni  fora  aotlon  ef 
the  lever  linkage  resulting  In  the  difference  between  counter  readings 
fer  a  given  gravity  difference  net  being  eenatant.  The  actual 
difference  In  oeunter  reedlngs  would  then  depend  en  the  starting 
position  ef  the  dial  fer  each  counter  reading. 

2*4.4  Hater  Sew  Lever  Llnkaee  Aaseablv 

The  lever  linkage  consists  ef  a  lower  lever,  connecting  linkage, 
upper  lever,  sere  length  spring,  beaa  and  beaa  weight  as  shewn  In 
Figure  4.  This  systaa  Is  the  heart  ef  the  gravity  Meter.  Many 
Individual  parts  aust  be  asseabled  to  create  this  delicate  systaa.  The 
connecting  linkage,  fer  ewaaple,  consists  ef  a  nuaber  ef  flat  springs 


screw  el 
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Plguro  4  —  tellMHtte  Illustration  of  tKo  Measuring  screw  usod  In  tho 
LoCosto  A  ftoaberg  "0"  gravity  Motor. 
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I: 

*  FI  guru  ft  —  Sehaaatle  Illustration  o-f  oonnaetlon  botuoan  tha  aaaaurlng 

aeraw  and  tha  lavar  llnkago  uaad  In  tha  LaCosta  ft  tombarg 
"O’*  gravity  aatar. 


a 
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Tli«  measuring  screw's  jewel  fa  designed  to  be  In  continuous  contact 
with  the  spherical  metal  bell  attached  to  the  lower  lever.  The 


spherical  eetal  ball  has  been  known  to  break  loose  from  its  support 
pest  in  which  case  the  reugh  edges  of  the  support  post  were  in  contact 
with  the  jewel.  This  causes  the  difference  in  readings  between  gravity 
stations  te  act  very  erratic  which  weuld  be  a  definite  indication  that 
the  gravity  eater  needed  te  be  repaired. 

The  beam  with  the  beam  weight  at  one  end  is  penal tted  to  wove  only  a 


few  thousands  ef  an  inch  in  the  lateral  and  horizontal  direction  before 
it  encounters  physical  steps  (Hemingaon,  1930,  private  eoaaaunication]. 
There  is  an  arrestment  kneb  which  penal ts  the  beaai  to  be  clamped 
against  the  stops  so  that  damage  te  the  systaai  can  be  minimized  during 
the  transportation  of  the  Instrument. 

The  beam  weight  in  addition  to  providing  necessary  mass  and  balance 
for  the  beam,  provides  the  means  of  calibrating  the  instrument.  How 
this  is  accomplished  is  explained  in  section  2.5. 

2.4.5  Meter  la*  PgUSll  tvstem 

The  link  between  the  lever  linkage  and  the  optical  system  is  by 
means  of  what  is  called  the  ladder,  which  is  suspended  from  the  bottom 
of  the  beam  near  the  beam  weight.  This  ladder  consists  of  two  posts 
with  a  number  of  thin  wire  steps  strung  between  the  posts  as  shown  in 
Figure  7. 

A  set  ef  prisms  direct  light  which  has  passed  through  the  ladder 
onto  an  etched  scale  mounted  on  the  meter  box.  The  eyepiece  is  focused 
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on  this  etched  seals  which  produces  an  imago  similar  to  tho  one  shown 
In  Figure  3. 

The  cross-hai r  viewed  through  the  eyepiece  is  actually  the  shadow  of 
one  of  the  thin  wire  steps  of  the  ladder.  Adjustments  are  made  so  only 
one  step  can  be  viewed  through  the  eyepiece  tHemingson,  1930,  private 
comaiunication].  The  step  that  is  actually  viewed  depends  on  how  the 
prisms  and  light  source  are  installed  and  the  ladder  is  constructed. 

To  Insure  that  a  step  can  always  be  viewed,  the  ladder  is  constructed 
with  many  steps. 

The  physical  stops  are  ed justed  so  the  cross— hair  will  move 
approximately  six  to  seven  scale  divisions  either  side  o-f  the  reading 
or  null  line.  The  reeding  line  for  the  Instrument  Is  determined  during 
the  construction  of  the  instrument  following  a  procedure  which  allows 
the  builder  to  deduce  when  the  beam  is  in  the  horizontal  or  null 
position.  Mhen  the  actual  null  position  does  not  correspond  exactly 
with  an  etched  scale  division,  the  nearest  etched  scale  division  is 
selected  es  the  reeding  line. 

Two  bubble  levels  mounted  perpendiculer  to  each  other  are  Installed 


on  the  meter  box.  The  levels  used  are  generally  AO  second  levels  but 
the  customer  cen  request  30  second  levels  be  Installed  [Hemingson, 
1930,  private  communication].  These  levels  are  adjusted  so  when  they 


ere  centered,  the  beam  is  horizontally  positioned  between  its  physical 
stops  when  it  is  in  its  null  position. 
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EXAMPLE 

READING  LINE  »  2.3 


r*  *  -*  •choaatic  viow  through  tho  ayapioco  of  tho  LaCoato  *  Hamburg 
"0"  gravity  motor. 
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2.5  Calibration  of  tha  "0"  Oravl  tv  Matar 

During  tha  construction  of  tha  gravity  matar.  tasts  ara  mada  on 
various  systams  and  assambllas  In  an  attempt  to  ansura  soma  typa  of 
uniformity  In  tha  operational  behavior  of  tha  Instrument  after  It  Is 
eoaiplatad.  tinea  minor  differences  In  the  parts  used  and  tha  assembly 
procedure  Mill  always  exist*  tha  behavior  of  each  Instrument  will  be 
uni  qua.  However*  there  does  exist  a  general  characteristic  behavior  of 
the  gravity  meter  caused  by  the  non— linearity  of  the  lever  linkage 
system  [Harrison  and  LaCoste*  197S]  which  can  be  Identified.  The 
general  characteristic  behavior  that  Is  sought  Is  how  differences  In 
counter  unit  readings  at  two  sites  are  related  to  the  gravity 
difference  between  the  two  sites.  The  method  that  enables  this 
relationship  to  be  deduced  Is  commonly  referred  to  as  the  calibration 
procedure. 

The  calibration  procedure  Is  a  two  step  process.  The  first  step* 
which  will  be  referred  to  as  the  factory  calibration  procedure* 
determines  the  generel  behavior  of  the  gravity  meter  over  Its  operating 
range  by  determining  what  will  be  called  relatl ve  scale  factors.  The 
second  step*  which  will  be  referred  to  as  the  field  calibration 
procedure*  relates  the  gravity  difference  between  two  stations  to  the 
difference  In  counter  unit  readings  between  the  stations  enabling  what 
will  be  called  the  absolute  scale  fact ors  for  the  gravity  meter  to  be 
determined.  The  end  product  of  the  calibration  procedure  Is  the 
Calibration  Table  i  which  relates  the  gravity  meter's  counter  readings 


to  their  relative  values  In  mllllgals 
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To  determine  the  general  characterise c  behavior  o f  tho  gravity 
■•tor  raquiroa  that  observations  bo  made  over  tho  ontiro  range  of  tho 
counter,  fro*  0  to  9999  counter  units.  One  way  to  accomplish  this 
would  bo  to  have  tho  gravity  motor  make  observations  at  various 
stations  with  known  gravity  values  distributed  over  tho  ontiro  range  of 
tho  gravity  motor.  This  of  course  would  bo  time  consuming,  expensive 
and  impractical.  Another  way  to  achiovo  tho  same  results  would  bo  to 
somehow  fool  tho  gravity  mater  into  thinking  that  tho  gravity  value  at 
a  site  had  changed. 

Since  tho  gravity  motor  works  on  tho  principle  that  tho  force  on  a 
spring  and  tho  resulting  ehango  in  tho  spring's  length  is  related  to 
tho  muss  it  is  supporting  and  the  acceleration  of  gravity  on  that  mass, 
one  way  of  achieving  a  change  in  the  force  without  changing  the 
acceleration  due  to  gravity  would  be  to  vary  the  mass  being  supported, 
to,  if  there  were  a  way  to  change  the  mass  of  the  gravity  meter's  beam, 
then  reedings  over  the  entire  range  of  the  gravity  meter  could  be 
obtained.  Due  to  the  construction  of  the  gravity  meter,  the  gravity 
meter's  beam  acts  as  a  lever.  The  effect  ef  changing  the  mass  of  the 
gravity  meter's  beam  can  be  achieved  by  changing  the  center  of  mass  of 
tho  boom  without  actually  changing  the  mass  ef  the  beam.  An  apparatus 
was  develeped  which  enables  bath  the  center  ef  mass  of  the  beam  to 
change  and  messes  te  be  added  end  remeved  from  the  beam.  The  name 
given  te  this  apparatus  is  Cloudcrof t .  Jr . 


2.5.1  Clouderoft .  Jr.  Apparatus 


Clouderoft,  Jr.  Is  a  device  that  consists  of  two  shafts  that  can  bo 
scrowad  into  tho  beam  weight  and  a  housing  that  is  attachad  to  tha 
meter  box.  Saa  Figura  9  for  a  schematic  view  of  tha  Clouderoft,  Jr. 
apparatus.  Ona  of  tha  shafts  is  threaded.  On  this  threaded  shaft  is 
placed  the  ranaa  adjustment  nut  which  can  be  positioned  anywhere  along 
tha  shaft  by  simply  rotating  that  nut.  Tha  range  adjustment  nut  acts 
as  a  counter-balance  and  by  moving  it  aithar  towards  or  away  from  tha 
beam  waight,  it  can  be  used  to  change  tha  canter  of  mass  of  tha  beam. 
Tha  other  shaft  has  an  opening  at  ona  and  into  which  a  pin  can  bo 
inserted.  On  tha  pertruding  and  df  tha  pin  is  attachad  with  a  drop  of 
glua  a  thin  wire  which  supports  a  mass  which  is  referred  to  as  a  bob. 
Two  thin  metal  wafers  called  weights  rest  on  tha  top  of  tha  bob. 

Tha  housing  that  is  attached  to  the  meter  bow  is  referred  to  as  tha 
bucket .  Zt  enables  the  weights  to  be  lifted  off  and  returned  to  the 
bob  by  a  process  of  raising  and  lowering  the  bucket.  The  bucket  is 
positioned  so  that  when  the  bucket  is  raised,  the  bob  will  descend 
freely  into  the  bucket.  As  tha  bucket  is  raised,  each  weight  comes  to 
rest  on  a  separate  ledge  of  the  bucket,  which  results  in  the  mass  being 
removed  from  the  beam.  Hlth  this  configuration,  either  both  weights 
are  resting  on  the  bob,  the  smaller  weight  is  resting  on  the  bob  or 
neither  weight  is  resting  on  the  bob. 

With  the  weights  being  either  on  or  off  the  bob,  the  terms  of 
weiaht-on  and  weiaht-off  are  used  respectively. 


Figure  9  —  Schoaatic  illustration  of  the  Cloudcroft  Jr.  calibration 
apparatus  used  in  the  calibration  LaCosto  4  Roaborg  "O’* 
gravity  aotor. 


The  wight*  uitd  are  vary  thin  with  vary  little  mi*.  Im  Figure  10 
far  tha  approximate  configuration  of  each  weight.  The  weights  are 
constructed  In  such  way  that  when  the  larger  weight  Is  added  to  the 
beaai  It  causes  a  change  of  approximately  Ut  counter  units.  When  the 
smaller  weight  Is  added,  It  causes  a  change  of  approximately  20  counter 
units,  tinea  a  change  ef  one  counter  unit  Is  approximately  equivalent 
to  a  change  of  one  mgal,  the  larger  weight  Is  referred  to  as  the  1»0 
moal  wal eht  and  the  smaller  weight  Is  referred  to  as  the  ZSL  fiflftl 
wal cht . 

With  this  system,  changes  of  approximately  20,  ItO  and  200  counter 
units  can  be  Induced  In  the  gravity  meter.  And  by  moving  the  range 
adjustment  nut  along  the  threaded  shaft  by  means  of  a  fork  shaped 
device,  the  gravity  meter*a  counter  can  be  made  to  read  any  value 
within  Its  range.  The  Coulderoft,  Jr.  apparatus  Is  Installed  to 
perform  the  factory  calibration  and  than  removod  prior  to  tho  field 
calibration. 

2.B.2  usiaex.  fiiiitecill an  rrmdwri 

The  factory  calibration  procedure  Is  quite  simple  to  perform  but 
somowhat  time  consuming.  The  procedure  normally  takes  from  10  to  12 
hours  to  obtain  the  approximate  120  observations  required  for 
determining  the  relative  scale  factors  over  the  operating  range  of  the 
gravity  meter.  The  steps  Involved  In  this  procedure  are: 

1)  the  Ceuldereft,  Jr.  apparatus  Is  Installed  and  the  meter  Is  put 


on— heat . 


ss 

2 >  t  ho  rang*  adjustment  nut  la  adjusted  in  tha  weight— off  position 
to  road  a  value  near  on*  and  oi  tha  counter  rang*. 

2)  null  the  instrument  with  the  weight— off  and  record  the  reading. 

4)  null  the  Instrument  with  the  weight— on  and  record  the  reading. 

B>  repeat  step  2). 

«>  repeat  step  ♦>. 

7)  adjust  the  range  adjustment  nut  so  the  weight— off  reading  is 

approximately  200  dial  units  larger  or  smaller  depending  on  which 
end  of  the  counter  range  first  set  of  observations  were  mad*. 

2>  repeat  steps  2>  through  7)  until  the  opposite  end  of  the  counter 
range  is  reached. 

1 

From  the  recorded  observations  a  set  of  differences,  with  each 
difference  being  the  average  difference  between  the  weight— on  and 
weight— off  observations  for  a  position  of  the  range  adjustment  nut.  is 
determined.  This  set  of  differences  is  then  divided  by  an  arbitrary 
value,  generally  around  200.  to  produce  a  set  of  relative  scale 
factors.  These  relative  scale  factors  relate  how  the  change  in  gravity 
resulting  from  the  addition  of  a  constant  mass  varies  over  the  range  of 
the  Instrument.  The  set  of  relative  scale  factors  obtained  is  assumed 
to  be  valid  at  the  average  of  the  weight-on  and  weight— off  readings  for 
each  position  of  the  range  adjustment  nut.  The  set  ef  relative  scale 
factors  Is  then  plotted  against  their  average  weight— on  and  weight— off 
readings.  Oenerally.  the  relative  scale  factors  over  the  entire  rang* 
ef  the  instrument  are  net  permitted  to  vary  by  more  than  2  parts  in 
1009.  This  Is  den*  because  the  graph  paper  used  for  plotting  these 
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rtUtivt  seals  factors  as  a  function  of  counter  units  at  tho  seals 
desired  does  not  peralt  the  relative  scale  factor  to  have  any  larger  of 
a  variation  (Perry,  1930,  private  eoaaunlcatlon] .  An  arbitrary  curve 
that  is  supposed  to  represent  the  data  points  Is  drawn  free  hand  or  by 
using  a  curve  template.  This  curve  Is  referred  to  as  the  calibration 
curve,  tee  figure  11  for  an  exaaple  of  the  plotted  data  used  to 
produce  the  calibration  curve.  The  calibration  curve  produced  will  not 
necessarily  go  through  all  the  relative  scale  factor  data  points.  The 
discrepancy  between  the  curve  and  the  data  points  oan  be  easily  as 
large  as  1  part  In  10900  ef  the  relative  scale  factor  value.  If  the 
resulting  calibration  curve  shows  any  unexpected  strange  behavior  such 
as  erratic  dips  or  heaps,  atteapts  are  aade  to  reaove  the  undesired 
behavior  by  changing  parts  of  the  Instruaents  such  as  the  gear  box 
and/or  the  aeasurlng  screw  (Perry,  1980,  private  eoaaunlcatlon] .  The 
last  resort  would  be  to  aedlfy  or  rebuild  the  lever  linkage  asseably. 

If  a  eoaponent  Is  replaced,  such  as  a  gear  box.  It  does  not  aean  that 
the  one  reaoved  Is  bad  and  eannet  be  used  again.  Many  tlaes, 
eoaponents  reaoved  frea  one  Instruaent  whose  the  calibration  curve  was 
not  satisfactory  will  not  produce  any  adverse  effects  In  the 
calibration  curve  when  re-installed  in  another  Instruaent  [Perry,  1980, 
private  eoaaunlcatlon]. 

The  Calibration  Table  1  for  a  gravity  aeter  Is  deteralned  prior  to 
the  delivery  of  thd  gravity  aeter  to  the  eustoaer  and  Is  not  altered 
unless  the  gravity  aeter  Is  returned  to  the  factory  and  a  aajor 
aodlf leatlon,  such  as,  replacing  the  aeasurlng  screw,  gear  box  or  long 


Figure  11  —  A  plot  of  relative  scale  factors  for  the  LeCeste  A  Romberg 
"0"  gravity  meter  "0-220". 
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levar  [Ferry,  lflt(  private  eoMMunloatlen] .  This  would  lapl/  that 
there  axlsts  only  and  Calibration  Table  1  'for  oaeh  gravity  ootar  which 
la  valid  at  any  tine.  Tho  Calibration  Tablo  1  la  tho  roault  of  factory 
calibration  procedures.  Tho  truth  of  tho  oattor  la  that  thla  fa  not 
alwaya  tho  eaao  which  loads  to  a  lot  of  confusion. 

It  appoars  that  tho  Boodotie  Survey  Squadron  out  of  F.  C.  War ran, 

AFS  In  Wyonlng,  which  Is  responsible  for  Halting  tho  Majority  of  tho 
gravity  aotor  obaor vat Ions  In  tho  United  States,  produces  their  own 
Calibration  Tablo  1.  Tho  difference  botwoon  their  Calibration  Tablo  1 
and  tho  one  provided  by  tho  Manufacturer  la  generally  just  a  constant 
scale  factor  applied  to  the  values  In  Ml  111  gal.  A  now  Calibration 
Table  1  Is  produced  periodically  because  It  la  believed  that  the 
gravity  Meter,s  calibration  changes  with  tloe  [Geruff,  19S0,  private 
eoaHaunleatlon]. 

When  the  Geodetic  Survey  Squadron  concludes  that  the  calibration  of 
the  gravity  Motor  has  changed.  It  detorMlnos  tho  scale  factor  that  It 
wishes  to  be  apply  to  the  value  In  Mllllgel  and  often  request  LaCoste  S 
Roaberg  to  produce  a  naw  Calibration  Table  1  for  thaw  using  tho  ssmo 
feroat  as  tho  original  Calibration  Table  1  [Ferry,  19S0,  private 
eoanunleatlon).  In  the  process,  due  to  round-off,  tho  now  Boodotie 
Survey  Squadron(a  Calibration  Table  1  values  In  Ml  111  gal  are  not  an 
exact  scale  factor  Multiple  of  the  original  Calibration  Tablo  1 
supplied  with  tho  gravity  aster.  This  aakes  It  very  difficult  to 
deteralne  which  Calibration  Tablo  l  should  be  used  because  there  no 
rework  the  Calibration  Table  1  to  Indicate  that  the  table  has  been 


Modified.  This  porpotuotos  the  notion  that  tho  coll  brat Ion  of  tha 
gravity  Motors  ehangas  with  tloo. 

Xt  Is  known  that  by  Making  changes  In  tho  levor  linkage  assoMbly, 
tho  general  characteristics  of  the  calibration  curve  can  be  changed. 
The  prlMe  example  of  this  Is  the  gravity  Meter  *G— 233*.  This  gravity 
Meter  was  especially  constructed  so  Its  calibration  curve  was  flat, 
that  Is  all  the  relative  scale  factors  had  the  same  value.  This  Is 
probably  not  exactly  true  but  to  within  2  to  3  parts  In  10000.  the 
relative  scale  factors  are  the  sane.  Given  enough  time  and  funds  any 
"0"  gravity  Meter  could  be  constructed  with  a  flat  calibration  curve 
[LaCoste,  1980,  private  coMMunlcatlon]. 

Whether  "0"  gravity  meters  with  flat  calibration  curves  are  better 
than  those  that  do  not  have  flat  calibration  curves  Is  hard  to  say 
because  only  one  "0"  gravity  Meter  Is  known  to  have  such  a 
characteristic. 

Once  an  acceptable  cell  brat Ion  curve  has  been  obtained,  the 
Instrument  Is  sent  for  Its  field  calibration. 

2.8.3  Lj.il Calibration  rrwitiwrt 

The  purpose  of  the  field  calibration  procedure  Is  to  enable  the 
absolute  scale  factors  to  be  determined.  The  absolute  scale  factors 
relate  the  counter  units  to  their  values  In  mllllgal.  This  Is 
aeeeMplIshed  by  taking  the  Instrument  to  an  area  near  Cloudcroft,  New 
Mexico  where  two  stations  exist,  Cloudcroft  and  La  Lux,  which  have  a 
gravity  difference  of  about  242  mgal.  A  number  of  repeated 
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observations  arc  node  between  those  two  stations.  Prom  those 
observations  on  overage  counter  difference  is  determined.  The  actual 
gravity  difference  between  Cloudcroft  and  La  Luc  is  critical  in  as  ouch 
as  the  better  the  value,  the  closer  the  value  in  milligal  found  in  the 
Calibration  Table  1  will  reflect  true  milligal  units.  This  is 
important  only  if  the  Calibration  Table  I's  values  in  milligal  are  to 
bo  used  without  being  adjusted. 

Using  an  assumed  value  for  the  gravity  difference  between  these  two 
stations,  a  field  scale  factor  is  computed.  It  relates  the  counter 
unit  difference  to  the  value  in  milligal  difference  by  dividing  the 
gravity  difference  by  the  average  counter  difference.  Although  the 
field  scale  factor  determined  in  this  manner  is  truly  only  volid  over 
the  range  of  the  readings  used  in  its  determination,  the  Calibration 
Table  1  is  assumed  to  be  velid  for  the  entire  range  of  the  gravity 
meter. 

*•*•*  Construction  Uu  fiiUJttlilgft  Table  1 

After  the  factory  and  field  calibration  procedures  have  been 
completed,  the  Calibration  Table  1  is  predueod.  tee  Table  2  for  an 
example  of  e  Calibration  Table  1  as  supplied  by  LaCeste  A  temberg,  Inc. 
It  is  very  important  to  understand  how  the  Calibration  Table  1  is 
produced  end  what  type  of  Information  this  table  does  and  does  not 
contain.  This  table  relates  counter  readings  to  volue  in  milligal.  ly 
reading  relative  scale  factor  values  off  of  the  plotted  calibration 
curve  at  Intervals  of  100  counter  units  and  starting  at  >0  counter 
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units,  a  table  of  counter  units  and  relative  scale  factors  Is  produced. 
These  relative  scale  factors  are  assumed  to  be  valid  for  plus/ml nus  50 
counter  units  from  the  point  on  tho  calibration  curve  that  the  reading 
was  made.  These  relative  scale  factors  are  then  all  scaled  by  the 
field  scale  factor  to  produce  what  Is  referred  to  as  the  factor  for 
Interval .  The  factors  for  Interval  are  assumed  to  be  valid  for 
plus/minus  50  counter  units.  From  this  Information,  the  Calibration 
Table  1  is  produced -which  relates  the  counter  readings  to  value  In 
mllllgal  via  the  factors  In  interval. 

In  the  Calibration  Table  1,  the  factor  for  interval  Is  assumed  to  be 
valid  for  a  range  of  80  counter  units  either  side  of  its  corresponding 
counter  reading.  The  value  In  mllllgal  for  a  counter  reading  Is 
obtained  by  multiplying  the  factor  for  interval  by  100  counter  units, 
which  Is  the  difference  between  two  consecutive  counter  readings,  and 
adding  It  to  the  previous  velue  for  the  value  In  milligal.  It  should 
be  noted  that  the  value  In  milligal  Is  derived  from  the  factor  for 
Interval  values  and  not  the  converse.  If  one  assumes  that  the  standard 
error  of  the  observed  difference  of  gravity  motor  readings  between 
Couldcroft  and  La  Luz  is  on  the  order  of  0.025  counter  units,  then  this 
Implies  that  the  field  scale  factor  determined  and  the  corresponding 
factor  for  Interval  ef  the  Calibration  Table  1  Is  accurate  to  about 


1  part  In  10000 
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!•  2 


Calibration  Table  1  for  tho  LaCoato  2  Romborg  "O’*  gravity 

■otor  "0-220". 


TABLE  1 


mUGAL  VALORS  FOR  LAC08TE  4  ROMBERG.  INC.  MOSEL  G  GRAVITY  METER  t G-  220 


comrrER 

VAUJR  IH 

FACTOR  FI 

READING* 

HXLLXGALS 

UREXVAL 

000 

000.00 

1.06106 

100 

106.11 

1.06094 

200 

212.20 

1.06083 

300 

318.28 

1.06074 

400 

424.36 

1.06065 

500 

530.42 

1.06060 

400 

636.48 

1.06057 

700 

742.54 

1.06057 

BOO 

84R.60 

1.06059 

non 

654.66 

1.06063 

1000 

1060.72 

1.06067 

lion 

1166.79 

1.06074 

1200 

1272.86 

1.06080 

1300 

1378.64 

1.06080 

1400 

1485.03 

1.06097 

1500 

1361.12 

1.06104 

1400 

1697.23 

1.06113 

1700 

1803.34 

1.06123 

1000 

1609.46 

1.06128 

1900 

2015.59 

1.06137 

2000 

2121.73 

1.06146 

2100 

2227.88 

1.06156 

2200 

2334.03 

1.06169 

2300 

2440.20 

1.06182 

2400 

2546.38 

1.06197 

2500 

2652.53 

1.06213 

2600 

2758.79 

1.06228 

2700 

2865.02 

1.06242 

2800 

2971.26 

1.06255 

2900 

3077.52 

1.06268 

3000 

3183.79 

1.06278 

3100 

3290.06 

1.06289 

3200 

3396.35 

1.06301 

3300 

3502.65 

1.06314 

3400 

3608.97 

1.06328 

3500 

•  3715.30 

1.06343 

*  Hotel 

U|At*tand  oSsal  on  eon 

COUKTER 

VALDE  IN 

FACTOR  FOR 

READING* 

MILL  I  GALS 

X2ITER.YAL 

3600 

3821.64 

1.06357 

3700 

3928.00 

1.06369 

3800 

4034.36 

1.06381 

3900 

4140.75 

1.06392 

4000 

4247.14 

1.06403 

4100 

4353.54 

1.06412 

4200 

4459.95 

1.06420 

4300 

4564.37 

1.06426 

4400 

4672.80 

1.06428 

4500 

4779.23 

1.06430 

4600 

4885.66 

1.06431 

4700 

4992.09 

1.06432 

4800 

5098.52 

1.06433 

4900 

5204.95 

1.06413 

5000 

5311.19 

1.06431 

5100 

5417.32 

1.06431 

5200 

.5524.25 

1.06430 

5300 

5630.68 

1.06427 

5400 

5737.11 

1.06423 

5500 

•  5843.53 

1.06418 

5600 

5949.95 

1.06412 

5700 

6056.36 

1.06403 

5800 

6162.76 

1.06391 

5900 

6269.15 

1.06376 

6000 

6375.53 

1.06360 

6100 

6481.89 

1.06343 

6200 

6583.23 

1.06324 

6300 

6694.56 

1.06304 

6400 

6800.86 

1.06284 

6500 

6907.14 

1.06261 

6600 

7013.41  . 

1.06239 

6700 

7119.64 

1.06212 

6800 

7225.86 

1.06181 

6900 

7333.04 

1.04146 

7000 

7438.18 

indicates  sFprnxlmtely  0.1  alllljtal. 


10-11-71 

«r  belta-osrrOa  i  .  a«*i« 
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2.C  Ing.trwiwrtu*  Error  Source 

Whenever  tho  "0"  gravity  mat ar a  behavior  differs  from  that  predicted 
by  a  linear  interpolation  Mithin  the  Calibration  Table  1,  there  are  two 
poasible  explanations.  One  explanation  is  that  the  anomalistic 
behavior  is  in  fact  present  but  the  Calibration  Table  1  does  not 
contain  thia  information.  In  this  context*  short  wave  length  refers  to 
mave  length  less  than  approximately  200  counter  units.  This  situation 
occurs  when  the  short  wave  length  behavior  cannot  be  represented  by  the 
long  uave  length  information  present  in  the  Calibration  Table  1.  This 
type  of  systematic  error  could  be  accounted  for  by  additional 
parameters  in  the  mathematical  model.  An  other  explanation  is  that  the 
anomalistic  behavior  is  erratic  and  random  in  nature  and  thus 
Impossible  to  be  modelled.  The  major  error  sources  that  fall  into 
either  of  these  two  categories  are  periodic  screw  effect .  tares  and 

2.«.i  £icigdifi  Isom  iiltei 

Due  to  the  construction  of  the  "0"  gravity  meter,  there  is  a 
possibility  that  an  angular  rotation  of  the  dial  will  not  produce  a 
strictly  linear  motion  of  the  measuring  screw.  The  departure  from  the 
linear  motion  could  be  due  to  periodic  errors  in  the  measuring  screw, 
eccentricity  in  the  measuring  screw  resulting  in  a  wobble  or  the 
non-linearity  of  the  lever  linkage  assembly  [Harrison  and  LaCoste, 
1972].  If  the  periodic  error  were  In  the  measuring  screw  system  and 
could  be  related  to  the  position  of  the  dial  and  the  counter,  then 
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there  Might  be  a  way  of  Modelling  this  effect. 

There  are  tMO  places  in  the  gravity  Meter  in  which  this  type  of 
effect  could  be  introduced.  One  place  is  in  the  gear  box  [KivinieMi, 
1974;  Harrison  and  LaCoste,  19781  and  the  other  is  at  the  point  of 
contact  between  the  Measuring  screw  and  the  lever  linkage. 

The  gear  box  could  introduce  a  periodic  effect  into  the  observations 
due  to  the  eccentricities  in  the  gears  of  the  gear  box  [Kivineioi, 
1974].  The  period  of  this  effect  would  depend  on  which  gear  box  was 
installed  in  the  gravity  Meter.  Xf  the  gravity  Meter  has  an  old  gear 
box  installed  in  it.  then  periods  of  1208.  1208/17,  184/17  and  1 
counter  units  could  be  present.  Xf  the  gravity  Meter  has  a  new  gear 
box  installed  in  it,  then  periods  of  220.  220/8,  22/8  and  1  counter 
units  could  be  present. 

The  other  place  that  a  periodic  effect  could  be  introduced  is  at  the 
point  of  contact  between  the  Measuring  screw  and  the  lever  linkage. 

This  results  when  the  ball  on  the  lever  linkage  and/or  the  hole  in  the 
hardened  jewel  is  net  spherieel  or  circular  in  shape.  Xf  this  were  the 
case,  then  each  rotation  on  the  Measuring  screw  would  produce  a  type  of 
period  effect.  Per  the  old  gear  box,  this  effect  would  occur  every 
1208/17  counter  units,  while  for  the  new  gear  box,  this  effect  would 
occur  every  220/8  counter  units.  Note  that  the  period  of  this  period 
effect  is  a  function  of  the  which  gear  box  is  Installed  in  the  gravity 
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*3.4.2  Hu 

Taro  Is  a  tarn  which  refers  to  unexplained  changes  in  the  reading 
level  of  the  "0”  gravity  eater.  A  tare  in  the  gravity  meter  is 
believed  to  be  the  result  of  small  shifts  of  the  components  lever 
linkage  that  are  screw  clamped  together  {Burris*  1930,  private 
communication].  Tares  by  nature  are  unpredictable  but  are  easily 
Introduced.  A  rapid  deceleration  or  acceleration  of  the  gravity  meter 
is  a  common  cause  that  will  introduce  a  tare.  This  occurs  when  the 
gravity  meter  is  dropped  or  jarred  especially  when  the  beam  is  not 
clamped.  Therefore,  it  is  very  Important  that  the  arrestment  knob  be 
turned  fully  clockwise,  so  that  the  beam  is  clamped  whenever  the 
gravity  meter  is  being  moved. 

Large  tares,  on  the  order  of  ISO  ygal  or  larger,  are  generally  easy 
to  detect.  But  smaller  tares  can  be  very  difficult  to  identify.  Any 
gravity  meter  tie  suspected  of  containing  a  large  tare  should  be 
removed  from  the  observation  set.  But  there  is  little  that  can  be  done 
for  the  ties  that  contain  the  undetected  small  tares. 

2.4.3  Instrumental  Drift 

The  drift  of  the  "0"  gravity  meter  is  not  totally  understood  at  this 
time.  There  appears  to  be  no  mechenieal  reason  why  readings  made  with 
a  properly  adjusted  "0"  gravity  meter  should  change  with  time  other 
than  as  a  result  of  tares  being  Introduced  [ferry,  1930,  private 
communication].  It  is  believed  that  the  so  called  instrumental  drift 
is  the  cumulative  result  of  a  number  of  small  tares  in  the  gravity 
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Mt*r  [Uotlla,  1974i  Burris,  1980,  private  communication].  Tht  tarsa 
occur  randoaily  rather  than  unlfonaly  which  aakos  Modelling  of  such  an 
offset  vary  difficult.  If  not  laposalblo. 


CHAPTER  THREE 


g&AYIT.T  flEIig  osservatiohs 

S.l  Oburvttiflni  Used 

For  this  study,  observations  wore  obtained  from  two  governmental 
organisations:  the  Hatlonal  Geodetic  Survey  of  the  Hational  Ocean 
Survey  of  the  Hatlonal  Oceanic  and  Atmospheric  Administration  of  the 
United  States  Department  of  Commerce  and  the  Beodetlc  Survey  Squadron 
stationed  at  F.f.  Warren  APS,  Wyoming  of  the  Defense  Happing  Agency 
Hydrographl e/Topegraphl c  Center.  The  data  obtained  consist  of  over 
4500  gravity  meter  observations  made  with  25  different  *8*  gravity 
meters  and  2  different  TDT  gravity  meters.  The  majority  of  the 
observations  were  made  along  the  United  States  Mid-Continent 
Calibration  Line  which  runs  along  the  Eastern  side  of  the  Rocky 
Mountains  with  stations  In  Texas,  Hew  Mexico,  Colorado,  Wyoming, 

Montana  and  North  Dakota.  All  observations  were  made  during  the  period 
1975— 19S0  and  were  received  In  the  form  of  copies  of  the  original 
observation  sheets.  See  Table  5  for  a  listing  of  the  gravity  meters 
used.  See  Figure  12  for  the  geographical  location  of  the  stations  In 
the  network  and  how  they  are  Interconnected. 

The  Information  supplied  en  the  observation  sheets  consisted  of  tho 
name  of  the  observer,  the  Instrument(s)  being  used,  the  station  name. 
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3  —  Summary  of 

the  gravity  meters  used 

In  study. 

Gravity 

Meter 

Date  on 
Calibration 

Number  of 

Observations  Made 

Number 

Table  1 

Observations 

From 

To 

0-10 

25/10/60 

164 

03/05/79 

29/10/50 

0-44 

25/04/63 

14 

10/02/75 

16/02/75 

0-47 

22/05/63 

15 

10/02/75 

16/02/75 

0-50 

15/06/63 

17 

11/05/75 

11/05/75 

0— 45 

23/ OS/64 

61 

05/10/76 

26/10/76 

0-51 

07/05/64 

427 

17/04/75 

26/10/76 

0— 51a 

17/10/77 

171 

11/01/75 

30/10/50 

0-1  OS 

07/09/65 

125 

11/01/75 

02/10/79 

0-111 

25/03/66 

427 

17/04/75 

26/10/74 

0-113 

25/03/64 

17 

11/05/75 

11/05/75 

0-1  IS 

09/05/46 

369 

17/04/75 

10/11/75 

0-1 15b 

14/02/75 

300 

11/01/75 

09/02/50 

0-123 

22/10/79 

39 

15/09/79 

02/1 ?/79 

0-125 

17/10/64 

140 

25/04/75 

05/02/7? 

0-130 

15/10/44 

19 

03/05/79 

03/05/79 

0-131 

15/05/75 

350 

27/03/75 

16/05/50 

0-140 

24/02/67 

17 

12/05/75 

12/05/75 

0-142 

14/03/67 

77 

07/01/75 

02/10/79 

0-157 

10/05/67 

435 

17/04/75 

26/10/76 

O— 157c 

25/01/75 

172 

11/01/75 

30/10/50 

0-175 

30/04/65 

17 

12/05/75 

12/05/75 

0-176 

19/04/65 

27 

09/02/75 

17/02/75 

0-191 

27/01/69 

16 

07/01/75 

25/04/75 

0-220 

11/10/75 

266 

13/05/75 

09/02/50 

0-253 

09/10/75 

114 

27/03/75 

15/11/79 

0-265 

15/05/75 

237 

27/03/75 

09/02/50 

0-269 

29/04/71 

147 

09/02/75 

02/10/79 

0-17 

» 

152 

19/20/77 

23/06/50 

0-43 

m 

126 

12/05/50 

23/06/50 

All  ditu  art  given  In  day.  Month,  year  order. 

•  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  IS  October  1977. 

b  —  Calibration  Table  1  changed  due  to  replaeoMent  of  long  lever  on 
27  October  1977. 

e  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  SO  Auguat  1977. 

«  —  No  Calibration  Table  1  la  provided  Mlth  *D*  gravity  Metero  since 
the  scale  factor  Is  aaauMed  to  be  a  constant  [LaCoste  4  Romberg, 
1979b]. 
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iMttim  •  station  Identl  float  Ion  number,  data  and  time  of  each 
observation  to  tha  minute,  obaarvod  dial  reading  In  dial  units,  height 
of  gravity  meter  above  or  below  the  station,  station  coordinates  and 
any  remarks  concerning  unusual  operating  conditions  or  Instrument 
behavior.  The  time  was  given  either  In  Universal  Coordinated  Time  or 
In  local  standard  time  with  the  correction  needed  to  obtain  Universal 
Coordinated  Time.  The  coordinates  of  the  station  were  given  as 
latitude,  longitude  and  elevation  with  the  latitude  and  longitude 
generally  given  to  the  nearest  0.1  minutes  and  the  elevation  given  to 
0.01  meters  or  equivalent.  See  Figure  13  and  Figure  14  for  sample  of 
observation  sheets. 

3.2  Observational  Procedure 

The  observational  procedure  recommended  Is  outlined  In  Land  Gravity 
Surveys,  DMAHTC/OSS— TM-9,  Preliminary  Edition,  October  1979  on  page  4—1 
which  basically  states  that 

1>A  valid  set  of  observations  consists  of  those  made  by  one 
observer.  This  Is  necessary  to  eliminate  parallax  and  other 
observer  peculiarities.  The  gravity  meter  must  have  been  at 
operating  temperature  for  at  least  4  hours  prior  to  beginning 
observations  and  during  the  observations  the  operating 
temperature  must  be  maintained. 

2)  The  gravity  meter  may  be  placed  directly  on  any  smooth,  hard, 
level  surface  for  observing.  If  any  of  these  conditions  are  not 
met,  then  the  gravity  meter  should  bo  placed  on  the  levelling 
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dish  which  Must  bs  firstly  seated  to  eliminate  any  movement  while 
the  observation  is  being  made. 

S>  The  gravity  meter  must  be  levelled  and  then  rough  nulled.  Rough 
nulling  is  accomplished  by  turning  the  arrestment  knob 
counter-clockwise  which  unelamps  the  balance  beam  and  then 
rotating  the  dial  until  the  beam  is  off  its  stops  but  not 
necessarily  to  the  reading  line.  The  rough  nulling  condition 
should  exist  for  approximately  5  minutes  before  an  observation  is 


made.  During  this  time,  the  observer  will  enter  station 
description  information.  The  observer  must  keep  the  sun  from 


shining  on  the  gravimeter  because  the  heat  might  cause  distortion 
in  the  level  vial  assembly.  Mhen  finished  with  the  observation, 
the  gravity  meter*s  balance  beam  must  be  clamped  by  turning  the 
arrestment  knob  clockwise  and  the  gravity  meter  returned  to  its 
carrying  case  with  the  carrying  ease  lid  closed  to  prevent  the 

gravity  meter  from  being  tipped  over  by  the  wind. 

/ 

4)  The  gravity  meter  is  nulled  by  approaching  the  reading  (nulling) 
line  from  the  down— scale  (left)  side  to  the  up-scale  (right) 
side.  The  null  position  is  the  coincidence  of  the  left  edge  of 
the  cross— heir  with  the  reeding  line.  If  the  observer  overshoots 
the  reeding  line,  the  dial  must  be  offset  ISO  degrees  down— scale 
and  the  reading  line  approached  again.  This  must  be  done  to 
eliminate  any  backlash  in  the  dial  gear  system. 

5)  A  valid  observation  at  a  station  consists  of  two  consecutive 
nulling,  ne  mere  than  4  minutes  apart,  that  agree  to  0.01  counter 
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units. 

Of  ths  sens  4800  gravity  astsr  observations  received,  how  many  of 
those  observations  were  obtained  following  the  procedure  outlined  above 
Is  unknown.  But  what  Is  known  Is  that  sons  of  the  observers  did  not 
follow  that  procedure.  Sees  of  the  data  received  consisted  of  a  single 
null  reading  per  observation.  This,  In  Itself.  Is  net  necessarily  bad. 
But,  the  reason  for  the  two  consecutive  nullings  Is  to  provide  a  Method 
of  detection  of  blunders  In  nulling,  reading,  and/or  recording. 

Without  the  second  nulling,  the  detection  of  these  types  of  blunders 
becomes  Impossible.  In  addition,  no  Information  Is  available 
concerning  the  repeatability  of  observations  made  with  the  Instrument. 

Another  practice,  known  to  occur  but  not  how  often.  Is  that  of  not 
actually  performing  the  second  nulling  [Beruff,  1951,  private 
communication].  Instead,  a  type  of  quasi— nulling  Is  performed.  After 
the  first  nulling  has  been  performed,  the  second  nulling  consisting  of 
the  dial  being  backed  off  at  least  the  required  159  degrees  and  then 
the  dial  being  set  back  at  the  first  nulling  position.  The  cross— hair 
Is.  checked  and  If  acceptable,  the  second  nulling  recorded  is  made 
identical  to  the  first  reading.  What  information  this  type  of 
procedure  provides.  If  any,  is  not  clear.  But  this  type  of  practice  is 
not  recommended  and  should  not  taka  plaoe. 

Another  practice  which  is  not  uncommon  Is  the  inconsistent  recording 
of  the  height  of  Instrument  [Wassails,  1950,  private  communication). 
This  occurs  when  for  some  visits  to  a  station,  the  levelling  disk  is 
used,  while  for  other  visits  to  the  same  station,  the  levelling  disk  Is 
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wot  used  and  no  height  of  Instrument  is  recorded.  The  levelling  disk 
can  elevate  the  gravity  meter  as  much  as  5  cm  which,  if  the  gradient  of 
gravity  is  assumed  to  be  that  of  normal  gravity,  would  result  is  a 
systematic  gravity  change  of  about  15  ygal.  For  this  reason,  the 
levelling  disk  should  always  be  used  when  making  gravity  meter 
observations  or  the  height  of  instrument  of  the  gravity  meter  should  be 
properly  recorded. 

3.3  Korkina  with  Gravi tv  Meter  Data 

All  gravity  meter  data  used  was  received  as  copies  of  the  original 
field  observation  sheets.  The  vital  information  necessary  to  perform 
the  gravity  base  station  network  adjustment  was  extracted  and  encoded 
for  use  in  the  computer. 

The  Information  encoded  was  the  station  information,  observation 
information,  and  instrument  information.  The  station  information 
included  the  station's  identification  code,  the  station's  name,  and  its 
location  given  by  its  latitude  and  longitude  to  0.1  minutes  and  its 
elevation  in  meters.  The  observation  information  included  the  recorded 
time  of  the  each  nulling  to  the  minute  and  its  correspondi ng  observed 
counter  reading.  The  time  recorded  was  either  in  Universal  Coordinated 
Time  (UTC)  or  in  local  standard  time  with  the  correction  needed  to 
obtain  UTC.  The  instrument  information  consisted  of  t;ie  identification 
number  of  the  gravity  meter  used  in  making  the  observations  along  with 
an  arbitrary  date  set  number  which  was  assigned  to  each  set  of 


observation  sheets  as  they  were  received. 


The  encoded  data  ware  than  visually  checked  -for  agreement  with  the 
iginal  data  and  any  errors  corrected.  Each  observation,  which 
lerally  consisted  of  two  separate  nullings,  is  the  average  of  their 
11  readings.  Ideally,  each  nulling  is  made  independently.  The 
rpose  of  the  two  independent  nullings  is  to  check  for  blunders  in  the 
Lling  process  and  in  the  recording  of  the  reading.  The  time  of  each 
tervation  is  the  average  of  the  time  for  the  nullings.  The  resulting 
Lues  were  used  as  the  observation  and  the  time  of  the  observation  in 
>  adjustment  program.  In  the  process,  the  time  which  was  originally 
/en  in  year,  month,  day,  hour  and  minute  format  was  converted  into  a 
re  convenient  form,  its  Julian  Date. 

The  way  the  time  of  the  nullings  is  actually  recorded  should  be  made 
iformly.  Either  all  times  are  given  in  local  standard  time  with  the 
rrection  needed  to  obtain  UTC,  or  they  are  given  in  UTC.  The  problem 
tin  recording  the  time  of  the  nulling  in  local  standard  time  is  that 
»  UTC  correction  varies  with  the  location  of  the  station  and  the  time 
the  year  that  the  nulling  is  made  in  the  United  States.  It  is 
sommended  that  each  gravity  meter  have  a  small  electronic  24  hour 
jital  display  clock  which  also  displays  the  current  date  a.fixed  to 
which  would  be  set  to  UTC.  Then  the  time  of  the  observation  could 
>ily  be  recorded  in  UTC  without  worrying  about  time  cone  or  seasonal 
inges  in  the  local  standard  time. 

•ut  by  far  the  most  confusing,  yet  very  important.  Information  is 
r  station  information  itself.  Cenerally,  complete  station 
formation  is  not  provided  on  the  field  observation  sf-eots.  Instead, 
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just  enough  station  Information  Is  provided  to  Identify  the  station  so 
Its  aravl tv  stat 1  on  descrl pt 1  on  or  site  description  can  be  located.  On 
the  gravity  station  description  form  Is  the  detailed  Information  about 
the  exact  location  of  the  gravity  station  which  gives  Its  latitude, 
longitude  and  elevation  with  a  word  description  of  its  location  plus  a 
diagram/photograph  of  the  station's  location.  See  Figure  15  and 
Figure  15  for  examples  of  gravity  station  description  forms. 

The  confusion  develops  when  the  information  on  the  field  observation 
sheet  is  not  sufficient  to  locate  Its  gravity  station  description,  if 
it  exists,  or  when  the  Information  of-  the  field  observation  sheets  does 
not  completely  agree  with  the  gravity  station  description  information. 
The  latter  problem  occurs  most  often  when  the  latitude,  longitude  and 
elevation  information  is  not  given  on  the  field  observation  sheet.  The 
question  is  not  which  latitude,  longitude  and  elevation  Information 
should  be  used  (that  is  clear;  the  gravity  station  description 
information  should  be  used),  but  rather  where  the  information  on  the 
field  observation  sheet  came  from  and  whether  the  station  is  really  the 
station  it  is  purported  to  be.  To  muddle  the -si tuati on  even  further, 
gravity  station  description  forms  for  the  same  station  have  been 
received  which  are  identical  in  description  and  date  except  for  a 
change  in  a  coordinate  of  the  station  and/or  station  designation.  See 
Figure  15  and  Figure  15  for  examples  of  this  situation. 

To  add  even  more  confusion,  some  stations  do  not  have  a  gravity 
station  description  form.  The  most  common  station  of  this  type  is 
commonly  referred  to  as  a  drift  station  [Spita,  1951,  private 


Figure  15  —  Gravity  station  description  -form  for  station  "Great  Falls 
0"  —  sample  no.  1. 
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communication] .  A  drift  station  Is  generally  established  when  the 
observer  Is  travelling  between  assigned  gravity  station  and  observer 
has  to  make  an  overnight  stop.  The  observer  In  this  situation  will 
establish  the  drift  station  at  this  overnight  stop,  making  an 
observation  that  night  and  again  the  following  morning,  before  starting 
off  for  the  assigned  gravity  station.  The  approximate  coordinates  of 
this  drift  station  can  bo  obtained  from  a  topographic  map  of  the  region 
but  generally  no  gravity  station  description  for  the  station  Is  ever 
produced  since  the  chances  are  the  station  will  never  be  reoeeupled. 

Xt  Is  recommended  that  for  every  station  where  a  gravity  meter 

\ 

observation  Is  made,  there  should  always  exist  a  gravity  station 
description  so  the  station  could  be  reoccupied. 

Another  problem  Is  that  of  the  station  designation  which  can  be  very 
misleading.  A  station  designation  of  just  a  name  Is  not  generally 
enough.  For  example,  there  are  two  stations  In  the  United  States  base 
Station  Network  by  the  name  of  Las  Vaoas  ft.  One  station  Is  In  Nevada 
and  the  other  Is  In  New  Mexico.  Generally,  In  addition  to  a  name  for  a 
station,  an  Identlf Iclatlon  code  such  as  an  International  Gravity 
Bureau  (XGB)  number  Is  associated  with  the  station.  Xn  the  ease  of  the 
station.  Las  Vegas  B,  In  Nevada,  the  identification  code  assigned  was 
120<SB  because  this  station  was  in  the  XGSN  71.  However,  the  Las  Vegas 
B  station  In  New  Mexico  had  no  identification  code  on  Its  gravity 
station  description  sheet.  But  on  some  of  the  field  observation 
shoots,  the  station  appears  with  an  Identification  code  of  11BBSB.  The 
matching  of  gravity  station  description  information  with  the 
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information  supplied  on  the  field  observation  sheets  can  become  very 
confusing.  This  leads  to  the  real  possibility  that  the  same  station 
could  be  Identified  In  the  adjustment  as  tMO  or  more  different 
stations.  This  possibility  would  result  In  a  weaker  network  adjustment 
and  some  confusion  since  there  would  be  two  adjusted  gravity  values  for 
the  same  station. 

s.s.i  aiaiiaa  l^aaliilsalLga 

In  order  to  uniquely  Identify  stations  for  the  Z6SN  71,  the  165, 
International  Gravity  tureau,  number  or  code  was  established  [Horelll, 
et  al,  1974).  The  main  feature  of  this  coding  system  was  that  It 
conveyed  Information  about  the  geographical  coordinates  of  the  station 
It  was  Identifying.  The  10ft  code  consists  of  five  digits  and  a  letter. 
The  first  three  digits  of  the  code  are  determined  from  the  geographical 
coordinates  of  a  station  using  the  following  relationships  given  In 
XGSN  71  [Morel 11 ,  et  al.,  1974].  5ee  Figure  17  for  how  the  three  digit 
code  Is  distributed  over  the  earth. 

The  other  two  digits  of  the  ZGB  code  are  the  units  of  the  latitude 
and  longitude  degrees  respectively.  In  the  formation  of  these  digits, 
no  rounding-off  Is  done. 

In  order  to  Identify  uniquely  stations  that  have  the  same  five  digit 
number,  a  -unique  letter  Is  attached  to  the  end  of  the  five  digit  code. 
This  permits  up  to  24  stations,  one  for  each  letter  of  the  alphabet,  to 
be  assigned  a  unique  code  for  every  lxl  degree  block  on  the  earth. 
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In  a  world-wide  gravity  basa  station  network  Ilka  ICSN  71,  26 
stations  par  lxl  dagraa  block  were  sufficient.  But,  whan  local  basa 
station  natMorks  ara  balng  astabllshad  or  densified,  the  number  of 
stations  In  a  lxl  degree  block  might  exceed  the  26  stations  permitted 
when  using  the  XGB  coda  system.  In  the  United  States,  the  U.  S. 
Interagency  Gravity  Standards  Committee  recommended  a  solution  to  this 
problem.  The  solution  mas  to  add  an  additional  letter  to  the  end  of 
the  XGB  code  which  would  permit  676  or  more  stations  to  be  uniquely 
Identified  in  any  lxl  degree  block.  The  exact  number  depended  on 
whether  blanks  and/or  numerals  were  considered  letters  [Uotila,  19S1, 
private  coaimunleatlon] . 

From  a  data  management  point  of  view,  this  change  might  necessitate 
a  modification  in  the  data  base  structure  for  gravity  meter 
observations  if  the  XGB  code  were  being  used  to  Identify  the  stations. 
This  Is  because  the  proposed  "modified”  XOB  code  would  require  seven 
characters  as  opposed  to  six  characters  presently  being  used. 

'  Xn  order  to  avoid  Increasing  of  the  number  of  characters  needed  to 
Identify  a  station  and  Its  location,  a  possible  solution  would  be  to 
have  kept  a  six  character  code  with  the  first  two  characters 
representing  the  10x10  degree  blocks  instead  of  the  first  three 
characters  as  Is  done  In  the  XGB  number.  This  could  easily  be 
accomplished  since  there  are  26  letters,  a— z,  and  10  numerals,  0—9,  for 
a  total  of  SB  characters  that  can  be  used.  Mlth  thirty— six  10  degree 
Intervals  in  longitude  and  only  eighteen  10  degree  intervals  in 
latitude,  each  10  degree  Interval  In  latitude  or  longitude  can  be 
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r«pra»ant«d  by  on*  of  the  36  character*  available.  This  would  pernlt 
the  geographical  lxl  degree  Information  to  be  represented  In  four 
characters  Instead  of  five  required  by  the  present  XQt  code.  This 
would  leave  a  set  of  two  characters  for  station  Identification  within 
each  lxl  degree  block.  This  coding  method  would  not  change  the  length 
of  the  code  and  relate  the  same  Information. 

However,  If  this  ceding  system  were  used,  an  existing  station  would 
have  a  new  Identification  cede  along  with  Its  181  cod*.  This  could 
lead  to  some  additional  confusion.  But,  If  the  fifth  character  of  the 
new  cod*  were  always  an  alphabetic  character,  then  the  two  codes  would 
be  unique  and  easily  distinguishable  because  the  fifth  character  of  the 
161  cod*  Is  always  numeric. 

The  real  concern  Is  not  which  coding  system  is  adopted,  but  whether 
the  system  will  be  used  by  all  organizations  that  collect  and 
distribute  gravity  base  station  information.  If  a  uniform  coding 
system  is  not  used,  the  confusion  that  can  occur  with  station 
identification  will  persist. 

3.3.2  gravity  Hater  Loops  and  Tries 

A  loom  Isa  set  of  gravity  meter  observations  that  starts  at  a 


station  and  generally  ends  on  the  same  station  after  a  number  of  other 
observations  have  been  made  at  other  stations.  A  loop  can  require  a 
few  hours  to  several  days. to  complete.  The  recommended  types  of  loops 


as  given  In  Land  Oravlty  Surveys  [DHAHTC/OSS,  1979]  are  referred  to  as 
ladder,  modified  ladder,  and  ] Ina  sequence.  A  ladder  and  modified 
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In  the  ladder 


sequence,  every  station  is  observed  as  If  it  were  a  rung  on  a  ladder 
and  the  observer  el labs  up  end  down  the  ladder  stopping  at  every  rung 
to  make  an  observation.  In  the  Modified  ladder  sequence,  rungs  are 
skipped  on  tho  way  up  or  down  the  ladder.  The  modified  ladder  sequence 
loop  is  used  when  difficult  field  conditions  are  encountered. 

Difficult  field  conditions,  for  example,  occur  when  the  station  to  be 
occupied  is  in  a  building  which  is  inaccessible  to  the  observer  on  the 
day  and/or  time  the  observer  tries  to  occupy  the  station.  The  line 
sequence  loop  does  not  end  at  the  starting  station  but  at  some  other 
stations.  See  Figure  IS  for  examples  of  the  types  of  loops. 

A  loop  in  which  the  first  and  last  station  observed  is  the  same  is 
often  referred  to  as  a  closed  loop.  A  closed  loop  can  be  used  to 
determine  if  a  linear  drift  exists  within  an  Instrument  [DMAHTC/GSS, 
1979].  It  is  based  on  the  premise  that  any  difference  between  repeated 
observations  at  the  same  station  is  a  result  of  a  linear  instrumental 
drift.  This,  of  course,  is  not  the  only  reason  why  a  difference 
between  repeated  observations  at  the  same  station  might  exist.  One 
possible  reason  is  that  one  or  mere  tares  could  have  occurred  between 
the  times  of  the  repeated  observations  which  could  result  in  a 
difference  in  the  observations.  Another  possible  reason  is  that  tho 
difference  is  due  to  observational  error  and  not  a  linear  drift. 

Determining  the  presence  of  a  linear  drift  in  an  instrument  does  not 
require  observations  to  be  made  following  a  closed  loop  structure  such 
as  the  ladder  or  the  modified  ladder  sequence.  The  line  sequence  loop 
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Figure  IS  —  Types  of  gravity  aster  loops  possible 


will  work  just  as  wall.  What  Is  nacassary  Is  that  tha  relationship  for 
tha  drift  of  an  instrument  should  bo  included  in  tha  mathematical  model 
usad  to  dascriba  tha  instrument's  behavior.  Tha  mathematical  model  for 
tha  instrument's  behavior  depends  on  tha  two  observations  usad  in 
forming  tha  aquation  and  not  on  what  type  of  loop  the  observations  were 
made. 

A  trio  will  be  defined  as  a  sat  of  consecutive  observations  for 
Mhich  tha  gravity  mater's  behavior  Is  assumed  to  remain  the  same.  A 
trip  could  involve  a  number  of  loops  or  just  part  of  a  loop  depending 
on  hoM  the  instrument  Is  used  and  Mhat  has  happened  to  it.  The 
question  arises  as  to  Mhen  a  trip  begins  and  Mhen  it  ends.  The  answer 
Involves  the  determination  of  Mhen  changes  in  the  behavior  of  the 
gravity  meter  can  be  expected.  There  are  two  good  examples  of  when  the 
behavior  of  the  gravity  meter  might  change.  First  and  the  most  common 
is  Mhen  a  tare  is  introducted  into  the  gravity  meter.  A  tare  reflects 
a 'discontinuity  in  the  behevior  of  a  gravity  meter  and  cannot  be 
modelled  into  the  observati onal  difference  between  two  consecutive 
observations  because  its  magnitude  and  direction  are  unknown.  A  number 
of  conditions  can  result  in  tares.  For  example,  tares  can  occur  when 
the  gravity  meter  is  taken  off-heat  and  put  badk  on— heat,  when  the 
gravity  meter  balance  beam  is  not  clamped  during  transported  on,  when 
the  gravity  meter  experiences  a  rapid  acceleration  and/or  deceleration 
such  as  when  the  meter  is  accidentally  jarred  even  if  the  balancr  beam 
is  properly  clamped  and  even  when  the  measuring  screw  is  being 
lubricated  [Ferry,  1980,  private  communication]. 
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The  ether  condition  under  which  a  new  trip  should  be  started  is  when 
the  activity  of  the  gravity  meter  is  unknown.  If  the  activity  of  the 
instrument  is  not  known  for  a  period  of  time,  too  many  possible  things 
could  have  occurred  to  make  the  gravity  meter's  behavior  inconsistent 
during  that  period  of  time.  For  example,  the  meter  could  have  been 
taken  off— heat. 

To  eliminate  personal  bias  from  being  introduced  when  the 
observations  are  being  divided  into  trips,  a  set  of  rules  was  used  to 
define  what  consitutes  a  trip  for  a  gravity  meter.  The  rules  followed 
were: 

1)  A  new  trip  will  begin  with  the  first  observation  after  a  tare  has 
been  detected. 

2)  A  new  trip  will  begin  when  the  time  between  consecutive 
observations  exceeds  T1  hours. 

3)  A  new  trip  will  begin  when  the  time  between  consecutive 
observations  at  the  same  station  exceeds  T2  hours. 

The  value  of  24  hours  was  assigned  to  T1  and  a  value  of  4  hours  was 
assigned  to  T2.  The  reason  for  choosing  the  value  for  T1  is  that  if  a 
gravity  meter  was  being  used  to  make  observations  in  a  loop,  one  would 
net  expect  the  time  between  observations  to  be  longer  than  a  day  before 
the  loops  was  completed.  The  choice  for  T2  is  based  on  the  premise 
that  repeated  observations  at  the  same  station  would  only  occur  when 
the  observations  made  were  part  of  the  same  loop.  This  occurs, 
generally,  at  overnight  stops.  Then  it  would  be  expected  that  the 
repeated  observations  occurred  the  night  before  and  the  next  morning. 
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which  generally  Meant  a  tine  interval  of  6  hours  or  more.  Therefore, 
repeated  observation  at  the  same  station  separated  by  more  than  6  hours 
would  indicate  that  the  latter  observation  would  be  the  first 
observation  of  a  new  loop. 

To  verify  this  reasoning,  an  extensive  study  was  performed  to  see 
what  the  effect  would  be  of  dividing  the  observations  into  trips  based 
on  various  length  of  time.  The  results  of  this  study  indicated  that 
when  the  value  of  T1  was  increased,  the  apparent  accuracy  of  the 
observed  counter  readings  decreased.  This  was  due  to  tares  being 
included  which  wore  present  in  the  increased  time  interval  previously 
ignored.  Conversely,  as  the  value  of  T1  decreased,  the  accuracy  of  the 
observed  counter  readings  in  some  cases  increased  to  unrealistic 
accuracies  of  less  than  10  ygal  when  T1  was  set  to  4  hours.  This 
occurred  because  generally  only  small  gravity  differences  could  be 
observed  in  that  time  Interval. 

The  study  also  looked  at  the  possibility  of  the  gravity  meter  having 
a  different  behavior  during  the  overnight  stops  as  opposed  to  the 
normal  observation  sequence.  However,  due  to  the  limited  number  of 
overnight  differences  available  for  each  gravity  meter,  there  was  not 
sufficient  evidences  that  any  change  in  the  gravity  meter's  behavior 
occurred. 

S.4  Observational  Errors 

As  mentioned  previously,  errors  such  as  tares  can  easily  be 
introducted  into  gravity  meter  observations.  There  is  little  that  the 


erver  can  do  to  guard  against  these  typo  of  •rrors,  except  handle 
i  gravity  aster  with  care  when  it  is  being  transported.  However, 
is  errors  can  be  Introduced  into  the  observations  no  matter  how  much 
e  is  taken  in  its  transportation.  These  errors  are  due  to  improper 
>ration  and  adjustment  of  the  gravity  meter  itself. 

The  gravity  meter  should  not  be  without  power  for  any  length  of  time 
:ause  without  the  power  the  operating  temperature  of  the  gravity 
er  cannot  be  maintained.  To  insure  a  stable  operating  environment, 

»  gravity  meter  must  be  kept  on— heat.  Therefore,  the  gravity  meter 
>uld  always  be  connected  to  either  a  battery  or  the 

irger— eliminator  except  for  the  short  time  required  to  change  between 
i  two  power  supplies. 

As  mentioned  previously,  transporting  the  gravity  meter  should  only 
ur  when  the  beam  is  clamped.  Failure  to  follow  this  simple  rule 
1,  almost  surely,  introduce  tares  into  the  gravity  meter. 

The  proper  adjustment  of  the  gravity  meter  is  also  very  important  to 
ure  the  stable  behavior  of  the  meter.  The  two  gravity  meter  levels, 
lone  level .  (parallel  to  the  counter),  and  the  cross  level . 
rpendicular  to  the  counter)  must  be  adjusted  according  to  the 
rating  manual  provided  with  each  gravity  meter  [LaCoste  A  Romberg, 

0 ] .  The  quality  of  the  observations  made  with  a  gravity  meter 
ends  to  a  large  part  on  these  two  levels  being  in  proper  adjustment, 
deviation  from  the  correct  position  of  the  levels  will  change  the 
/ity  meter's  sensitivity  and  reading  line. 
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The  sensivitiy  should  bo  checked  periodically  to  assure  that  the 
recommended  sensitivity  is  kept  at  6  to  12  eyepiece  divisions  -for  every 
dial  revolution  [LaCoste  4  Romberg,  I960].  If  the  sensitivity  is  not 
kept  within  this  range,  the  reading  line  for  the  instrument  will 
change.  For  a  decrease  in  the  sensitivity,  the  reading  line  will  shift 
up  the  scale.  Conversely,  for  an  increase  in  the  sensitivity,  the 
reading  line  will  shift  down  the  scale  [DMAHTC/GSS,  1979]. 

Even  if  the  instrument  is  in  perfect  adjustment,  observational 
errors  can  be  introduced  by  the  observer  through  the  nulling  of  the 
instrument  and  the  reading  and  recording  of  the  observation. 

3.5  Honkasalo  Correction  Term 

The  so  called  Honkasalo  correction  term  \s  a  latitude  dependent 
correction  which  was  applied  to  all  absolute  gravity  sites  used  in  the 
I6SN  71  adjustment  [Morelli ,  et  al.,  1974].  The  result  was  that  all 
adjusted  gravity  station  values  published  for  the  IGSN  71  included  the 
Honkasalo  correction  term.  The  correction  is  based  on  the  premi se  that 
the  earth  tide  correction  applied  to  measured  gravity  given  by  equation 
(1)  in  Honkasalo  [1964]  is  only  zero  when  summed  over  the  whole  earth's 
surface  and  not  zero  when  summed  over  a  particular  latitude.  This 
systematic  effect  according  to  Honkasalo  [1964]  should  be  removed  if 
the  earth  tide  correction  summed  over  a  particular  latitude  is  to  be 
zero.  The  amount  to  be  removed  is  given  by  equation  (5)  in  Honkasalo 


[1964]. 
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At  the  XVII  General  Assembly  of  tho  International  Union  of  Geodesy 
an  Geophysics  at  Canberra.  2—14  December,  1979,  Resolution  No.  15  was 
passed  by  the  International  Association  of  Geodesy  which  resolved  that 
the  Honkasalo  correction  should  not  be  applied  to  observed  gravity 
[Uotila,  1950].  Therefore,  all  published  IGSN  71  station  values  should 
have  the  Honkasalo  correction  removed  according  to  the  method  given  by 
Uotila  [1950]. 

In  this  regard,  care  must  be  taken  to  assure  that  any  station  value 
used  as  control  in  an  adjustment  does  not  include  the  Honkasalo  term. 
The  Italian  absolute  station  determinations  as  given  in  Marson  and 
Alaaia  [1975]  include  the  Honkasalo  correction  term.  In  order  to  have 
a  consistent  set  of  absolute  stations,  tho  Honkasalo  term  must  be 
removed  from  those  absolute  station  values  since  the  later  Italian 
determinations  [Marson  and  Alasia,  1950]  and  all  U.  S.  determinations 
made  by  Hammond  do  not  Include  the  Honkasalo  terir  [Hammond,  letter  to 
Uotila,  1951]. 

3-4  t ihx  Ihs.  gfflttltY  ttllwt  Changes 

Environmental  and  geophysical  changes  can  result  in  the  actual 
changing  of  the  value  of  gravity  at  a  site.  These  changes  can  be 
classified  as  either  long  term  which  tend  to  be  of  a  permanent  nature 
and  short  term  which  tend  to  be  of  a  temporary  nature. 
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3.6.1  Lena  Tirn  Effects 

Long  term  variations  in  tha  value  of  gravity  at  a  station  result 
from  geophysical  changes  within  tha  earth.  These  changes  are  caused  by 
a  redistribution  of  the  earth's  masses  resulting  in  a  change  in  the 
value  of  gravity  at  a  station.  Some  examples  of  effects  that  can  cause 
these  long  term  variations  are:  displacement  of  the  earth's  core,  mass 
redi stri buti on  in  tha  crust  and/or  mantle,  changes  in  the  position  of 
the  station,  changes  in  tha  earth's  rotation  and/or  figure  and  changes 
in  the  gravitational  constant  [Boedecker,  1981].  These  long  term 
variations,  generally,  are  caused  by  geophysical  events  having 
unpredictable  effect  on  the  value  of  gravity  at  a  station.  For  this 
reason,  the  modelling  of  these  effects  in  a  gravity  station  network  is 
presently  not  feasible.  The  magnitude  for  these  types  of  effects  may 
be  on  the  order  of  tens  of  ygal/year  [Boulanger,  1979]. 

3.‘6.2  Short  Term  Effects 

Short  term  variations  in  the  value  of  gravity  at  a  station  rasults 
from  such  things  as  earth  tide,  variation  in  the  level  of  the 
groundwater,  and  changes  in  the  distribution  of  the  atmospheric  masses 
[Boedecker,  1981]  are  more  predictable  than  the  long  term  variations 
previously  mentioned,  provided  sufficient  data  is  available.  The  earth 
tide  which  Is  caused  by  the  gravitational  pull  of  the  moon  and  sun  on 
the  earth  can  be  theoretically  be  modelled  to  an  accuracy  of  better 
than  0.81  Ugal  [Heikkinen,  1978]  provided  adequate  information  is 
available  about  the  location  of  the  station  and  the  epoch  of  the 


observation.  Tho  change  in  tho  laval  of  groundwater  below  a  station 
results  in  a  predictable  change  in  the  value  of  gravity.  This  is  the 
same  as  changing  the  mass  distribution  below  the  station.  To  determine 
the  effect  that  changes  in  the  groundwater  levels  has  on  the  value  of 
gravity  at  a  station  requires  detailed  information  about  the  extent  and 
level  of  the  groundwater  at  the  epoch  of  the  gravity  observation.  This 
information  is  generally  not  available.  This  makes  modelling  the 
effect  of  the  changes  in  the  groundwater  level  not  feasible.  The 
magnitude  of  the  change  in  gravity  at  a  station  caused  by  a  change  in 
the  ground  water  level  is  generally  in  the  range  of  10—20  ygal  but  has 
been  reported  to  be  more  than  100  ygal  [Boulanger,  1979]. 

The  short  term  variations  in  the  value  of  gravity  caused  by  changes 
in  the  atmospheric  masses  above  a  station  is  seasonal  in  naturw  with 
estimates  of  this  variation  being  as  large  as  20—30  ygal  [Boulanger, 
1979].  Even  though  this  variation  might  be  able  to  be  modelled,  its 
effect  on  the  gravity  difference  between  two  consecutively  observed 
stations  could  probably  not  be  detected  since  the  time  between  gravity 
observations  is,  generally,  less  than  six  hours  while  the  change  in  the 
mass  distribution  of  the  b.mosphere  is  assumed  to  be  more  gradual, 
taking  on  the  order  of  days  to  weeks  before  changes  in  the  value  of 
gravity  can  be  observed.  This  effect  should  be  considered  when 
absolute  gravity  measurements  are  being  made. 

Besides  short  term  variations  in  the  value  of  gravity  at  a  station, 
there  are  short  term  effects  that  influence  the  observations  made  with 
a  gravity  meter.  These  are  caused  by  variations  in  the  voltage. 
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temperature,  atmospheric  air  prassura,  and  magnatic  field.  The  effect 
thasa  variations  hava  on  the  gravity  meter  observations  depends  on  the 
individual  instrument  being  used.  Under  controlled  condi tions.  the 
effect  of  each  appaar  to  ba  predictable  (Kiviniemi,  1974].  However, 
the  magnitude  and  direction  of  tha  affect  requires  additional 
information,  such  as  tha  voltage  of  tha  power  supply,  atmospheric 
pressure,  temperature,  and  alignment  of  the  gravity  meter  relative  to 
magnetic  north,  which  is  not  generally  recorded  when  the  gravity  meter 
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4.X  Introduction 

A  Gravity  Raso  Station  Network  conaiata  of  a  aot  of  rocovorablo 
gravity  atationa  distributed  ovar  a  goographical  rogion  for  which  tho 
value  of  gravity  haa  baan  determined.  Thara  are  two  methods  available 
for  determining  the  value  of  gravity  of  a  atation  in  a  network.  One 
method  of  determining  a  station's  gravity  value  is  by  direct 
measurement  of  gravity.  This  is  done  by  Instruments  referred  to  as 
absolute  gravity  measuring  apparatuses  or  absolute  gravity  meters. 

There  are  two  types  of  absolute  gravity  maters,  permanent  and 
transportable  (portable)  which  either  employ  the  free  fall  or  the 
symmetrical  free  rise  and  fall  technique  [Sakuma,  1974].  With  the 
permanent  absolute  gravity  maters,  claims  for  their  accuracies  or 
precision  on  the  order  of  a  few  ygal  are  made;  while  with  the  portable 
absolute  gravity  maters,  accuracies  or  precision  in  the  neighborhood  of 
10  ygal  are  obtainable  {Marson  and  Alasia,  1974;  Marson  and  Alasia,  , 
1980;  Milcox,  1980]. 

The  other  mathod  available  is  by  making  relative  gravity  meter  ties 
from  stations  of  known  gravity  values  to  other  stations.  The  most 
common  gravity  mater  used  in  tho  geodetic  community  for  this  purpose 
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are  tha  LaCoste  5  Romberg  * 0 '  and  '0*  gravity  meters.  See  Section  1.2 
-for  a  description  of  the  basic  difference  between  the  LaCoste  A  Romberg 
*0*  and  *0*  gravity  meters. 

The  portable  gravity  meters,  such  as  LaCoste  A  Romberg  •  G*  gravity 
meters,  provide  only  relative  gravity  difference  information  about  a 
network  and  nothing  concerning  the  actual  gravity  value  of  a  station. 
The  absolute  gravity  meters  on  the  other  hand  provide  information  about 
the  actual  value  of  a  station.  Once  a  station's  gravity  value  is 
known,  portable  gravity  meters  can  be  used  to  make  ties  between 
stations  with  known  gravity  values  and  stations  with  unknown  gravity 
values.  However,  the  gravity  difference  can  only  be  deduced  if  the 
relationship  between  the  counter  unit  difference  and  their  actual 
gravity  difference  is  known.  The  relationship  is  purported  to  be  given 
by  the  Calibration  Table  1  supplied  with  each  gravity  meter.  As  it  has 
been  mentioned,  the  scale  for  this  calibration  table  comes  from  an 
assumed  gravity  difference  between  two  stations  in  New  Mexico, 
Cloudcroft  and  La  Lux.  Zf  the  assumed  gravity  difference  between  these 
two  stations  is  in  error,  then  the  scale  factor  determined  from  their 
assumed  difference  would  cause  the  Calibration  Table  1  to  be  off  by  a 
scale  factor. 

To  determine  the  scale  factor  that  is  to  be  applied  to  Calibration 
Table  1  requires  the  knowledge  of  at  least  one  gravity  difference. 

This  can  only  come  from  the  difference  between  two  stations  of  known 
gravity  value.  A  station  of  known  gravity  value  is  often  referred  to 
as  an  absolute  station.  Having  more  than  two  absolute  stations  does 


net  necessarily  provide  more  information  concerning  the  linear  scale 
factor  that  is  being  attempted  to  be  determined.  The  distribution  of 
the  known  gravity  stations  over  the  range  of  the  network  plays  a  very 
Important  role  in  how  well  the  scale  factor  can  be  determined  [Uotila, 
1976].  Xn  order  to  obtain  the  best  determi nati on  of  the  linear  scale 
factor  requires  that  Interpolation  instead  of  extrapolati on  be  done. 

Xn  addition,  the  shorter  the  Interpolation  interval,  the  better.  Xt  is 
clear  that  for  a  network  having  a  gravity  difference  between  the 
gravity  station  with  the  largest  gravity  value  and  the  gravity  station 
with  the  smallest  gravity  value  of  X  mgal  and  consisting  of  n  known 
gravity  stations,  the  scale  factors  are  best  determined  when  the 
gravity  difference  between  known  gravity  stations  is  approximately 
X/(n— 1)  mgal  with  n>2. 

4.2  fiaoir&i  ai  Ptfr-wrA 

If  the  control  for  a  network  is  not  good,  then  the  results  of  the 
network  adjustment  can  not  be  expected  to  be  good.  The  absolute 
gravity  stations  in  a  network  provide  the  control  for  the  network.  Xn 
the  U.  S.  Gravity  Gass  Station  Network,  two  different  absolute  gravity 
measuring  devices  were  used.  One  was  from  Italy  [Marson  and  Alasia, 
1976;  Harson  and  Alasia,  1960]  and  the  other  was  from  the  United  States 
(Hammond  and  Xliff,  1976].  The  accuracy  of  the  determination  made  with 
each  of  the  Instruments  was  purported  to  be  in  the  neighborhood  of 
10  ugal.  But  the  difference  between  the  values  of  gravity  determined 
at  the  same  station  by  these  Instruments  has  on  occasion  been  as  large 


as  100  Ugal  or  mo rm  [Marson  and  Alasia,  1970  and  1900;  Hammond,  1901, 
letter  to  Uot  1  la].  Tha  reason  -for  such  large  differences  between  the 
determinations  made  by  these  instruments  is  unknown.  Even  more 
puxxling  is  why  repeated  determi nati ons  made  with  the  same  instrument 
do  not  agree  very  well,  fty  comparing  the  values  in  Table  4  and 
Table  5,  examples  of  differences  between  gravity  values  at  the  same 
site  can  be  seen,  situations. 

A  variety  of  reasons  can  be  hypothesised  why  the  two  absolute 
gravity  measuring  devices  give  different  values  of  gravity  at  the  same 
station.  Some  reasons  ara:  there  is  a  systematic  difference  between 
the  two  instruments;  there  is  a  scale  problem  with  the  timing  and/or 
distance  required  for  the  determination;  the  gravity  value  at  the 
station  actually  changed;  and  external  forces  influenced  the 
determination. 

Of  these  possible  reasons,  only  the  last  one  has  been  proven  to  be  a 
real  cause.  Mhen  the  Italian  apparatus  made  measurements  at  the 
Holloman,  AFB  in  New  Mexico  in  May  and  June  of  19S0,  difference  between 
the  two  determinations  of  the  value  of  gravity  at  the  site  of  SO  pgal 
was  noticed.  Further  investigation  as  to  the  reason  for  this 
difference  revealed  that  the  gravity  value  obtained  depended  on  whether 
a  gyro  testing  system  in  a  near  by  building  was  operating.  Mhen  the 
gyro  testing  system  was  not  in  operation,  the  gravity  values  determined 
by  the  Italian  and  the  United  States  absolute  gravity  apparatus  agreed 
very  well.  However,  when  the  gyro  testing  system  was  operating,  the 
values  determined  disagreed  by  approximately  80  ugal.  It  seems  that 
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Table  4  -  Listing  of  values  of  absolute  sites  determined  by  Harson  and 
Alasla. 

Standard 


IG4/GSS 

Station 

Gravl ty 

Error 

Date 

Code 

Name 

In  /| gal 

In  /U gal 

Determl nod 

1S221A 

Boston, HA 

940374659* 

11 

10/  4-11/77 

U994H 

Denver , CO 

979594267* 

10 

10/16-19/77 

15560A 

Bismarck , ND 

940612442* 

10 

10/25-27/77 

11926A 

Hoi 1 Oman. NM 

979139513* 

10 

11/  3-7  /77 

12172A 

San  Francisco, CA 

979972060* 

10 

11/15-17/77 

0415  OC 

Ml  ami , FL 

979004319* 

10 

11/21-26/77 

155 OSD 

Boulder, CO 

979604494 

11 

5/26-27/40 

11926A 

Holloman.NM 

979139544 

12 

6/  2-3  /40 

119S04 

McDonald  Obs.,TX 

974420097 

11 

6/  6-7  /40 

155V04 

Sheri  dan, MY 

940209007 

11 

6/12-14/40 

156E05 

Great  Falls, MO 

940497412 

10 

6/17-14/40 

231A01 

Anchorage, AK 

941924994 

10 

6/27-24/40 

1CB  —  International  Gravity  ftureau 
GSS  —  Geodetic  Survey  Squadron 

*  denotes  Honkasalo  correction  removed  from  published  value. 
Information  obtained  from  (Marson  and  Alasla.  1974  and  1940]. 


Tabla  5  -  Listing  of  values  of  absolute  sites  determined  by  Hammond 


Standard 


XGB/GSS 

Station 

Oravity 

Error 

Date 

Code 

Name 

in  yifgal 

1  n  yffgal 

Determi ned 

11994H 

Denver, CO 

979598277 

10 

3/27-29/79 

119S03 

McDonald  Obs.,TX 

978828655 

8 

7/  3-4  /79 

U92SA 

Holloman, NM 

979139600 

10 

7/  6-7  /79 

119C01 

Tri ni dad, CO 

979330370 

10 

7/10-11/79 

119C03 

Ht.  Evans, CO 

979256059 

8 

7/12-13/79 

155V01 

Casper, MY 

979947244 

25 

7/15-17/79 

135V03 

Sheri  dan, MY 

980208912 

10 

7/18-19/79 

156E05 

Great  Falla, MT 

980497311 

10 

7/ 21-22/79 

U926A 

Hoi 1 oman, NM 

979139600 

8 

5/14,31/80 

15221 A 

Boston, MA 

980378681 

10 

7/  7  /80 

156E05 

Great  Falls, MT 

980497367 

10 

10/  9-11/80 

155V03 

Sheri  dan, MY 

980208964 

10 

10/13-16/80 

15505D 

Boulder , CO 

979608601 

10 

10/18-23/80 

1 19C01 

Trinidad,  CO 

979330393 

10 

10/25-26/80 

119S04 

McDonald  Obs.,TX 

978820087 

10 

10/28-29/80 

15221 A 

Boston, MA 

980378768 

10 

2/  -  /SI 

XGB  —  International  Gravity  Bureau 
GSS  —  Geodetic  Survey  Squadron 

Information  obtained  from  [Hammond,  1*81,  letter  to  Uotila]. 


so 


the  electronics  of  the  Italian  apparatus  wars  affected  by  tha  gyre 
tasting  system  while  tha  Unitad  Statas  apparatus*  alaetrenics  wara  net 
a* feet ad  [Wilcox,  1980]. 


Any  arrers  that  exist  in  tha  valua  of  gravity  at  tha  abseluta 
stations  Mill  ba  raflaetbd  diraetly  into  any  seala  factors  for  tha 
gravity  meters  determined  whan  tha  netMork  adjustment  is  performed.  As 
a  result,  tha  valua  of  gravity  determined  for  tha  other  stations  in  the 
adjustment  Mill  ba  affected. 

Tha  quality  ef  a  gravity  netMork  not  only  depends  on  tha  accuracy  of 
tha  absolute  stations'  gravity  values  but  also  on  tha  distribution  of 
the  absolute  stations  in  tha  netMork  and  tha  gravity  mater  ties  made 
betMaen  gravity  stations.  Tha  question  that  comas  up  is  Mhat  is  th« 
bast  network  conf iguratian  far  a  sot  of  gravity  stations. 


4.2.1  Criteria  for  tha  Bast  Network 


In  order  to  say  that  ana  network  is  batter  than  another,  a  criteria 
must  ba  established  which  will  enable  this  decision  to  ba  made. 

Assuming  there  are  two  networks.  Kj  and  each  containing  tha  same 
stations  but  with  different  gravity  mater  ties  made  between  tha 
stations  and  possibly  different  abseluta  gravity  station,  a  decision  as 
to  which  one  is  preferred,  based  an  the  variance-covariance  matrices 
for  their  adjusted  station  values,  can  be  made  using  one  of  the 
properties  described  by  Fedorov  [1972]. 

Uetila  [1978]  points  out  that  the  most  appropriate  criteria  for 


comparii 


gravity 


invol' 


SI 

ver 1 anca- covari anca  matrix  for  tho  adjusted  station  values  which  has 

the  minimum  trace.  The  trace  of  a  matrix  is  the  sum  of  its  diagonal 

elements.  Therefore,  one  can  say  that  network,  E  ^ ,  is  preferred  to 

network,  E.,  if  the  trace  is  less  than  the  trace  . 

L  t2 

Using  this  criteria,  Uotila  [1978]  describes  a  method  for 
determining  at  which  station  in  a  network  an  absolute  gravity 
measurement  should  be  made  in  order  to  Improve  the  network  the  most. 
This  method  can  be  used  with  a  slight  modification  for  determining 
which  gravity  meter  tie  would  improve  the  network  the  -tost. 

4.2.2  Selection  of  the  Bravi tv  Meter  Tie  to  Improve  a  Network 

Since  a  gravity  meter  tie  provides  Information  about  the  gravity 
difference  between  existing  stations  in  the  network,  the  best  gravity 
meter  tie  to  make  would  be  the  one  which  results  in  the  biggest 
Improvement  in  the  variances  for  the  station  gravity  values.  The 
selection  can  be  made  by  using  a  slight  modification  of  the  method 
Uotila  [1978]  described  for  selection  of  absolute  gravity  sites. 

Assuming  that  the  minimum  variance  solution  for  a  set  of  equations 
for  a  gravity  base  station  network  Is  given  by 

X  -  N^U  «*•!> 

and 

X  -  X.  +  X  <«•*> 


where 


’1 
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Xq  —  Initial  estimated  value  of  parameters, 

X  —  correction  to  X  , 

o 


Xa  -  adjusted  parameter,  values. 

U  —  constant  vector  of  the  normal  aquations, 

..-1 


—  variance— covariance  matrix  of  the  adjusted  parameters, 


ien,  if  a  set  of  observations,  L*,  with  their  vari ance-covari ance 

b 

itrix,  7,2  -  P_ ^ ,  and  a  mathematical  model 
Lh  2 


Lb  ■  F(V 


(4.3) 


e  added  to  the  original  solution,  the  combined  solution  for  the 
rameters  is  given  by 


x;  -  xo +  x2 


(4.4) 


ere 


X2  * 


*(N1  +  A2P2A2)1(U  +  A2P2L2) 


(4.5) 
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As  shewn  by  Uotila  [1978b],  the  vert ance-co variance  aatrix  ef  the 


adjusted  paraaeters  for  the  coabined  solution  Is  given  by 

Ex2  -  Nj1  -  M^A^A^A*  ♦  (4.9) 

3 

The  change  In  the  trace  of  the  variance-covariance  aatrix  for  the 
adjusted  paraaeters  resulting  frea  the  including  of  the  new  gravity 
difference  inforaatien  will  always  be  negative  and  its  aagnitude  is 
gi  van  by 

N^A^(A2N^aJ  *  P21)  lA2Nt!  <4*10> 


What  is  desired  is  the  gravity  difference  which  will  aake  (4.10)  the 


largest.  Assuaing  that  a  single  uncerrelated  gravity  difference,  g... 


can  be  observed  with  a  variance  ef  a *  between  two  station  whose 


’U 


gravity  values  are  given  by  g^  and  g^  where 


91j  "  9,  *  9j 


(4.11) 


then  A^  aatrix  will  be  a  row  aatrix  ef  sere  eleaents  except  for  a  «1  in 
the  1— th  eeluan  and  a  —1  in  the  j— th  celuan.  The  value  ef 


Vi’**  *  fi 


(4.12) 


can  be  shewn  to  be  the  value  ef  the  variance  for  the  i— th  and  j-th 
station  ainus  twice  their  covariance  plus  the  variance  ef  the  new 
observed  gravity  difference.  The  Inverse  ef  (4.12)  is  the  reciprocal 
ef  that  sue.  The  aatrix  product  ef  AgN~^  will  result  in  a  aatrix 
feraed  frea  two  eeluans  ef  one  being  the  i— th  celuan  and  the  ether 


1 


baing  t ho  ntfttlvi  of  the  J— th  column.  If  th«  two  columns  of  tho 


matrix,  »  arm  thought  of  aa  tha  vectors,  Yj  and  Yj , 

given  by  (4.10)  can  bo  shown  to  bo  tho  value  of 


VV 


°i  *  °j  *  laij  *  <4 


(4. IS) 


whore 


•  —  represents  the  dot  product  of  two  vectors, 

—  covariance  between  station  g^  and  g^» 

(j2  —  variance  of  the  new  observed  gravity  difference. 

Therefore,  the  gravity  difference  between  stations  gg  and  gj  which 
maximises  tha  value  of  (4.1S>  would  Indicate  where  a  gravity  tie  should 
be  made  to  Improve  the  network  the  most.  This  method  Is  similar  to  the 
one  described  by  llotlla  [1970b]  for  selection  of  absolute  gravity 
sites. 


In  a  network  of  n  stations,  there  Is  a  possibility  of  n(n~l>/2 
different  gravity  differences  with  at  least  n— l  of  these  differences 
already  existing  In  the  network,  ly  using  this  method,  the  effect  on 
the  trace  of  the  stat1en*s  variance-covariance  matrix  of  adding  a  new 
and/or  existing  gravity  meter  ties  with  a  certain  accuracy  can  be  seen. 

For  every  station  added  to  a  network  of  n  stations,  n  more  possible 
gravity  differences  are  Introduced.  Many  times  these  stations  which 
are  added  to  a  network  are  very  close  to  an  existing  station.  The 
added  stations  will  generally  have  gravity  values  very  close  to  that  of 
the  neighboring  station  that  had  been  previously  established.  As  a 


result,  the  saell  gravity  differences,  generally  less  than  10  egal, 
provide  little  eddltlonal  Information  that  can  be  used  In  Improving  the 
determination  of  the  gravity  meter  scale  factor.  These  types  of 
stations  are  referred  to  as  eccentric  stations.  If  thare  Isa  large 
number  of  eccentric  atattone  In  an  area,  any  Improvement  In  the  area 
Mill  result  In  a  large  laproveawnt  In  the  trace  of  the  stations* 
variances,  but  In  reality  the  actual  Improvement  of  the  variances  for 
stations  outside  the  local  area  could  be  small.  In  such  a  situation,  a 
local  Improvement  results  rather  than  the  desired  over  all  network 
Improvement . 

Therefore,  when  attempting  to  select  where  new  gravity  meter  ties 
should  be  mode.  If  possible,  eccentric  stations  should'not  be 
considered.  Zf  possible,  only  one  station  in  a  local  area  should  be 
used.  The  station  with  the  most  gravity  ties  to  stations  outside  the 
local  area  should  be  used,  provided  that  all  the  eccentric  stations  In 
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CHAPTER  FIVE 
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MATHEMATICAL  MODEL t 


S.l  Tvdm  o f  Obaarvahloa 


A  gravity  meter  raiding  Involves  tha  absarvar  racardlng  tha  valua  of 
tha  count ar  and  tha  position  af  tha  dial  aftar  tha  gravity  meter  has 
baan  nullad,  remembering  that  l  count ar  unit  Is  aqulvalant  to  100  dial 
units.  Tha  obaarvatl an  Is  than  ganarally  cenvartad  ta  Its  valua  In 
mllllnal  by  Interpolating  within  tha  factory  supplied  Calibration 
Table  1.  Thb  manufacturer  suggests  that  a  linear  Interpolation  bo  uada 
within  tha  Calibration  Table  1  to  obtain  tha  valua  In  mllllgal  for  an 
observed  reading  [LaCaSta  I  Romberg,  1*00].  This  can  bo  accomplished 
by  following  tha  siaqilo  procedure  outlined.  First,  tha  valua  In 
mllllgal  far  tha  mauntar  roadl no  which  Is  nearest  ta,  but  lass  than, 
tha  observed  reading  Is  deter ml  nod.  Let  tha  counter  reading  valua  used 
be  Y  and  Its  correspond ng  value  in  mllllgal  bd  X.  Tha  difference,  Z, 
between  tha  observed  reading. and  tha  counter  reading  used,  Y,  is  than 
obtained.  Meat,  multiply  the  difference,  Z,  by  the  factor  for  interval 
for  the  counter  reading,  Y,  and  add  the  result  to  the  value  of 
.mllllgal,  X.  The  result  Is  the  value  of  mllllgal  far  the  observed 
readl ng. 


as 


Using  a  llnaar  Interpolation  proeoduro  within  tho  Calibration 
Table  1  assume  that  tho  relationship  between  the  counter  readings  and 
the  values  in  siilllgal  can  be  adequately  represented  by  a  pioce-wlsa 
linear  function,  tinea  the  Calibration  Table  1  Information  is  produced 

t 

by  reading  the  value  of  factor  for  interval  off  of  a  continuous  curve 
of  factory  calibration  scale  factor  data.  It  Is  possible  that  a 
piece-wise  continuous  relationship,  such  as  a  cubic  splino  [Spath, 
1964],  should  be  used.  When  a  comparison  was  done  to  see  the 
difference  between  the  use  of  a  linear  or  cubic  spline  Interpolation 
procedure  within  the  Calibration  Table  1,  the  maximum  observed 
difference  between  the  two  methods  was  on  the  order  of  1  to  S  figal. 

This  shows  that  the  behavior  of  the  Calibration  Tablo  1  Is  relatively 
smooth. 

Since  the  Calibration  Table  1  attempts  to  represent  a  continuous 
type  of  continuous  Information,  the  cubic  spline  method  of 
Interpolation  would  be  preferred  to  the  linear  method  of  Interpolation 
even  though  the  difference  between  the  two  methods  Is  small. 

As  mentioned  previously,  the  values  In  mllligal  are  derived  from  the 
factor  for  Interval  values  based  en  the  assumption  that  the  factor  for 
Interval  ever  a  particular  interval,  usually  100  counter  units.  Is 
constant.  2n  reality,  this  Is  net  exactly  true.  The  factor  for 
Interval  values  represent  a  continuous  function  as  opposed  to  a 
piecewise  linear  function.  Therefore,  the  values  In  mllligal  should 
be  determined  from  the  Integration  of  this  continuous  function  or  an 
approximation  of  It,  such  as  a  cubic  spline.  The  difference  between 


tli*  TiMf  I  values  In  el 111  gal  and  thn*  computed  fro*  i 

pife<«  hMm  rapPMWtitlM  pf  the  factor  far  Intarvat  can  be  as  large 
as  Itwif  ufal  avar  a  •••  ays!  Interval*  The  largo  differences  can  lead 
fa  a  syatsmatle  aarar  being  Introduced  In  large  gravity  differences, 
aver  1M  agai,  which  weald  net  N  evident  In  sealler  differences.  This 


eeans  that  as  the  gravity  difference  Increase*  in  size  the  variance  af 
the  difference  adept  the  eeen  ef  the  gravity  difference  that  is 
acceptable  should  slse  be  increased.  This  makes  the  detection  ef  tare 


ef  e  given  sip#  ears  difficult  far  gravity  difference#  that  are  large 
as  appeaed  te  these  that  ere  smell. 

Once  the  value  in  eilllgel  far  all  readings  has  been  detereined,  the 
gravity  dlffereepe  between  twe  stations  can  be  detereined  by  computing 
the  difference  between  their  values  In  allligal.  However,  the 
resulting  difference  will  net  be  Independent  ef  the  tine  ef  the 
observation  and  eight  have  te  be  sealed  by  a  factor  te  obtain  the 
proper  units  ef  gravity.  The  difference  oan  be  eade  tiee  Independent 
by  rpeeving  any  known  tiee  dependent  effect  such  a*  the  earth  tide 
effect.  In  order  te  ebtele  a  gravity  difference  in  the  proper  units, 
the  absolute  seals  faster  applied  te  create  the  Calibration  Table  1 
scat  be  correct  and  valid  for  thp  range  in  which  the  gravity  meter  was 
being  used.  It  oust  be  rsmswbered  that  the  absolute  scale  factor 
applied  te  the  relative  scale  factor  during  the  factory  calibration 
procedure  was  truly  valid  only  over  a  range  of  approximately  242  mgals 
for  gravity  values  In  the  region  ef  t7f  gal.  Assuming  it  to  bo  valid 
for  any  range  over  which  the  gravity  meter  is  being  used,  might  not  bo 


•9 


eorraet . 

It  la  oUar  that  tha  only  true  ebsarvad  quantity  for  a  gravity  Motor 
la  tho  obaorvod  oountor  and  dial  reading.  From  tho  obaarvod  eountor 
readings,  their  values  In  allllgal  are  derived.  And  finally,  the 
gravity  difference  between  two  values  In  nllllgal  can  ba  obtained. 

S.1.1  Observed  Counter  Raadlno 

The  quantity  that  Is  observed  on  a  gravity  oeter  Is  the  position  of 
the  counter  and  the  dial  which  when  eeablned  yields  the  observed 
counter  reading.  The  problea  with  using  this  quantity  as  the 
observable  In  a  Mathematical  model  is  finding  an  analytical 
relationship  that  will  tranform  the  counter  readings  Into  units  of 
gravity.  An  empirical  relationship  exists  In  the  form  of  the 
Calibration  Table  1  which  Is  supplied  with  each  LaCoste  A  Romberg  »G» 
gravity  meter.  It  Is  necessary  to  determine  If  there  Is  an  analytical 
expression  for  the  amplrlelal  relationship  expressed  In  the  Calibration 
Table  1  which  can  be  used  as  a  functional  relationship  In  a  least 
squares  adjustment  medal.  If  net.  It  Is  necessary  to  determine  If 
there  Is  some  analytical  function  which  approximates  that  empirical 
relationship.  If  the  relationship  can  be  expressed  adequately,  then 
the  observed  counter  readings  can  be  used  directly  as  observables  In  a 
least  squares  adjustment  model. 

Assuming  that  such  an  analytical  relationship  can  be  found  and  that 
It  Is  a  simple  function,  then  solving  for  periodic  screw  effects 
becomes  possible  without  getting  Into  the  problem  of  having  to  assign 
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Heights  to  two  dHfsrsnt  observables*  counter  readings,  and  their 
values  in  alUlgal  Mhloh  represents  the  a  see  quantity. 

1.1. S  Value  In  Ml  11 Inal 

It  the  observed  counter  reading  cannot  be  used  as  observables  In  the 
adjustnent  due  to  the  lack  of  an  adequate  Model  ter  the  Calibration 
Table  1,  then  the  next  beat  quantity  that  could  be  used  would  be  the 
values  In  nllllgal  tor  the  observed  counter  readings.  A  Model  which 
uses  ae  Its  observable*  the  value  In  etlllgal,  has  been  developed  by 
Uetlla  (1974]  and  used  with  sows  success.  The  Model  which  Uotlla 
(1974]  proposed  used  an  equation  Involving  the  difference  between  two 
gravity  Meter  observations  which  laplles  the  following  expression  for 
oaeh  gravity  water  observation 

n  , 

I  0.x  +  k(t  •  T  )  ♦  S  •  6  -  a  *».i> 

fio  1  ° 

where 

ft  I  —  coefficient  of  the  1— th  order  scale  factor  tera* 

»  —  value  In  Ml  111  gal  of  the  observed  counter  reeding  corrected 

for  all  known  systOMatle  effects*  such  ai*  Mirth  tides  an 
height  of  instruMent  above  the  station* 
h  —  eeefflctent  of  the  drift  tero* 
t  -  epoch  of  the  observation* 

T  -  sane  arbitrary  Initial  epoch  associated  with  the  set  of 
o 

gravity  rat  or  oboorvatlora. 


S  —  an  unknown  of foot  valuo  that  auit  bo  added  to  obtain  tho 
oorract  absolute  gravity  valuo  for  tho  station, 

0  —  gravity  valuo  of  tho  station, 

n  —  largest  ordor  scale  factor  taro  to  bo  Included  In  tho  aodel. 

•oeauso  tho  linear  scale  factor  taro  In  addition  to  scaling  tho 

valuo  In  mllllgal  would  also  scale  all  of  tho  systematic  effects,  this 

model  Is  theoretically  Ineorroet  ovon  when  n  =  1 .  However,  this  effect 

would  be  minimal  because  the  value  of  the  linear  scale  factor  tore  Is 

very  close  to  1.  Per  n  >l,  the  question  of  what  the  Model  represents 

become  very  confusing.  As  noted  by  (loti la  [1974],  higher  order  scale 

factor  terMS  eight  not  have  any  physical  Justification,  but,  even  If 

they  were  physically  Justified,  a  similar  problem  exists  when  n > 1  as 

when  nsl.  If  x  Is  expressed  as  *  ♦  C  where  i  Is  the  value  In 

ollllgal  of  the  observed  counter  reading  obtained  free  the  Calibration 

Table  1  and  C  Is  all  associated  systematic  effects,  then  when  n  s  2, 

2  2 

for  example,  terms  Involving  D  C  and  C  enter  Into  the  model. 

Theoretically,  these  terms  are  net  appropriate;  however,  their  affect 

„  2 

would  be  small  compared  to  the  value  of  c 

A  more  appropriate  expression  for  equation  S.l  might  be 

n  . 

^DjZ1  ♦  k(t  -  To)  +  S  +  C  -  6  -  0  <*•*> 

where  the  meaning  of  the  symbols  has  baen  previously  defined.  Hlth  the 
expression,  the  problem  of  scaling  the  systematic  effects  Is  removed. 


5.1.3  to MXlS*  aifliCMfiM 


Tho  aest  eoaaonly  used*  and  probably  nest  Incorrectly  usod 
obsorvablo  Is  tho  gravity  difference.  Many  different  Medals  have  boon 
proposed  using  tho  gravity  difference  as  tho  obsorvablo  [Mhalon,  1974; 
HcConnoll  and  Oantar*  1974;  Torgo  and  Kannglosor.  19301.  A  typical 
Medal  Might  bo  expressed  as 

lAfljj  +  k(tj  -  tj)  +  C,  -  Cj  -  Gj  +  6j  -  0  <5.3> 

where 

1  —  coefficient  of  the  linear  scale  factor  tore, 

-  ebsorved  gravity  aeter  difference  In  Mllllgal  between  station 
1  and  ji 

h  —  coefficient  of  the  drift  tore* 

t  —  tlae  of  observation  at  station  n, 
fl 

0n  —  gravity  value  at  station  n, 

Cn  —  all  systeMatle  effects  for  the  observation  associated  with 
the  observation  at  n. 

equation  (5.3)  can  be  rewritten  Into  true  observation  fora  as 


the  observation  fora  exists  when  the  observed  quantity  Is  strictly  a 
function  of  parameters 
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The  lttit  squares  solution  of  this  system  of  aquations  roquiros  tho 
var 1 anco— covarl anco  Matrix  of  tho  gravity  differences  computed  from 
observations  to  be  known.  Generally,  this  variance-covariance  matrix 
is  assumed  to  be  diagonal  which  neglects  the  correlation  that  exist 
between  successive  gravity  differences  as  stated  by  Uotlla  [1974]. 

Zf  the  correlation  In  the  variance-covariance  matrix  of  tho  gravity 
differences  Is  not  included,  the  effect  on  tho  value  of  the  adjusted 
parameter  values  will  probably  bo  vary  small.  However,  the 


variance-covariance  matrix  for  the  adjusted  parameters  will  be  greatly 
affected.  As  a  result,  the  estimate  of  the  variances  for  the  adjusted 
parameters  will  be  too  low.  This  results  In  an  over  optimistic 


estimates  of  the  standard  error  of  all  parameters. 

5.2  Characteristics  a±  Calibration  HJfeli  1 

Kach  "0"  gravity  meter  produced  Is  a  hand  built  product  which 

/ 

results  In  Its  own  characteristic  behavior.  Even  though  this  behavior 
Is  different  for  every  gravity  meter,  the  Calibration  Table  1  which 
represents  this  behavier  does  exhibit  some  common  characteristics.  The 
functional  relationship  between  the  value  in  mllligal  and  the  counter 
readings,  although  unknown.  Is  nearly  linear.  Zf  the  linear  trend  is 
removed  from  the  Calibration  Table  1,  the  higher  order  trends  can  be 
seen  mere  readily.  Attempts  have  been  made  to  represent  the  entire 
Calibration  Table  1  by  using  polynomial  functions  up  to  the  fifth  order 
with  some  success  [Uotlla,  1975b].  Problems  arise  in  the  separation  of 
the  eoefficents  of  these  higher  order  polynomials,  since  the 
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Mfflemits  tend  to  bo  highly  correlated  with  each  other.  Correlations 
k  largo  as  0.999  and  largor  aro  not  uncommon.  With  those  high 
arrelations,  tho  least  squares  polynomial  fit  tends  to  be  very 
■nsltlve  to  the  data  used  In  Making  the  fit.  Since  the  actual 
unetlonal  behavior  Is  unknown,  using  these  higher  order  polynomials  to 
•present  the  relationship  eould  very  easily  lead  to  fitting  the  data 
nstead  of  the  Identification  of  any  general  behavior  characteristics. 

.2.1  Possible  Models 

After  the  linear  trend  has  bean  removed  from  the  Calibration 
ible  1,  It  la  apparent  that  the  residuals  of  the  value  In  mllllgal 
xhlblt  a  definite  trend,  a  type  of  wave  or  periodic  function.  To 
ttermlne  the  characteristics  of  the  trend  exhibited,  a  spectral 
talysls  of  the  values  In  mllllgal  versus  the  counter  readings  for  each 
•avity  meter's  Calibration  Table  1  was  performed.  The  analysis 
idleated  that  the  only  sinusoidal  term  present  was  a  low  frequency 
trm.  The  wave  length  ef  this  term  was  In  all  cases  larger  than  4000 
•unter  units.  The  median  value  and  mode  of  the  wave  length  for  all 
te  gravity  meter  was  around  S000  counter  units.  See  Figure  19  for 
temples  of  this  trend.  One  possible  relation  other  than  a  higher 
der  polynomial  that  eould  be  used  te  model  this  type  of  trend  Is  a 
nusoidal  relation  given  by 

Bj  sin(  M  )  +  B2  cos (  )  CS.S) 


ere 
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d  —  th«  counter  reading, 

T  -  the  period  of  the  sinusoid, 

B1'B2  ~  th*  •■plitude  of  the  sinusoid. 

This  relationship  can  bo  expressed  as  a  single  function  ty  using  a 
basic  trigonometric  identity  giving 

A  COs(  +  u>  )  (5.5) 

where 


A  —  amplitude  of  the  sinusoid, 

d  —  observed  counter  reading, 

T  -  period  of  that  sinusoid  in  counter  units, 
a)  —  phase  angle  -for  the  sinusoid. 


To  see  hew  well  the  Calibration  Table  l*s  values  in  mi  111  gal  can  bo 
represented  as  a  function  of  the  counter  readings,  a  least  squares 
adjustment  was  performed  using  a  model  containing  a  polynomial  and 
sinusoid  terms  as  given  by 

l  l,d*  +  l  A.  cos (  ~~  +  w.  )  -  y  C5.7) 

1-0  '  J-1  J  Tj  J 

where 

y  —  value  In  mi 111  gal  from  Calibration  Table  1, 
d  —  counter  reading  corresponding  to  the  value  in  mi  Hi  gal, 

Aj  —  amplitude  of  the  1— th  sinusoid  in  mi  111  gal, 

Tj  —  period  of  the  1— th  sinusoid  in  units  of  d, 

«j  -  phase  angle  of  the  1-th  sinusoid  in  radians. 
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n  —  tht  order  of  tho  polynomial, 

k  —  tho  numbor  of  linuioldt. 

Various  valuos  of  n  and  k  mops  soloctod  for  a  numbor  of  counter 
readings  rangos  to  sea  how  well  tho  Calibration  Tablo  l*a  valuos  In 
mllllgal  could  be  represented  analytically  as  a  function  of  the  counter 
readings.  A  sample  of  the  results  for  various  models  tested  Is  given 
In  Table  A  and  Table  7.  Table  A  represents  the  situation  in  which  a 
gravity  meter  Is  used  over  a  world  wide  range  of  5.4  gal  with 
observations  made  within  a  range  of  the  counter  readings  from  700  to 
S300  counter  units.  The  uses  of  a  linear  model,  n  *  1  and  k  *  0, 
results  in  root  mean  square  (RMS)  fit  of  the  57  residuals  for  each 
gravity  meter  ranging  from  approximately  200  to  1200  ugal.  However, 
when  a  sinusoid  was  Included  In  the  model  along  with  the  linear  term, 
n  s  1  and  k  e  1,  the  RMS  fit  of  the  57  residuals  was  reduced  to  a  range 
of  approximately  IS  te  55  ygal.  But  for  the  data  used  In  this  study 
which  represents  only  a  limited  range  of  counter,  the  counter  readings 
from  about  2000  to  4400  counter  units  were  used.  Zf  only  that  part  of 
the  Calibration  Tablo  1  Is  used,  then  the  situation  gets  better  as  is 
shown  In  Table  7.  The  number  ef  residuals  used  to  produce  the  results 
given  In  Table  7  was  25. 

Table  7  compares  four  different  models,  three  different  polynomial 
models  and  a  sinusoid  or  often  referred  te  as  periodic  model.  As  can 
be  seen  In  Table  7,  for  every  Calibration  Table  1  selected,  except  for 
gravity  metar  *5-157*,  the  sinusoid  model,  n  a  1  and  k  a  1,  resulted  in 
the  beat  RMS  fit  ef  the  data.  Zn  most  cases,  the  difference  between  the 


Tabl*  ft  -  SuMiry  of  the  requite  of  Modelling  tho  Calibration  Table  1 

count or  roadingo  fro*  700  to  4300  using  a  linaar  function  and 
a  linaar  function  plus  a  par iodic  tera. 


Gravity 

Calibration 

Root  Moan 

Square  Pit  in  ugal 

Natar 

Data 

Linaar 

Periodic 

0-10 

28/10/40 

1192.28 

41.02  (7480) 

8-44 

25/04/ftl 

433.31 

32.83  (4381) 

6—48 

23/03/44 

768.22 

31.00(10222) 

8-81 

07/08/44 

1134.44 

14.48  (8430) 

6-103 

07/09/48 

468.47 

14.44  (4832) 

6-103 

14/11/79 

731.74 

18.23  (7018) 

6-111 

28/03/44 

904.01 

12.43  (7924) 

6-1  IS 

14/08/78 

212.04 

80.70  (4048) 

6—126 

24/10/78 

1014.92 

14.23(10288) 

6-131 

18/06/78 

320.29 

84.42  (4380) 

6-142 

14/03/47 

288.27 

42.47  (4283) 

6— 1S7 

10/08/47 

1029.44 

48.1.4  (  8428) 

6-220 

14/11/79 

1011.04 

18.34  (8036) 

6— 2«8 

18/06/78 

902.27 

12.41  (8493) 

6—240 

11/10/78 

344.88 

28.81  (6924) 

Tha  quant It in*  in 
eountar  unita 

paranthama*  danota  tha  period  of  tha  periodic  tara 

Modal  usad  la  givan  by  aquation 

(8.7)  where 

PERIODIC  implies  n*l  and  kal. 


Table  7  —  Summary  of  the  results  of  Modelling  the  Calibration  Table  1 
counter  readings  -from  2000  to  4000  using  various  order 
polynomials  and  a  linear  function  plus  a  periodic  term. 


Gravity 

Calibration 

Root 

Meter 

Date 

Linear 

G— 10 

25/10/40 

352.31 

0-44 

25/04/43 

209.54 

G— 40 

23/03/44 

131.03 

0— 01 

07/03/44 

293.43 

0-103 

07/09/45 

193.91 

£-103 

14/11/79 

214.31 

G-lll 

25/03/44 

243.94 

G— 115 

14/05/73 

34.23 

0-1 23 

24/10/75 

243.53 

0—131 

15/05/73 

125.31 

0-142 

14/03/47 

123.35 

0-157 

10/03/47 

229.13 

0-220 

14/11/79 

244.25 

0-243 

15/05/73 

243.40 

0-249 

11/10/73 

121.34 

The  quantities 

denote 

counter  units. 

Model  used  Is  given  by  equation  (5.7) 
SECOND  Implies  nsl  and  ksO*  THIRD  Imp 
Implies  nsl  and  k*l. 


Mean  Square  Pit  Inx/gal 


Second 

Third 

Perl odl  e 

27.41 

23.93 

7.44 

(4431) 

11.49 

11.59 

3.49 

(5104) 

14.34 

14.23 

13.40(10494) 

14.19 

9.85 

9.83(24155) 

13.43 

9.21 

5.17 

(5920) 

13.43 

3.90 

0.99 

(8447) 

14.24 

4.13 

3.30(11244) 

4.24 

5.39 

1.41 

(4754) 

14.30 

9.04 

8.54(10349) 

14.23 

4.48 

0.75 

(4049) 

21.34 

9.94 

4.03 

(4721) 

14.34 

4.73 

7.32 

(3248) 

9.99 

5.05 

2.45 

(9014) 

30.49 

4.12 

2.11 

(7551) 

14.23 

4.49 

2.42 

(4442) 

the  period  of  the  periodic  tern  In 

where  LINEAR  Implies  nsl  and  ksOt 
lies  nsS  and  ksO  and  PERIODIC 


Srd  order  polynoalal  and  the  sinusoid  aodel  mss  relatively  saall. 
HoMever,  there  mss.  generally,  s  notlcoablo  laprovaaent  In  the  RMS  fit 
as  the  degree  ef  the  polynomial  Increased. 


tt  Is  knoMn  that  the  Calibration  Table  1  for  each  gravity  aeter  mss 
produced  by  an  arbitrary  aethed  ef  curve  fitting  ef  the  relative  scale 
factors  frea  the  factory  calibration  procedure.  Also,  In  general,  the 
least  count  ef  the  counter  reading  recorded  Is  equivalent  to 
approxlaataly  18  ygal.  Frea  this  Inferaatlen.  It  Mould  appear  that 
either  a  sinusoid  or  Srd  order  polynoalal  aedal  Mould  be  a  good 
analytic  representation  ef  the  Calibration  Table  1  Inferaatlen  In  the 
range  ef  the  counter  readings  observed  along  the  United  States 
Mid-Continent  Calibration  Lina.  Tho  question  ceaes  up  as  to  Mhlch 
aedel  should  be  used. 

5.2.2  Sinusoidal  Medal 

The  reason  a  sinusoidal  aedal  Mas  selected  ever  a  Srd  order 
polynoalal  aodel  Is  based  on  the  fact  that,  after  the  linear  trend  is 
reaeved,  the  slnuseld  represents  a  staple  curve  Mhere  as  the  Srd  order 
polynoalal  Is  a  eaaposlto  of  tMO  staple  curves,  one  a  function  ef  x2 
and  the  ether  a  function  ef  k *.  This  ashes  the  polynoalal  aedel 
soaoMhat  undesirable  because  ehanges  in  the  coefficients  ef  the  2nd  and 
Srd  order  teras  are  difficult  to  relate  to  ehanges  in  the  basic 
characteristic  behavior  ef  the  curve  It  Is  representing.  The 
sinusoidal  aedel  Is  easier  to  interpret  Mhat  changes  in  the  value  ef 
the  aaplitude,  period  or  phase  angle  represent,  ly  using  the 


sinusoidal  nodal,  tho  eharactorl stl c  behavior  of  the  Calibration 


Tablo  1  can  bo  classlflod  by  tho  period  of  tho  sinusoid  Mhlch  boat 
roprosonts  tho  Information. 

Tho  determination  of  tho  coefficients  of  tho  sinusoidal  modal  given 
by  equation  <5.7>  whore  n  a  1  and  k  s  1  was  dona  by  a  least  squares 
adjustment  of  tho  sot  of  non— linear  observation  equations.  For  tho 
adjustment,  tho  counter  readings  wore  assumed  to  bo  without  error  and 
tho  variance  associated  with  each  Independent  and  uncorrolatod 
observation  (value  In  mllllgal)  was  assumed  to  be  of  tho  same  constant 
value.  The  solution  was  Iterated  until  the  change  In  V*FV  was  less 
thajt  1  part  In  10  billion.  To  assure  that  the  adjusted  value  of  the 
period  of  the  sinusoid  would  represent  a  long  wave  length,  the  period 
was  weight  constrained  to  permit  a  proper  solution  to  bo  achieved.  Th 
actual  value  of  the  weight  needed  depended  on  tho  initial  estimates  of 
tho  coefficients.  The  determination  of  tho  weights  and  the  good 
estimates  for  the  coefficients  Involved  a  lot  of  trial  and  error.  It 
was  not  uncommon  to  perform  180-S80  Iterations  before  a  solution  was 
achieved  based  on  the  desired  accuracy  of  V*FV.  The  large  number  of 
Iterations  required  Indicated  there  was  a  very  strong  dependency  that 
exists  between  the  period  and  phase  angle  In  the  sinusoidal  model. 

With  the  period  and  phase  angle  being  so  highly  correlated.  It  is  not 
feasible  to  Include  both  quantities  as  parameters  In  an  adjustment 
model  and  expect  It  to  converge  rapidly.  Therefore,  either  the  period 
or  phase  angle  must  bo  fixed  or  heavily  constrained  within  an 
adjustment. 
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Whan  the  par  lad  fa  heavily  constrained  which  assantlally  fix  Its 
value*  a  solution  can  easily  ho  achieved  In  a  few  Iterations,  by 
obtaining  solutions  for  a  nuaber  of  different  periods*  graphs  similar 
to  Figure  20  or  Figure  21  can  be  produced  for  each  gravity  ooter  which 
show  the  approximate  relationship  between  values  of  the  eoeffl cents  of 
the  sinusoidal  nodal  and  the  value  of  tho  RMS  of  the  residuals.  As  can 
bo  soon  In  either  Figure  20  or  Figure  21*  changes  In  tho  value  of 
par led  near  tho  solution  with  tho  nlnlnun  OHS  of  tho  residuals  results 
In  little  change  In  tha  RMS  of  tho  residuals  but  naans  largo  changes  In 
the  phase  angle*  aaplltuda  and  linear  scale  factor  toms.  Using 
graphs*  slollar  to  Figure  21*  eat 1 nates  for  tha  value  of  tho  sinusoidal 


nodal  coefficients  can  bo  determined  for  different  period  of  tho 
sinusoid.  From  those  graphs*  a  period  near  tho  nlnlnun  RMS  of  tho 


residuals  was  selected  for  each  gravity  motor.  Since  tho  value  of  tho 
ported  has  little  effaet  on  tho  RMS  of  tho  residuals*  all  periods  used 
id  tho  final  adjustment  were  rounded  off  to  tho  nearest  thousand  of 
counter  units.  Te  see  If  the  rounding  off  would  hove  any  effect  an  th< 
final  adjustments*  a  number  of  preliminary  adjustments  were  done  for 
various  values  of  tho  periods  near  tho  rounded  off  valuo.  Tho  rosults 
Indicated  That  tho  least  squares  seultlon  for  tho  model  used  would  not 
bo  effected  by  using  the  rounded  off  periods. 


S.S  Mathematical  Modnla  for  Brevity  Maters 


The  basic  equation  usod  for  the  VBV  gravity  meter  is  of  tho  fora 
given  by  equation  C2.1X  where  C  Includes  such  things  as  the  correction 


Flfun  20  —  Tha  result!  sf  a  nusbsr  af  laast  squares  fits  a*  a  linaar 
funettan  plus  a  par ladle  tans  to  tha  data  from  tha 
Call  Prat  tea  Tabla  1  datad  17  Oetabar  1977  far  gravity  aatar 
0—41  far  aalaetad  valuas  af  tha  par lad. 


«  Plguro  21  -  Tho  PMultt  of  a  mmbtr  of  laast  squares  fits  of  a  Knur 

—  function  plus  a  periodic  torn  to  tho  data  fro*  tho 

Calibration  Tablo  1  da tod  14  Novoobor  1979  for  gravity 
•otor  0-220  for  aoloetad  valuoa  of  tho  porlod. 
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for  oarth  tides,  the  correction  for  height  of  Instrument  above  the 
station,  the  correction  for  a  possible  time  dependent  drift  of  the 
Instrument  and  It  could  Include  the  correction  for  other  environmental 
or  geophysical  effects  such  as  changes  In  the  ground  mater  level  and 
changes  In  the  atmosphere  above  the  station. 

Zn  order  to  make  the  gravity  meter  observations  Independent  of  the 
epoch  of  the  observation,  any  time  dependent  effects  must  be  Included 
In  C.  One  time  dependent  effect  that  must  be  removed  Is  caused  by 
gravitational  attraction  of  the  sun  and  the  moon,  commonly  referred  to 
as  the  earth  tide. 

S.S.l  larth  Tide 

The  determination  of  the  earth  tide  Is  accomplished  by  computing  tho 
effects  of  the  sun  and  the  moon  on  the  gravity  value  of  tho  station  at 
the  epoch  of  the  observation. 

The  program  used  to  compute  the  earth  tide  was  obtained  from  the 
Finnish  Seodetlc  Institute.  The  program  computes  the  vertical 
component  of  the  tidal  force  to  an  accuracy  of  S.01  ygal  basod  on  tho 
assuaptlon  of  a  completely  rigid  homogeneous  ellipsoid  of  revolution 
(Melkklnen,  1975].  Since  the  earth  Is  not  really  a  completely  rigid 
homogeneous  ellipsoid  of  revolution.  In  order  to  obtain  a  realistic 
value  for  the  vertical  component  of  gravity,  the  tidal  force  obtained 
Is  sealed.  The  value  ef  this  scale  factor  depends  on  regional 


conditions,  but  a  world-mi  do  conventional  value  of  1.15  was  recommended 


Canberra,  Australia*  lHt  [ZAO*  1980].  Zn  tha  eoaputatlon  a f  tha  aarth 
tide  ‘far  aaeh  ebservatl on,  a  seal  a  factor  af  1.11  mss  applla«l  ta  tha 
eeapitad  vartloal  oaapanant . 

Tha  laplaaontatlon  af  tha  Finnish  aarth  tlda  prograa  as  raealvad  and 
given  In  Molfcklnen  [1970]  proved  te  ha  vary  Inefficient  as  far  as 
Input/eutput  <Z/0>  was  eeneernad.  Ay  aaklng  soaa  staple  Z/0  changes 
which  basically  Involved  aedlfylng  tha  subroutine  OETREC  sa  that  tha 
Z/0  was  dona  unforaatted  and  anly  once  during  tha  execution-  of  tha 
pragraa*  tha  eeoputatlon  tlaa  required  ta  eaaputa  tha  aarth  tide  waa 
reduced  by  a  factor  of  approxlaataly  10.  After  tha  aodlfl cations  ware 
aada*  It  took  approxlaataly  0.S  aeeonds  par  aarth  tide  correction  on  an 
Aadhal  A70  V/A-ZZ  eoaputar  operated  by  tha  Znstruetlonal  and  Rasaarch 
Coaputar  Canter  (ZRCC)  at  Tha  Ohio  State  University. 

One  other  alner  prablaa  Involving  tha  eoaputar  prograa  waa  noticed 
and  corrected  In  tha  version  af  the  pragraa  used  In  this  study  which 
should  ba  aentlened.  Tha  pragraa  requires  inforeatlon  about  the 
difference  between  ophoaorls  and  universal  tlaa  but  tha  pragraa  only 
knows  tha  differences  far  the  years  197S*  197S-1979.  If  aarth  tide 
corrections  waaa  ta  ba  eeaputed  for  observations  aada  during  any  other 
year*  It  would  require  appropriate  addition  af  ceding  In  tha  MAIN 
pragraa.  tinea  tha  currant  difference  Is  alaast  ana  alnute*  by  net 
Including  this  difference  could  result  In  a  systeeatle  error  of 
0.8  Vtgal  [Melkklnen*  1978].  Since  the  tlae  far  aaeh  observation  used 
was  recorded  ta  the  nearest  alnute*  It  can  ba  expected  that  the 
accuracy  af  tha  theoretical  aarth  tide  eeaputed  would  ba  an  tha  order 


107 


of  0.5  ygal. 

S.S.2  Correction  ±2£  Ull  Hoieht  s±  Instrument 

Zn  order  to  essure  thet  the  observations  at  a  site  refer  to  the  same 
point*  any  correction  for  the  height  of  Instrument  must  be  made.  The 
correction  for  the  height  of  instruawnt  above  the  station  is  made  by 
assuming  a  local  gradient  of  gravity  for  the  station.  Since  this 
Information  is  not  generally  available*  the  value  of  normal  gravity* 
•S(ti  eStvSs  (10  ^ygal/m),  Is  used.  Thus  the  correction  to  be  made  is 
equal  holoht  of  Instrument  in  meters  times  — 300.4  which  gives  the 
correction  needed  in  units  of  (igal.  This  correction  could  be  the 
source  error  if  there  were  a  large  varietion  in  the  height  of  an 
instrument  at  a  station  because  the  local  gravity  gradient  can  easily 
differ  from  the  normal  gravity  gradient  by  20  to  30  per  cent  or  more 
[Marson  and  Alasia*  1900].  This  error  could  be  on  the  order  of  1  ygal 
per  cm  of  elevation  difference. 

*•>•3  inatrimnUH  Britt 

Zn  order  to  check  for  the  possibility  of  drift  in  the  *6*  gravity 
meter*  a  modal  must  be  assumed  which  will  represent  this  behavior. 

Since  the  drift  is  believed  te  be  the  accumulative  effect  of  a  number 
of  small  tares  [Ourris*  1900*  private  communication],  the  hypothesis 
that  these  tares  occur  uniformly  could  be  adopted  which  could  be 
represented  by  a  simple  linear  model  such  as 


10S 

k(t  -  T  )  ts.s> 

o 

ihtra 

k  —  the  drift  rat*  of  tho  instrument, 

t  —  tho  opoeh  of  tho  observati on, 

T  —  bom*  initial  opoeh. 

o 

However,  tho  appropriatonoas  of  including  such  a  tern  is 
luostionablo  bocauso  of  tho  seemly  unpredictable  and  erratic  behavior 
if  tho  tares.  Snail  taros  can  occur  whenever  there  is  any  change  in 
:he  operating  condition  of  tho  gravity  aster,  such  as  would  be  caused 
ly  vibrations  during  transportation,  changes  in  the  aabient  temperature 
ihich  night  causes  tho  hoator  in  the  instruaont  to  cycle  on  and  off 
'opoatodly,  or  changos  in  the  ataospheric  pressure. 

A  number  of  least  squares  adjustments  were  performed  using 
ibservations  aado  along  tho  United  States  Mid-Continent  Calibration 
in*  and  a  linear  model  siallar  to  that  expressed  in  equation  (5.2) 
ihich  incorporated  a  drift  tern  for  each  gravity  motor.  A  number  of 
djustnents  were  made  involving  different  combinations  of  tho  various 
ypos  of  drift  rate  terns.  The  drift  terns  tried  involved  solving  for 
linear  drift  rat*  for  each  instrument,  a  linear  drift  rate  for  each 

] 

1 

rip,  and  oven  a  linear  drift  rate  involving  only  ties  between  the  same 
tations.  Tho  latter  drift  rate  term  could  bo  referred  to  as  night 
rift  because  the  ties  were  made  between  the  same  stations  which. 


enerally,  occurred  as  a  result  of  overnight  stops. 
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For  example.  whan  a  linear  drift  rata  for  aaeh  instrument  was 
Included  in  the  modal,  there  was  no  evidence  of  any  linear  drift  rate 
for  any  Instrument  Mhich  had  observations  made  over  a  period  of  a  month 
or  more.  Observations  for  instruments  Mhich  Involved  a  more  limited 
period  did  indicate  a  possible  linear  drift,  but  these  drifts  Mere 
discarded  because  there  Mere  too  few  observations  or  the  time  period  of 
the  observations  Mas  too  short  to  draM  tha  conclusion  that  the  drift 
Mas  the  behavior  of  the  gravity  mater.  The  adjustments  performed 
indicated  that  the  linear  drift  rates  of  the  type  mentioned  Mere  not 
representati ve  of  the  characteristics  of  any  of  the  gravity  meters 
used  in  this  study. 

This  fact  became  more  apparent  Mhan  linear  drift  rates  for  each  trip 
Mere  Included  in  the  model.  The  linear  drift  rates  determined  for  a 
gravity  meter  Mere  inconsistent  both  in  magnitude  and  direction.  This 
Mould  lead  one  to  conclude  that  if  there  Mere  a  linear  drift  rate  for  a 
gravity  meter,  it  is  not  of  an  Instrumental  nature  but  possibly  a 
function  of  the  instruments*  handling,  mode  of  transportati on  or  some 
other  unknown  cause. 

No  conclusion  could  be  dramn  concerning  the  existence  of  a  linear 
night  drift  rate  due  to  the  lack  of  a  suff leant  number  of  'night* 
observations  for  any  gravity  meter. 

It  is  worth  mentioning  that  although  the  gravity  meter  observations 
used  in  this  study  indicated  no  predictable  linear  drift  rate  behavior, 
this  does  not  mean  that  a  linear  drift  rate  for  the  gravity  meters  does 
not  exist.  Zt  is  quite  possible  that  if  all  the  observations  Mere  made 


1X0 


so  the  tios  formed  ladder  loops  instead  of  the  Mdiffcd  ladder  and  line 
loops  which  typifies  tha  observations  In  this  study,  the  cencluslon 
eight  be  different. 

The  apparent  erratic  behavior  of  the  linear  drift  rate  Is 
deeonstrated  latter  on  In  section  (.1. 

To  deeostrate  tha  apparent  erratic  behavior  of  a  possible  linear 
drift,  an  adjustment  was  perforeed  using  the  saee  nodal  as  used  In 
ADJUSTMENT  B,  see  section  4.S,  with  T1  and  T2  both  having  the  value  of 
24  hours.  The  residuals  obtained  were  then  plotted  against  tine,  fart 
of  the  plot  for  gravity  eater  *0-01*  can  be  seen  In  Figure  22.  As 
aentioned.  Figure  22  shews  no  Indication  of  any  uniform  linear 
instruaental  drift  rate. 

B.  .4  ffMftMUmil  Mrttll  ttliJMI  YlllUl  lit  PIUlQll 

The  values  in  ailllgal  must  be  used  as  the  observables  In  an 
adjustnent  whenever  It  Is  unreasonable  to  use  counter  readings  as  the 
observables.  This  occurs  whan  a  gravity  eater  has  been  used  over  only 
a  limited  range  of  gravity  values,  generally  less  than  400  egal.  In 
order  to  use  the  counter  readings  as  the  observables  requires  solving 
for  a  long  wave  length  sinusoidal  term  whose  period  Is  always  larger 
than  4000  counter  units  whleh  is  equivalent  to  approximately  4000  egal. 
Trying  to  solve  for  the  coefficients  of  the  sinusoidal  term  using  a 
llm'iis*'  range  of  data  would  lesd  to  poorly  determined  values  of  the 
coefficients  which  could  be  reflected  in  incorrect  values  for  ether 
parameters  of  the  adjustment,  such  as,  linear  scale  factor  teres  and 


gravity  station  valuos.  Thor of or a,  whenever  a  gravity  Motor  Is  used 
which  has  observations  spanning  a  llMltod  gravity  range,  the 
appropriate  Model  to  be  used  weuld  be  ene  using  the  value  In  mIIII  gal 


The  actual  observatlen  Made  with  the  gravity  Meter  Is  the  reading 
Made  of  the  position  of  counter  and  the  dial  after  the  Instrueent  has 
been  nulled.  The  reading  la  generally  recorded  In  counter  units.  The 
readings  In  counter  units  are  to  be  related  to  a  value  In  Mllllgal 
through  the  use  of  the  Calibration  Table  1.  If  one  believes  that  the 
gravity  Meter's  behavior  reawlns  constant  with  ties  and  Is  reflected 
via  the  Calibration  Table  1.  then  the  value  In  Mllllgal  obtained  free 
the  Calibration  Table  1  would  be  correct  except  for  a  possible 
correction  to  the  Initial  offset  and  to  a  scale  factor.  These 
correction  teras  could  be  considered  as  part  of  a  required 
transforwatlon  that  Must  be  applied  when  using  the  Calibration  Table  1. 
For  a  sleple  transforwetl on.  ene  could  consider  a  linear  transforwatlon 
which  would  laply 

z  ■  ly  +  R  (5.*) 

where 

R  -  somo  offset  value  te  be  applied  to  Calibration  Table  1. 

1  —  scale  factor  required  to  be  applied  to  Calibration  Table  1. 

y  —  value  In  Mllllgal  frea  Calibration  Table  1. 

Then  by  rewriting  equation  (2.1)  where  f(0>  »  z,  C  * 
end  z  a  ly  ♦  2  results  in 
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z  +  C  +  S-  G-  0 


(5.10) 


where 

—  earth  tide  correction. 

Cj  •  height  of  instruaont  correction. 

Cj  —  other  car r act Ions  dua  ta  environaental  factors. 

If  ana  assuaes  that  station  ft.  is  observed  at  tiae  1  and  station  ft. 

1  J 

is  obaarvad  at  tiae  j,  than  by  differencing  tha  two  consecutive 
aquations  of  tha  fora  (ft. 10)  for  the  two  different  stations,  one 
obtains 


1(V|  -  Yj)  +  c,  -  Cj  -  G,  ♦  Gj  -  0 


(8.11) 


Since  it  is  asauaad  that  tha  Calibration  Table  l  supplied  is  relatively 
good,  the  value  of  tha  scale  factor  1  will  be  vary  close  to  one.  and  we 
can  rewrite  1  as  1  *  1*  giving 


Y|  -yJ  *  ,,(yl  *V  +C,  -  Cj  -G,  ♦Gj  -0 


(5.12) 


Tha  siailar  tranafaraatian  needed  far  tha  Calibration  Table  1  could 
be  expressed  as  a  higher  order  polyneaial  transforaation  giving 


»l  *  yJ  Vj)  ♦  Cj  *  Cj  -  Gj  +  Gj  ■  ° 


(8.13) 


where 


1 •  —  tha  scale  factor  corresponding  to  the  i— th  degree. 
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Than  If  two  auceeaaive  ebaer vat Iona  war*  differenced,  lattlng 


•y,  -  yj> 


y I j ,  tha  reaulting  equation  would  be  the  baaie 


■adal  far  gravity  differeneea.  Zt  ia  apparent  that  if  n>  1,  gravity 


diffaraneaa  ean  na  langar  ba  ui 


il  gi  van 


by  equation  (8. IS)  ainea  Cy”  —  y”  >  a  (y^  —  y^)0  anly  whan  n  *  1.  Par 
exawple,  if  n  a  2.  tha  right  hand  aida  af  tha  pravioua  aquation  baeowaa 


-  *y}yj  ♦  y*  ii 

eaaa  whan  y.  *  y.. 


y|  —  yj  axeapt  in  tha  trivial 


5.S.B 


Ofidtl  lfi£  Count ar  Eliding 


Tha  funetianal  relationahip  af  tha  obeorvable,  f(0>,  aa  givan  in 
aquation  (2.1)  ean  ba  axpraaaad  approximately  by  aquation  (5.7)  with 
n  a  1  and  k  a  1.  Uaing  tha  information  in  aquation  (5.7),  aquation 
(2.1)  ean  bo  rowrittan  aa 


n  ,  k 
r  i  jl  .  r 


l  1  ,d'  +  y  A.  cos  (  +  u>,  )  +  S  +  C-  G*»0 

-1  1  j-t  j  TJ  j 


(5.14) 


and  aiailarly  aquation  (2.2)  ean  ba  rowrittan  aa 


i i.(dj -dh ♦  i  i (-i)"*'* , «>( ^♦u,) ♦ 

lit  1  '  Z  j£|  n-l  J  Ti  •> 


k  2 


(5.15) 


den  (5.15) 


Cv  -  C2  -  G,  +  G2  -  0 


leaaaa  tha  baaie  aathaaatieal  aodel  for  tha 


rQf  gravity  aatar  whi 


Th*  us*  of  aquation  (8.18)  laplles  there  oust  b*  a  sufficient 
distribution  of  th*  observations  oado  with  a  gravity  Motor  to  solvo  for 
tho  eoofflelonts  of  tho  sinusoid.  Hhat  actually  constitutes  a 
sufficient  distribution  Is  difficult  to  state.  But  cloarly,  certain 


If  all  th*  observations  were  oado  at  stations  having  approx1»at*ly  tho 
saao  gravity  value,  then  th*  aedel  expressed  by  equation  (8.18)  Mould 
b*  Inappropriate. 
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XIIZZJ1S  a t  IHi  mathihatical  modil» 

4.1  tolut t on  of  tho  ttfttfll 

The  watheMatleal  — 4»1  uuA  for  th«  gravity  adjustwont  program  can 
Im  expressed  by  F(X*D  a  •  *  where  F  la  a  vaetor  of  nonlintar  functions 
of  thM  paraweters  X  and  tho  oPsorvablos  L.  Tho  problaw  Is  to  find  tho 
vectors  X  and  L  such  that  vSv  Is  a  nlnlnun  subjoet  to  tho  constraint, 
F(X,L)  a  o.  Ivon  though  tho  weight  Matrix  for  tho  obsorvat 1 ons ,  F*  Is 
not  a  unit  Matrix*  tho  Method  used  to  solve  this  problow  Is  eoaaaenly 
reforrod  to  as  "least  squares"  ovon  though  tho  More  appropriate  naoe 
would  be  "weighted  least  squares". 

The  solution  of  this  type  of  nonlinear  least  squares  problaw  Is  by  a 
Howton-Causs  Iteration*  which  expresses  the  nonlinear  function  F  as  a 
Taylor's  series  about  sows  initial  values  of  the  peraweters  and 
observables,  by  ignoring  the  second  and  higher  order  tense  of  the 
Taylor's  series*  It  redeeee  to  solving  a  linear  least  squares  problaa. 
If  the  function  F  is  net  satisfied  by  the  adjusted  parowotors  and 
observables  obtained*  the  solution  Is  Iterated  by  using  the  adjusted 
paraweters  and  observables  as  the  new  points  of  expansion  for  the 
Taylor's  series  {Fops*  1*7*1. 


Tho  aultlvarlate  Taylor’s  series  expansion  of  F(X,L>  about  X  ,L 

o  o 

with  the  condition  P(X.L)  a  |  boeoaos 


(x-x)  ♦ 

V-o 


-  L0>  ♦  F(X0.L0)  -  0 

Xo-Lo 


(6.1) 


VLo* 


To  slapllfy  tho  notation,  lot  X  ■  J  -  X  ,  A  •  SF/Slf  I  and 

0  I  X,  L 

«  a  o  o 

B  “  3F/3L  J  ^  Tho  residuals,  V,  of  tho  obsorvablos  aro  given  by 

V  -  L  -  L.  where  U  Is  tho  adjusted  value  of  tho  observable  and  L.  Is 
D  b 


tho  observed  value.  This  relation  can  be  expressed  as 


L  -  Lc  -  »  ♦  Lb  -  Lc 


(6.2) 


which  allows  equation  (6.1)  to  be  written  as 


w  (I '  ♦§  »  ♦§  (S  -  Lo)  ♦  F<VLo>  - 0 

xo,Lo  VL.  Xo'Lo 


(6.3) 


Using  the  notation  previously  aentloned  yields 


AX  f-  BV  +  B(Lb  -  L0)  F(X0,L0)  -  0 


(6.4) 


and  by  letting  H  a  F<XQ,L0>  ♦  XL^—  LQ),  equation  (6.4)  boeoaos 


AX  ♦  BV  +  W  -  0 


(6.5) 


With  T.  being  the  covariance  aatrlx  of  tho  obsorvablos  and  O*  being 
Lb 

soae  constant,  the  weight  aatrlx,  F,  for  tho  obsorvablos  Is  related  to 


(«.«> 


The  derivation  of  tho  standard  linear  least  squares  problem  Is  well 
documented  (Hamilton,  lH4t  Mikhail  and  Ackermann,  1976]  and  yields 


X  -  -(At(BP“tBt)"1A)"1At(BP“,Bt)"1W 


<6.7> 


V  -  -  P’1Bt(BP*1Bt)"1(AX  ♦  W) 


<6.*> 


with  the  best  estimate  of  the  vectors  X  and  L  being  given  by 


X  «  X  ♦  X 
o  o 


(6.9) 


ro  -  Lb  +  5 


(ft. 10) 


of  the  pa 


satisfy  the  function  P,  the  solution  needs  to  be  Iterated  with  the 


point  of  expansion  of  the  Taylor's  series  new  being  This 


now  evaluated  at  the  new  point  of  expansion  and 


the  M  matrix  Is  recomputed.  If  V  denotes  the  vector  of  residuals  for 


the  observables  from  the  previews  iteration. 


w  ■ 


(6.11) 


noting  that  for  the  seroth  Iteration  V  will  be  a  sere  vector. 


The  Iterative  process  continues  until  the  mathematical  model  Is 


satisfied.  Due  to  round-off  error  In  the  computer's  representation  of 


values  and  coaputatlon  of  tha  Matrix  operations,  the  chances  are  the 
■odel  will  never  truly  be  satisfied.  However,  It  can  generally  be 
satisfied  so  F(X,L)<e  ,  where  e  Is  a  vector  as  saall  as  desired  within 
the  Halts  of  the  round-off  error  of  the  coaputatlons.  Therefore,  the 
Iteration  should  continue  until  F(X,U)<e  Is  satisfied.  The  nuaber  of 
Iterations  required  will  depend  on  how  close  the  Initial  paraaeter 
values  and  observables  are  to  the  values  required  to  satisfy 
FCX,D  *  0. 

The  aatheaatleal  aodel  F(X,L>  a  B  Is  really  the  cowposlte  of  two 
functions,  F^CX.Lj)  a  0  and  F2<X,L2>  a  0,  where  the  vector  of 
paraaeters,  X,  Is  the  saae  for  F^  and  F^  Further,  L ^  Is  the  vector  of 
observables  froa  gravity  asters  and  L2  Is  a  vector  of  derived  absolute 
gravity  station  values  froa  either  peraanent  or  portable  absolute 
gravity  aeasurlng  devices.  This  laplles  that  there  are  "observed** 
gravity  station  values.  These  values  are  "observed"  In  the  sense  that 
a  variance  for  the  stat1on*s  gravity  value  exists.  This  type  of  aodel 
has  been  referred  to  as  the  "eoablned  observation  and  paraaeter  aodel 
with  weighted  paraaeters"  (Uotlla,  1*67]. 

If  one  partitions  the  A,  t,  F  and  H  Matrices  Into  two  parts;  part  1 

for  F  <X,L . >  ■  0  and  part  2  for  F  <X,L>  ■  0,  and  letting  M  s  IF-1** 

11  2  2 

where  the  observed  values  ef  L.  and  L  are  not  correlated  then 
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H1  0 


0  M„ 


B1  0 


0  B« 


0 


0  P. 


B?  0 


0  B 


(ft. 12) 


St  can  than  be  shewn  that  equations  (ft. 7),  (i.t>  and  (6.10)  can  be 

expressed  as 


mv.  ofrMr’fMV  °riM 


0  M, 


0 


(6.13) 


M  [M  “lf'f  °1h  0 ff[M x * K 

»2  0  M  0  »'  0  H2  1  A2  w 


(ft. 14) 


•  SI  W  N 

which  can  be  -further  slnpll-fled  te 


(ft. IS) 


X  -  -(AfM“lA1  +  aJm^Aj)-1  (aJm^Wj  +  AjM"'W2) 


(ft. lft) 


p^bJh'^a^x  +  w,) 
P21b2M21(A2X  +  W2^ 


(ft. 17) 


L°‘  r. 


S  * v, 

Lh  +  V2 

L  b2  2J 


(ft. IS) 


2  t 

Xn  addition.  It  ean  be  shewn  that  a‘  a  V  pv/D.F.  where  D.F.  represents 
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th«  degrees  of  freedom  of  the  adjustment  which  Is  equal  to  the  number 
of  equations  minus  the  number  of  parameters  In  the  medal  and 

V*PV  -  (AX  ♦  VOV1!!  («.!*> 


or 

VCPV  -  vfp1V1  ♦  v|p2V2  -  (x'aJ  ♦  +  (X*A*  +  W2)M21W2  U.20) 

Additional  examination  of  the  structure  of  the  A,  B,  P  and  H 
matrices  reveals  that  A^.  Bj,  t  ^  and  M  j  matrices  can  be  further 
partitioned  by  trips  and  Aj.  »2'  f  2  *nd  M2  Mtr1ca*  can  be  partitioned 
by  each  set  of  absolute  gravity  observation  equations.  Applying  this 
Information,  equations  <6. It),  ((.17)  and  (A. 20)  can  be  written  as 


■{,£,  j£,AuH->u}  {,*,  j^AuM.Jw.j} 


V*PV 


plKiHn<An*  +  w„> 

+  V 

P21B21H21*A21X  +  W21* 
P2tB2tM2t^A2tX  +  W2t*J 
2  rl 


(6.21) 


<*. 22) 


(6.2S) 


where  s  *  r^  and  t  a  r2  with  r^  being  the  nusiber  of  gravity  meter  trips 
made  and  r2  being  the  number  ef  sets  ef  abselute  gravity  observation 
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iquitieni. 

The  B^and  ^matrices  have  a  vary  wall  defined  pattern  that*  can  ba 
exploited  In  forming  the  Mj  matrices.  Assuming  that  all  observations 
aada  with  a  particular  gravity  motor  over  a  trip  have  the  same 
variance*  a j ,  where  1  denotes  the  trip  number*  then  the  covariance 
aatrlx  of  the  gravity  meter  observation  for  trip  1  Is  given  by 

P11  “  al  I  IS. 24) 


vhere  I  Is  and  nxn  Identity  aatrlx  and  n  Is  the  number  of  observations 
in  the  trip. 

The  elements  of  the  •  matrix  for  the  gravity  meter  observations  for 
trip  1  are  given  by 


3F 

w. 

3F 

1 < 1 <n— 1  and  j«l+1 

(6.25) 

0 

j>i»1  and  j<} 

there  lj  is  the  j—th  observation  in  trip  1.  The  patterned  structure  of 
:he  B  matrices  for  the  gravity  meter  observations  is  simlliar  to 


s 


MIDI 
0  «  «  0  0 
0  0  »  *  0 


0  0 
0  0 
0  0 


(6.26) 


0  0  0  0  0  ...»  0 

0  0  0  0  0  ...»  « 


here  »  represents  a  non-zero  element  and  0  represents  a  zero  element 


of  the  natrlx  B. 


Since 


M11  "  B1 1 P1 1B1 1 


(6.27) 


by  substituting  aquation  (6.24)  and  (6.2S)  into  (6.27)  gives 


7H11  -81lBfl 


(6.26) 


with  the  value  of  the  elements  of  M  .  given  by 


mii  “  bH  +  bl+1,|+1 


1 < I <n- 1 


mIj  "  mjl  "  “bIj 


j"i+1  and  1 < I <n— 1 


(6.29) 


"jj  "  mji  "  ° 


j>?+1  and  1 < i <n— 1 


There  are  times  when  the  B  matrix  has  its  elements  satisfying 


bii  •-‘u 


j«l+1  and  1 < i <n-l 


(6. SO) 


This  occurs  during  the  initial  adjustment  iteration  when  the  parameters 
other  than  the  gravity  station  values  have  an  initial  value  of  zero. 

In  Hhich  case  the  inverse  of  matrix  can  be  obtained  directly 
[Oergory  and  Karney,  1969,  pp  4B— 46]  from  the  expression 


°i  *  s4t  c 


(6.31) 


Mhere  n  is  the  order  of  the  matrix  M, , ,  k  is  the  value  of  a  non— zero 


eleaient  of  the  B  matrix  squared  and  the  elements  of  C  are  given  by 
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C.j  -  Cjj  -  l(n  -  j  +1)  i<j  “•**> 

c.2  computer  freanaw 

The  computer  program  developed  to  handle  the  adjustment  was  written 
In  PL/1  and  designed  to  run  en  the  Amdahl  470  V/6— XX  computer  using  a 
virtual  cere  region  of  2043  Kbytes.  The  reason  the  program  was  written 
In  PL/1  was  to  allow  for  dynamic  allocation  of  arrays  and  the  ease  of 
handling  output  formatting. 

The  PL/1  program  written  to  perform  the  gravity  network  adjustment 
was  designed  to  meet  the  following  requirements: 

a)  the  program  should  be  able  to  handle  any  number  of  unknowns  less 
than  32767. 

b)  the  program  should  be  able  to  determine  what  stations  are  in  the 
network,  what  instruments  are  being  used  and  should  be  able  to 
break  the  observation  set  into  trips  with  the  capability  of 
handling  instrumental  tares. 

c>  the  program  should  be  able  to  use  different  models  for  each 
instrument. 

# 

d)  the  program  should  allow  for  full  variance-covariance  weighting 
of  any  parameter. 

e>  the  program  should  be  able  to  do  post  adjustment  analysis. 

f)  the  program  should  be  able  to  iterate  on  the  solution. 

g>  the  program  should  havo  as  few  fixed  values  and  limits  as 


possible. 


For  tho  Most  part*  all  tho  above  requirements  wore  satisfied  except 
for  one  Minor  restriction.  Sue  to  a  PL/1  limitation  which  restricts 
the  sice  of  an  array  subscript  to  a  Maxi  mum  of  *27*7,  the  present  liMlt 
on  the  nuMber  o-f  observations  allowed  In  any  trip  is  101  observations. 
This  restriction  Is  necessary  because  when  the  weight  Matrix,  H  \  for 
the  equations  o-f  a  trip  Is  full  and  It  Is  required  that  all  eleMents  of 
the  Matrix  to  be  In  core  at  the  same  tloe.  This  llMltation  has  not 
caused  any  problees  In  this  study. 

The  prograe  can  be  thought  of  as  consisting  of  three  sections: 
Pre-Processor,  Adjustment  and  Post— Analysis.  The  options  for  each 
section  are  controlled  by  values  assigned  to  a  fixed  binary  external 
array  of  ISO  eleMents. 

*•2.1  rrt-f»-ggtaagr  HfiUgfi 

•y  prior  analysis,  the  Input  data  set  of  the  observations  had  all 
detectable  blunders  rewoved  and  any  large  suspected  tares  flagged.  The 
Method  used  to  detect  the  tares  and  blunders  Involved  determining  the 
Mean  gravity  difference  between  all  of  the  stations  which  were  tied. 

The  gravity  difference  for  each  gravity  meter  tie  was  determined  using 
the  the  values  In  Milllgal  from  the  gravity  Meter's  Calibration  Table  1 
which  had  been  Interpolated  from  their  counter  reading  observations. 

Any  large  difference  which  was  greater  than  100  ygal  from  the  Mean 
gravity  difference  between  the  two  stations  tied  was  flagged  so  the 
observations  Involved  could  be  Investigated.  In  the  majority  of  the 
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eases,  the  flagged  observations  had  natations  an  the  original 
ebsarvatien  sheets  which  explained  why  a  large  difference  might  has 
been  encountered.  The  most  common  notations  were  "beam  vibrations'*  and 
"meter  off— heat".  Using  this  information,  it  was  possible  to  determine 
which  observations  should  be  treated  as  a  tare  and  which  ones  were 
blunders.  A  blunder  is  indicated  when  two  consecutive  gravity  meter 
ties  are  flagged  as  suspect.  The  common  observation  involved  in  the 
two  ties  would  bo  labelled  as  a  blunder  and  removed  from  the  data  set. 

The  purpose  of  this  section  is  to  determine  the  make-up  of  the  input 
data  set  of  observations.  It  determines  how  many  different  stations 
and  instruments  there  are  and  breaks  the  observations  into  trips  based 
on  a  number  of  criteria.  If  it  has  been  determined  that  a  tare  has 
occurred  between  observations,  then  the  current  trip  ends  with  the  last 
observation  before  the  tare  and  a  new  trip  begins  with  the  first 
observation  after  the  tare,  based  on  the  time  interval  between 
observations,  a  new  trip  is  begun  if  the  time  interval  exceeds  a 
maximum  time  specified  which  depends  on  whether  the  time  interval  is 
between  the  same  station  or  different  stations.  After  a  tentative  trip 
has  been  identified,  it  is  checked  to  see  that  there  are  at  least  a 
minimum  number  of  observations  in  the  trip.  Zf  not,  then  the 
observations  in  that  trip  are  deleted  and  are  not  included  in  the 
adjustment . 

After  the  observations  to  be  used  have  been  identified  and  broken 
into  trips,  the  stations  are  assigned  parameter  numbers  based  on 
ascending  order  of  their  assigned  identification  code.  The 
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Identification  cod*  Is  either  the  IBS  (International  Gravity  Bureau) 
lumber r  if  the  station  was  pert  of  the  XGBN  71,  or  some  other  arbitrary 
assigned  code  for  any  other  station.  The  instruments  are  assigned  a 


Any  instrument  requires  a  model  different  from  the  default  model. 


which  was  the  model  using  the  value  in  milliga 


code,  input,  saved,  and  the  remaining  model  parameters  were  assigned. 
This  is  followed  by  the  assignment  of  initial  values  to  all  the 
parameters  based  on  input  information.  A  final  check  is  made  to  be 
sure  that  there  are  sufficient  redundancies  for  each  instrument  and 
that  all  the  stations  are  inter-connected  by  observational  ties  to 
guard  against  trying  to  solve  a  singular  system. 


B.2.2  MlgatffWrt 


The  purpose  of  this  section  is  to  actually  perform  the  weighted 
least  squares  adjustment.  For  each  Iteration  performed,  all  the 


non— sere  parti  els  with  respect 


computed  along  with  the  medal  miselesures  and  the  full  weight  matrix, 
H*1,  for  each  trip  is  formed.  See  Appendix  A  for  e  detailed  discussion 
of  how  this  is  accomplished.  The  eentributio.*  to  the  normal  equations 
for  each  trip  is  then  computed  by  performing  only  non— sere 
multiplication.  See  Appendix  B  far  a  detailed  discussion  of  how  this 


is  done 
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The  solution  Is  Iterated  until  tho  model  Is  satisfied  or  tho 
Iteration  count  exceeds  tho  maximum  allowed.  Tho  aodol  P(X,L>  Is 
satisfied  for  soao  prodotoralnod  value  of  e»  uhon  tho  condition 
pcx, lx  e  is  fulfilled. 

«.2.S  Post-Analysis  fffiUttn 

Tho  purpose  of  this  section  Is  to  do  a  post— ad justment  analysis  of 
the  results  of  the  adjustment.  Zn  this  section  the  variance  of  each 
residual  Is  computed,  and  the  noraallsed  residuals  are  foraed.  For 
details  on  how  this  Is  done,  see  Appendix  C.  Any  large  normalized 
residual  Is  listed  to  be  checked  for  possible  blunder  or  tare.  The 
aeanlng  of  large  Is  based  on  the  Tau-erlterla  as  outlined  by  Pope 
[1974]  and  depends  on  the  nuaber  of  observations  and  the  significance 
level  selected.  Zn  this  section,  the  root  aean  square  of  the  residuals 
for  each  instrument  is  coaputed.  The  program  also  computes  an  estimate 
of  the  aposteriorl  aceuarcy  of  the  observable  by  trip  and  instrument 
based  on  the  estimated  initial  variance  of  tho  observable,  its 
residual,  and  the  aposteriorl  variance  of  the  residuals. 

Twttna  Ytr.l ma  Haiti* 

Per  the  comparison  of  the  various  possible  models,  three  different 
adjustment  mere  performed  using  the  same  set  of  observation  data.  All 
adjustments  were  made  using  the  same  set  of  absolute  sites  as  control. 
The  estimated  stand  errors  assigned  to  the  absolute  sites  ara  given  in 


Table  S 


1 


1 


umption  that  the 


value  of  gravity  determined  at  a  alt*  using  tha  absoluta  gravity  meter 
davalepsd  by  either  Harsen  and  Alasla  or  Namnond  had  tha  saaa  accuracy 
and  gravity  valuas  did  not  Includa  any  bad  dot oral  nations.  By 
ravlawlng  tha  gravity  stations  In  Tabla  4  and  Table  5,  saa  Section  4.2, 
which  had  Multiple  deteralnatl ons.  It  can  be  seen  that  156E05,  Great 
Falls,  valuas  differ  by  101  Ugsl  snd  that  15211A,  Boston,  valuas  differ 
by  109  ygal.  Zt  follows  that  with  a  probability  of  0.99,  a  single 
deterialnatl on  wade  with  aithar  apparatus  would  11a  within  2.57a  of  the 
Mean  gravity  value  for  that  station.  Xf  it  ware  assumed  that  the  naan 
gravity  value  at  a  station  wars  located  at  tha  middle  of  tha  two 
extrema  gravity  valuas,  than  tha  accuracy  of  a  determination  could  be 
estimated  to  be  batter  than  s  s  ((x2  —  xl)/2)/2.57  where  x2  Is  tha 
largest  gravity  value  determined,  xl  is  the  smallest  gravity  value 
determined,  and  s  Is  the  estimate  of  the  accuracy  of  a  determination. 
Using  this  relationship,  the  accuracy  of  the  absolute  gravity  meter 
determinations  which  could  account  for  the  different  gravity  values 
reported  would  be  approximately  20  ygal  Inetead  of  the  reported  10  ugal 
accuracies.  The  accuracy  of  115V01,  Casper,  was  kept  at  25  ygal 
because  there  was  good  no  reason  to  lower  Its  reported  accuracy. 

The  variance  for  each  gravity  meter  observation  was  assumed  to  be 
the  same  for  all  gravity  meters.  The  variance  assigned  to  the 

2 

observations  made  with  the  gravity  meters  were  0.0004  counter  units  or 
400  ygal2,  depending  en  which  type  of  observable  was  used  in  the  model 
for  each  gravity  meter.  The  above  values  were  based  on  preliminary 
adjustments  which  Indicated  that  no  large  differences  existed  between 
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tha  (posteriori  estimates  of  tho  accuracy  of  the  observations  for  any 
gravity  actor.  A  total  of  4332  gravity  observations  wore  consi dared  of 
which  41  was  rejected  because  of  detected  tares  or  blunders.  The 
reaaining  4434  observations,  aada  with  27  different  gravity  aeters, 
resulted  in  337  different  trips  based  on  the  value  of  24  hours  for  T1 
and  4  hours  for  T2  (see  Section  3.3.2).  The  ainiaua  nuaber  of 
observations  in  any  trip  is  2  while  the  largest  trip  consisted  of  23 
observations.  See  Figure  23  for  a  histograa  of  the  frequency  of  the 
trips. 

The  three  adjustaents  will  be  referred  to  as  ADJUSTMENT 
ADJUSTMENT  *,  and  ADJUSTMENT  £.  All  adjustaent  aodals  included  a 
linear  scale  factor  tara  but  no  drift  teras  for  each  gravity  aeters. 
ADJUSTMENT  A  was  perforaed  using  as  observables  the  observed  counter 
units  transforaed  to  their  values  in  ailligal  by  using  the  factory 
supplied  Calibration  Table  1.  ADJUSTMENT  3  and  C  were  perforaed  using 
as'  the  observable  tha  observed  counter  readings  whenever  the  range  and 
distribution  of  tha  data  justified  their  use.  Otherwise,  the  values  in 
ailligal  were  used.  Seth  ADJUSTMENT  S  and  C  aodals  Included  a  long 
wave  sinusoidal  tera  to  represent  the  Calibration  Table  1  inforaation. 
ADJUSTMENT  S’s  aodel  contained  ne  periodic  screw  error  teras  whereas 
ADJUSTMENT  C*s  aodel  included  periodic  screw  error  teras  having  a 
period  of  1284/17  counter  units  for  selected  gravity  aeters.  All 
gravity  aeters  used  in  the  adjustaents  had  "old*  gear  boxes  installed 
in  them.  The  gravity  aeters  whieh  included  a  periodic  screw  error  tera 
were  selected  based  on  the  nuaber  and  distribution  of  their 
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Table  S  —  List  o-f  absolute  stations  used  In  ADJUSTMENTS  A,  B, 

2SG/GSS 

Code 

Station 

Nano 

Gravity  Value 
in  ^Bgal 

Estimated 

Standard 

Error 

In  ;Ugal 

ba 

U9S04 

McDonald  Obs..TX 

97BB20097 

20 

119S03 

McDonald  Obs.»TX 

978328455 

20 

08150C 

Mlaal ,FL 

979004319 

20 

9 

11924A 

Moll omen, NM 

979139592 

20 

■  -  ;  . 

119C03 

Mt.  Evans ,C0 

979254059 

20 

1 ISCSI 

Trinidad' CO 

979330370 

20 

"SI 

U994H 

Denver 'CO 

97959B272 

20 

1S50SD 

Boulder, CO 

979405498 

20 

» 

1S5V01 

Casper ,MY 

979947244 

25 

12172A 

San  FrancIscO'CA 

979972040 

20 

15221 A 

Boston,MA 

9B037B459 

20 

1S5B0A 

Bismarck ,ND 

980412882 

20 

Kt  —  International  Gravity  Kuroau 

[  '  GSS  —  Gaodatle  Survey  Squadron 


and  C. 
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observatl ons.  Whenever  periodic  toms  wore  Included,  their  periods 

were  held  fixed  and  their  phase  angles  Mere  assumed  to  have  a  variance 

2  2 

of  0.01  radians  or  «32.3  degrees  .  The  periods  Mere  held  fixed  to 
eliminate  the  high  correlation  Mhleh  exists  betMeen  the  period  and  the 
phase  angle  In  the  periodic  term.  The  phase  angles  Mere  Meakly 
constrained  to  reflect  the  accuracy  of  the  estimate  of  the  value  of  the 
phase  angles  determined  from  the  residuals  of  the  observations  from 
ADJUSTMENT  A.  An  estimate  of  the  phase  angles  for  each  gravity  meter 
Mas  determined  by  using  a  least  squares  adjustment  Mhleh  used  only  the 
residuals  associated  Mlth  the  observations  from  ADJUSTMENT  A.  The 
model  used  Mas  a  sinusoidal  tern  having  a  period  of  1208/17  counter 
units.  From  these  adjustments,  the  accuracy  of  the  phase  angle  terms 
mss  almays  less  than  0.1  radians.  All  the  accuracies  for  the  phase 
angles  were  rounded  off  to  0.1  radians  and  used  In  ADJUSTMENT  C. 

The  Initial  estimates  for  the  amplitude,  phase  angle  and  period  of 
the  long  rnave  length  sinusoidal  term  for  the  gravity  meters  Mere  based 
on  the  analysis  of  each  gravity  meter's  Calibration  Table  1  data.  The 
Initial  estimates  for  the  amplitude  and  phase  angle  of  the  periodic 
screw  error  terms  Mere  determined  from  a  least  squares  fit  of  the 
residuals  from  ADJUSTMENT  •  for  each  gravity  meter.  A  summary  of  the 
results  of  the  adjustments  can  be  found  in  Tables  9—12. 

tevloMlng  Table  10,  one  finds  that  the  aposterlorl  estimate  of  the 
accuracy  ef  the  observations  for  each  gravity  meter  decreased  slightly 
from  ADJUSTMENT  A  values  to  ADJUSTMENT  •  except  '0-81'  Mhleh  shoMed  a 
large  change.  The  reason  for  this  change  Is  unknoMn  at  this  time  and 


Tabic  9  -  Summary  of  the  results  of  ADJUSTMENTS  A,  B  and  C. 


ADJUSTMENT 

A 

ADJUSTMENT 

5 

ADJUSTMENT 

C 

Number 

of  Observations 

4444 

4444 

4444 

Total  Number 
of  Parameters 

245 

317 

351 

Number  of 

Gravity  Stations 

254 

254 

254 

Number  of  Seale 
Factor  Terms 

29 

29 

29 

Number  of  Periodic 
Amplitude  Terms 

0 

14 

33 

Number  of  Periodic 
Phase  Angle  Terms 

0 

14 

33 

Number  of 

IS 

24 

45 

Number  of  Trips 

437 

437 

437 

Number  of  Iterations 

r 

V  PV 

1 

5025.319 

2 

3444.095 

2 

3433.410 

Degrees  of  Freedom 

3374 

3354 

3341 

Aposteriori  Variance 
of  Unit  Weight 

1.44942 

1.04520 

1.02774 

Tabic  10  —  Summary  of  aposteriori  estimates  for  the  accuracy  of 

observations  made  with  the  various  gravity  meters  used  in 


ADJUSTMENTS  A,  B.  and  C 

GRAVITY 

ADJUSTMENT 

ADJUSTMENT 

ADJUSTMENT 

METER 

A 

B 

C 

0-10 

19.78 

14.78m 

17.03m 

D— 17 

17.04 

13.74m 

13.73m 

D— 43 

30.97 

25.85m 

25.47m 

0—44 

14.49 

14.38 

15.44 

0—47 

15.11 

14.84 

15.40 

0-50 

29.35 

29.29 

29.32 

0-48 

31.97 

29.71m 

28.52m 

0-81 

31.21 

18.58m 

18.22m 

0-81a 

18.72 

17.93m 

17.47m 

C-103 

22.03 

20.49m 

20.20m 

0-111 

22.44 

20.81m 

19.54m 

G— 113 

21.94 

21.70 

21.81 

0-115 

19.42 

17.92m 

17.24m 

0-1 15b 

21.22 

19.59m 

19.08m 

0-123 

21.45 

22.42 

22.04 

0-125 

25.47 

23.82M 

22.41m 

G— 130 

14.82 

14.84 

14.87 

0-131 

18.44 

17.42M 

17.41m 

0-140 

18.13 

18.22 

18.18 

0-142 

23.83 

21.54M 

21.04m 

0-157 

23.47 

20.29m 

20.09m 

0— 157c 

22.74 

21 .57m 

21.33m 

0-175 

21.74 

21.74 

21.75 

0-174 

29.14 

27.94 

28.01 

0-191 

18.28 

18.04 

18.19 

0-220 

18.09 

17.05M 

16.42m 

0-253 

24.54 

23.72m 

23.62m 

0-248 

21.94 

20.71m 

20.01m 

0-249 

18.58 

17.47m 

15 .60m 

Average 

Oravity 

for  All 

Meters  22.71 

19.73 

19.24 

a  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  18  October  1977. 

b  —  Calibration  Table  1  changed  due  to  replacement  of  long  lever  on 
27  October  1977. 

c  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  30  August  1977. 

*  —  indicates  that  the  standard 
stanard  errors  are  in^gal. 


error  is  in  dial  units,  all  others 


Table  IX  —  Summary  of  the  amplitudes  of. the  long  wave  sinusoidal  terms 
for  various  gravity  meters  used  in  ADJUSTMENTS  B  and  C. 


Gravi ty 

Initial 

Adjustment 

Adjustment 

Fixed 

Mater 

Value 

•  Value 

C  Value 

Period 

G— 10 

1024 

1576 

1584 

6000 

G-68 

3742 

3617 

3507 

11000 

G— 81 

5502 

5339 

5347 

10000 

G— 81a 

417 

382 

388 

5000 

G— 1 03 

1028 

984 

1010 

7000 

6—111 

2943 

2851 

2858 

9000 

G— 11S 

120 

110 

103 

3000 

G— 115b 

300 

267 

266 

5000 

G— 125 

2174 

1860 

2182 

10000 

G— 131 

410 

367 

373 

5000 

G-142 

586 

551 

563 

5000 

G— 157 

2760 

2646 

2638 

10000 

G— 157c 

1827 

1716 

1757 

8000 

G— 220 

2158 

2J41 

2075 

8000 

G-268 

4874 

4574 

4527 

11000 

G— 209 

685 

758 

655 

7000 

ill  values  are  in  units  of^jgal  except  the  periods  which  are  in  counter 
mi  ts. 


i 

» 

! 


Calibration  Table  1 
on  IS  October  1977. 
Calibration  Table  1 
27  October  1977. 
Calibration  Table  1 
on  30  August  1977. 


changed  due  to 
changed  due  to 
changed  due  to 


addition  of 
replacement 
addition  of 


electronic  readout 
of  long  lever  on 
electronic  readout 
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Table  12  —  Summary  of  tha  amplitudes  of  the  periodic  screw  error  terms 
for  various  gravity  meters  used  in  ADJUSTMENT  C. 


Qravity 

Number  of 

Initial 

Adjusted 

Meter 

Observations 

Value 

Value 

0-10 

164 

8 

7.74 

(3.99) 

G— 68 

61 

19 

18.54 

(6.12) 

0-81 

427 

9 

8.69 

(2.09) 

0-81  a 

171 

3 

2.40 

(4.56) 

0-103 

125 

7 

8.28 

(3.80) 

0-111 

427 

18 

17.22 

(2.45) 

0-115 

369 

7 

6.95 

(2.54) 

0-1 15b 

300 

5 

4.81 

(2.19) 

0-125 

140 

13 

12.31 

(3.10) 

0-131 

350 

4 

4.43 

(2.39) 

0-142 

77 

9 

7.19 

(3.87) 

0-157 

435 

4 

4.19 

(2.54) 

O— 157c 

172 

9 

8.17 

(3.84) 

0-220 

266 

9 

7.35 

(2.95) 

G— 253 

114 

3 

5.40 

(3.56) 

0-268 

237 

9 

10.06 

(2.79) 

0-269 

147 

16 

18.52 

(3.50) 

All  values  are  in  units  of  yifgal  and  the  values  in  parentheses  represent 
the  estimated  aposteriori  standard  error. 

a  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  18  October  1977. 

b  —  Calibration  Table  1  changed  due  to  replacement  of  long  lever  on 
27  October  1977. 

c  —  Calibration  Table  1  changed  due  to  addition  of  electronic  readout 
on  30  August  1977. 


s 
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was  unexpected. 

Tha  valua  of  the  aposteriori  Htiaata  of  tho  accuracies  of  tha 
obsarvatlons  found  in  Table  10  ia  baaad  on  tha  following  aaaunptions. 
Firit.  tha  estimated  varianoa  of  tha  adjusted  observations  and  thoir 
true  variance  for  each  gravity  meter  are  nearly  equal  and  will  bo 
assumed  to  be  equal.  Secondly,  tha  normalized  residuals  for  each 
gravity  aatar  are  assuaad  to  be  neraally  distributed  and  are  effected 
vary  little  by  changes  in  their  assuaad  observational  accuracy.  It  is 
known  that  for  a  least  square  adjustment,  the  variance  of  the  adjusted 
observations  is  a  function  of  the  variance  of  the  observation  and  the 
variance  of  the  residual  for  the  observation.  The  variance  of  an 
adjusted  observation  is  given  by 


(ft. 33) 


where 

O*  —  variance  of  adjusted  observations, 

3 

of  —  variance  of  the  observations, 

D 

aj  ~  variance  of  the  residuals. 

Vi 

A  sinilar  expression  exists  for  the  estimate  of  the  variance  of  an 
adjusted  observation.  It  is  given  in  the  following  equation  where  A 
eeans  the  values  are  estimated  quantities  derived  free  the  adjustment. 


(ft. 34) 


If  it  were  assumed  that  the  true  and  estimated  variance  of  adjusted 
observations  were  equal,  then  the  following  expression  is  obtained  by 
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Since  that**  was  no  reliable  information  about  the  accuracy  of  tha 
observations  Made  with  tha  LaCoste  4  Romberg  gravity  meters,  an 
accuracy  of  an  observation  had  to  be  estimated.  The  value  selected  mas 
0.02  counter  units  which  is  approximately  equivalent  to  20  ygal.  This 
value  was  based  on  the  results  from  preliminary  adjustments  which 
tested  various  estimates  for  the  observational  accuracy.  The 
aposteriori  estimate  of  the  accuracy  of  the  observations  from  these 
adjustments  indicated  accuracies  around  0.02  counter  units  for  each 
gravity  meter  would  be  appropriate.  Using  the  relationship  expressed 
in  equation  <6. SO),  the  values  found  in  Table  10  were  computed  after 
the  final  adjustments  had  been  don*.  These  values  indicate  that  the 
gravity  meters  used  in  this  study  had  approximately  the  same 
observational  accuracy.  It  is  in  the  range  of  0.010  to  0.025  counter 
units  which  is  equivalent  to  an  accuracy  of  about  15  to  25  ygal. 
counter  units. 

To  determine  if  there  is  a  significant  difference  between  the 
adjustment  medel  used,  it  would  be  convenient  if  a  statistical  test 
existed  to  test  the  significance  ef  the  nonlinear  models  used. 
Unfortunately,  the  tests  proposed  for  testing  nonlinear  hypothoses  are 
often  too  complex  to  bo  used  in  practice  and  depend  on  the  nonlinearity 
of  the  problem  (Hamilton,  1944,  p  147).  2n  practice,  tests  for  the 
linear  ease  are  performed.  Considering  that  the  tests  will  bo  inexact. 
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a  griitir  significance  laval  for  the  tasts  Might  ba  warranted. 

AssuMing  tha  linaar  Modal  for  tha  two  adjustments,  ADJUSTMENT  X  and 
ADJUSTMENT  XX *  usad  tha  taaa  aat  of  observations,  then  one  can  test  if 
a  sat  of  parameters*  9,  has  any  significant  affect  on  tha  adjustment  by 
parforaing  tha  adjustment  with  and  without  tha  Q— parameters  and  form  a 
test  statistic*  F ,  based  on  tha  degrees  of  freedom  of  tha  adjustments, 
tha  number  of  additional  constraints*  and  tha  values  of  the  V^PV. 

Let  ADJUSTMENT  X  ba  performed  including  tha  Q— paraaatars  and 
ADJUSTMENT  XX  be  performed  excluding  the  9— parameters.  Then  the  null 
hypothesis  that  all  Q— parameters  could  be  set  to  xero  is  rejected  if 


VtPV 


11 


-  VtPV 


1  —  >  F 

c  m 


(6.39) 


where  r  is  the  degree  of  freedom  for  ADJUSTMENT  X  and  s  is  the  number 
of  Q-parameters. 

Using  this  concept,  a  test  was  performed  to  see  if  there  was  a 

significant  difference  between  ADJUSTMENT  B  and  ADJUSTMENT  C  at  the  5X 

significance  level.  This  test  would  indicate  if  the  periodic  screw 

error  terms  included  in  ADJUSTMENT  C  were  significant.  From 

information  contained  in  Table  S*  the  test  statistic*  F*  was  calculated 

based  on  vVVj  a  3433. SIS*  a  3644. 99S,  r  a  3341  and  s  a  34  which 

results  in  F  a  6.03  with  F_.  ....  _  __  -  1.43*  Therefore,  at  the  5X 

i  ,0.35 

significance  level*  the  null  hypothesis  that  the  periodic  screw  error 
terms  are  not  significant  is  rejected. 


e 


Unfortunately*  there  appears  to  be  no  known  statisieal  test  which 
be  used  to  test  ADJUSTMENT  A  against  ADJUSTMENT  B  or  ADJUSTMENT  C 
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because  ADJUSTMENT  A  involves  e  nonlinear  relationship  of  the  true 
observables,  counter  readings. 

For  very  large  degrees  of  -free don  as  in  these  adjustment,  F— tables 
are  generelly  not  aveilable,  however,  e  fairly  good  approximati on  to 
the  F— distribution  percentile  when  the  degrees  of  freedom  are  larger 
than  30  according  to  Kossack  and  Hensehke  [1975]  can  be  obtained  from 

1°9  Fa  (v,,V2)  *  (  /-R-TT  )  "  eg  (4.40) 

whore  h  *  ^V1V2^V1  ♦  Vj)  end  g  ■  (v2  •  Vj)/Vj\>2  the  w*luM  of  •» 

b  end  c  being  a  function  of  O  .  See  Table  13  for  some  selected  values. 
It  should  bo  noted  that  Vj  and  \>2  represent  the  degrees  of  freedom 
of  the  numerator  and  denominator  respectively. 

It  is  clear  that,  whenever  possible,  the  original  counter  readings 
should  be  used  as  the  observable  in  the  model  for  a  gravity  base 
station  network  adjustment.  One  possible  way  of  accomplishing  this  is 
to  model  the  Calibration  Table  1  information  by  means  of  a  long  wave 
sinusoidal  term.  The  use  of  such  e  model  requires  the  observations  be 
made  over  e  sufficiently  large  range  and  be  adequately  distributed 
throughout  that  range.  Sn  addition,  the  periodic  screw  error  terms 
should  also  be  include  in  the  model  provided  it  can  be  justified  by  the 
distribution  of  the  data.  In  order  to  cheek  whether  a  periodic  screw 
error  term  for  a  particular  gravity  meter  should  be  included  in  the 
model,  graphs  similar  to  those  in  Figure  24  and  Figure  25  were  produced 
for  each  possible  period  of  the  periodic  screw  error  term.  If  the 
distribution  of  the  observations  over  the  period  selected  was  lacking. 
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Table  13  —  Values  of  the  parameters  used  in  the  approximation  of  the 
7— distribution  for  large  degrees  of  freedom. 


Value  of  Cl 

Value  of  a 

Value  of  b 

Value  of 

0.50 

0.0 

— 

290.000 

0.75 

0.5559 

0.55 

0.355 

0.90 

1.1131 

0.77 

0.527 

0.95 

1.4257 

0.95 

0.451 

0.975 

1.7023 

1.14 

0.544 

0.99 

2.0204 

1.40 

1.073 

0.995 

2.2373 

1.41 

1.250 

0.999 

2.4541 

2.09 

1.472 

0.9995 

2.5550 

2.30 

1.557 

as  it  tha  case  In  Figure  28,  than  Including  tha  periodic  acrew  error 
tarn  for  that  period  In  the  pedal  night  net  be  appropriate.  However, 

If  the  distribution  Is  elnllar  to  that  given  In  Figure  24,  then  It 
would  be  appropriate  to  Include  a  periodic  screw  error  tern  for  that 
period. In  the  nedel  far  that  Instrunent. 

4.4  Consistency  gf,  AfcMJJtl *  STM*!!*  YtlWtl 

As  can  be  seen  In  Table  4  and  Table  8,  see  section  4.2,  sone  of  the 
values  of  gravity  dateralnad  by  the  Italians,  Marson  and  Alasia,  differ 
considerably  f ran  the  value  deternlned  by  Hannond  at  the  sane  site. 

The  reason  for  these  differences  Is  not  known  at  this  tine  but  Is 
presently  b«1ng  Investigated  by  Hassaond.  To  show  the  effect  on  the 
station  values  of  the  central  selected  for  a  gravity  base  station 
network  adjustnent.  three  additional  adjustnents,  £, 

I.  •"*  APVmittHT  P»r*orned  using  ADJUSTMENT  C'S 

nodal.  ADJUSTMENTS  D,  I  and  F  used  the  sane  nodal  but  with  difference 
absolute  stations  constrained.  Zn  ADJUSTMENT  D,  the  absolute  stations 
constrained  were  these  deternlned  by  Nannend.  Xn  ADJUSTMENT  E,  the 
absolute  stations  deternlned  by  the  Italians,  Marson  and  Alasia,  were 
constrained.  Zn  ADJUSTMENT  F,  only  four  absolute  stations  were 
constrained.  The  four  stations  selected  to  be  constrained  were  the 
ones  whose  deternlned  values  were  In  good  agreenent  with  both  Hannond's 
and  the  Italian's  values.  A  sunnary  of  the  adjusted  absolute  station 
values  can  be  found  tn  Table  14.  A  sunnary  of  how  the  adjusted 
absolute  station  values  differed  fron  their  Initial  values  for 


Figure  24  —  Distribution  of  •  sot  of  observations  for  B— 131  assuming  a 
period  of  1204/17  counter  units  —  set  no.  1 
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Table  14  —  List  of  the  values  tha  absolute  stations  used  and  the 
results  of  ADJUSTMENTS  D,  E,  and  F. 


XBG/GSS 

Code 

Station 

Name 

ADJ.  D 

Value  Used 

ADJ.  E 

Value  Used 

ADJ.  F 

Value  Used 

119S04 

119S03 

0B1S0C 

1 1926A 

119C03 

119C01 

11994H 

15S05D 

15SV01 

12172A 

155V03 

15221 A 
156E05 
15560A 

McDonald  Obs. 
McDonald  Obs. 

Ml  ami 

Holloman 

Mt.  Evans 
Trinidad 

Denver 

Boulder 

ClBnAM 

97SS200B7 

978828655 

978820097 

978820092b 

979004319 

979139584 

979139600 

979256059 

979330382b 

979598277 

979608601 

979947244a 

979139592b 

979598267 

979608498 

979598272b 

wGspvr 

San  Francisco 

Sheridan 

Boston 

Great  Falls 

B1 smarck 

979972060 

980209007 

980378659 

980497412 

980612882 

980208938b 

980378673 

980497339 

980378666b 

Number 

of  Observat 1 ons 

4484 

4484 

4484 

Number 

of  Parameters 

351 

351 

351 

Number  of  Weighted 
Parameters 

44 

43 

37 

Number 

of  Trips 

837 

837 

837 

Number 

of  Iterations 

2 

2 

2 

V  PV 

3517.475 

3578.143 

3337.026 

Degrees 

of  Freedom 

3340 

3339 

3333 

Aposterlorl  Variance 
of  Unit  Weight 

1.05314 

1.07162 

1.00121 

XQB  —  International  Gravity  Bureau 
OSS  —  Geodetic  Survey  Squadron 

NOTE':  The  variance  -for  each  gravity  value  Is  assumed  to  be  400/lgal 
except  those  marked  with  an  *a'  which  have  a  variance  of  625 /(gal  and 
those  marked  with  a  'b*  which  have  a  variance  of  200  /{gal  .  The  values 
of  gravity  are  given  In  units  of  gal. 
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Table  15  —  Summary  of  the  difference  between  the  adjusted  absolute 

station  values  used  as  control  in  the  various  adjustments 
and  their  initial  values. 


1BG/GSS 

Code 

Station 

Name 

The  value  in /(gal  to  be  added  to 
the  initial  station  value  to  arri 
at  the  adjusted  station  values 

ADJ.  A  ADJ.  B  ADJ.  C 

119S04 

McDonald  Obs. 

-31 

-7 

-7 

U9S03 

McDonald  Obs. 

-30 

-13 

-15 

0S150C 

Mi  ami 

18 

4 

4 

21926A 

Holloman 

-2 

7 

4 

119C03 

Mt .  Evans 

40 

18 

21 

119C01 

Tri ni dad 

-27 

-31 

-26 

11994M 

Denver 

35 

29 

27 

15505D 

Boulder 

29 

21 

20 

155V01 

Casper 

-22 

-27 

-26 

12172A 

San  Francisco 

1 

-5 

-7 

15221 A 

Boston 

-11 

-6 

-6 

15560A 

Bismarck 

-14 

0 

0 

ADJ.  D 

ADJ.  E 

ADJ.  F 

119S04 

McDonald  Obs. 

-1 

-29 

-4 

119S03 

McDonald  Obs. 

-14 

---  - 

— 

08150C 

Mi  ami 

— ■ 

3 

. 

1 1 926 A 

Holloman 

11 

-9 

-3 

119C03 

Mt.  Evans 

29 

— 

— 

119C01 

Trinidad 

-21 

— 

— 

11994H 

Denver 

55 

36 

12 

155 05D 

Boulder 

-42 

24 

— 

155V01 

Casper 

9 

— 

— 

12172A 

San  Francisco 

— 

21 

— 

155V03 

Sheridan 

14 

-23 

— 

15221 A 

Boston 

1 

4 

-5 

156E05 

Great  Falls 

-15 

-33 

■  ■■ 

15560A 

Bismarck 

!■-  — 

8 

...i  ii 

1GB  —  International  Gravity  Bureau 
OSS  —  Geodetic  Survey  Squadron 

NOTE:  The  initial  station  values  for  ADJ.  A,  ADJ.  B,  and  ADJ.  C  can  be 
found  in  Table  8.  The  initial  station  values  for  ADJ.  D,  ADJ.  E,  and 
ADJ.  F  can  be  found  in  Table  14.  Values  are  given  for  only  those 
stations  that  were  used  as  control  in  the  adjustments. 


ADJUSTMENTS  A,  B,  C,  D,  E,  and  P  can  ba  saan  in  Tabla  15. 

Tasta  similar  to  tha  onaa  uaad  to  compare  ADJUSTMENTS  B  and  C  ware 

uaad  to  eoa^are  ADJUSTMENTS  D  and  F  and  to  compara  ADJUSTMENTS  E  and  P. 

Tha  taata  Mara  uaad  to  indieata  if  tha  addad  constraints  introduced  by 

tha  additional  abaoluta  sites,  ware  conaiatant  with  the  four  common 

abaoluta  aitaa  uaad  in  ADJUSTMENT  P.  Uaing  tha  information  in 

Tabla  14,  at  tha  5Z  significance  laval,  tha  test  statistics,  PI,  for 

comparison  of  ADJUSTMENT  D  and  P  and  tho  tost  statistic,  F2,  for 

comparison  of  ADJUSTMENT  E  and  P,  ware  computed.  The  values  obtained 

mere  PI  s  25.747  and  F2  s  40.138  with  F^  ^333  q  c.r  =  2*°  and 

F,  n  ac  a  2,1.  These  results  indicate  that  tha  additional 

o ,  3333,0*95 

absolute  sites  constrained  in  ADJUSTMENT  D  and  E  were  not  consistent 
with  tha  four  absolute  sites  conatrainad  in  ADJUSTMENT  P. 

Tha  affect  on  tha  adjusted  station  values  becomes  apparent  when  the 
difference  between  tha  three  adjustments  are  compared.  Figure  26  shows 
tha  difference  in  tha  station  values  of  ADJUSTMENT  E  -  ADJUSTMENT  D. 
Similarly,  Figure  27  and  Figure  2S  shows  tha  difference  in  station 
value  of  ADJUSTMENT  D  -  ADJUSTMENT  P  and  ADJUSTMENT  E  -  ADJUSTMENT  F 
respectively.  It  is  apparent  from  Figures  26—28,  that  some  sort  of 
linear  trend  with  respect  to  tha  value  of  gravity  exist  for  these 
differences.  The  linear  trend  could  be  caused  by  some  nonlinear  scale 
factor  relationship  f«r  the  gravity  meter  which  is  a  function  of  the 
gravity  at  the  station  or,  more  likely,  by  an  inconsistent  set  of 
absolute  station  values.  As  with  any  gravity  base  station  network,  any 
error  in  the  value  of  the  absolute  stations  is  absorbed  directly  into 
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he  scale  factor  terms  for  the  gravity  meters  used.  The  detection  of 
tad  or  inconsistent  absolute  values  is  very  difficult ,  especially  when 
here  might  be  more  than  one  errant  station  value.  For  example,  as  can 
te  seen  in  Figure  24,  Miami  and  Bsimarck,  two  stations  whose  values 
tere  determined  only  by  the  Italians,  appear  to  be  inconsistent  with 
he  other  absolute  sites.  But  from  the  information  in  Table  IS  for 
ADJUSTMENT  E,  these  two  stations  appear  to  be  consistent  with  the  other 
stations  determined  by  the  Italians.  The  apparent  inconsistency  of  the 
ibsolute  determinations  makes  it  very  difficult  to  check  how  well 
>arious  mathematical  models  for  the  gravity  meter  behavior  perform 
tecause  the  control  for  the  gravity  network  is  provided  by  the 
nconsi stent  absolute  gravity  station  values. 
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27  —  Plot  of  th#  difference  between  Adjusted  station  values  for 
ADJUSTMENT  D  -  ADJUSTMENT  F 


CHAPTER  SEVEN 


C0HCLU8I0NS 

To  obtain  tha  valua  of  gravity  for  stations  In  a  gravity  baso 
station  network  ostabllshod  by  gravity  tlas  Made  with  LaCosta  A  Romberg 
*6*  gravity  meters  requires  that  an  adequate  nodal  be  produced  that  can 
be  used  to  represent  the  gravity  meter's  behavior.  For  networks  having 
a  sufficiently  large  nunber  of  stations  spanning  a  wide  range  of 
gravity  values,  such  as  the  United  States  Gravity  Ease  Station  Network 
which  includes  the  Mid-Continent  Calibration  Line,  a  nodal  has  been 
developed  which  will  pernlt  the  actual  counter  readings  to  be  used 
directly  as  the  observables.  This  nodal  is  based  on  the  ability  to 
approximate  the  gravity  meter's  Calibration  Table  1  Information  by 
naans  of  a  first  order  polynoni al  tern  plus  sinusoidal  terms.  The 
first  order  polynoni al  tern  and  long  wave  length  sinusoidal  tern,  when 
deternined  for  the  range  of  the  observations  nade  for  the  United  States 
Gravity  tase  Station  Network,  resulted  In  a  RMS  of  the  least  square  fit 
to  the  gravity  meter  Calibration  Table  l's  values  in  nilllgal  of  better 
than  10  ygal  for  nost  of  the  gravity  meters  used  in  the  network.  The 
computer  program  used  in  the  adjustments  was  developed  so  a  linear 
drift  rate  term  could  be  Included  in  the  mathematical  model.  However, 
since  physical  explanation  for  the  drift  in  a  gravity  meter  is  a  series 


of  small  tores  securing  at  random  times,  the  drift  term  was  not 
included  in  the  adjustment  model  used. 

Mith  this  type  of  model ,  including  additional  sinusoidal  terms  is 
relatively  easy.  These  additional  terms  are  used  to  represent  any 
periodic  screw  error  effect  that  might  be  present  in  a  gravity  meter. 
Due  to  the  construction  of  the  gravity  meter's  gear  box,  there  is  the 
possiblity  that  sinusoidal  terms  having  periods  of  1206,  1206/17, 
124/17  and  1  counter  units  could  exist  for  gravity  meters  having  the 
old  gear  box.  For  those  gravity  meters  having  the  now  gear  box, 
periods  of  220,  220/S,  22/S  and  1  counter  units  are  possible.  In 
addition,  periods  of  half  those  values  are  also  possible.  However,  it 
became  apparent  that  the  moat  likely  cause  of  a  periodic  screw  error 
effect  in  the  instrument  results  from  the  contact  between  the  jewel 
press  fit  in  the  measuring  screw  and  the  spherical  metal  ball  at  the 
end  of  the  lever  linkage  assembly.  If  the  fit  between  those  two  parts 
is  not  perfect,  a  type  of  periodic  effect  could  oasily  be  introduced. 
The  period  of  this  effect  would  be  1206/17  or  220/3  counter  units 
depending  on  which  gear  box  was  Installed  in  the  gravity  meter. 

This  study  indicated  that  when  a  periodic  term  having  a  period  of 
1206/17  counter  units  was  Included  in  the  model  for  the  gravity  motor, 
amplitudes  as  large  as  IS  ygal  were  found.  Howovor,  most  amplitudes 
were  below  10  ygal.  This  is  far  less  than  the  values  predicted  by 
Harrison  and  LaCoate  [1978]  of  35  to  §0  ygal.  This  does  not  mean  that 
periodic  screw  error  effects  as  large  as  35  to  50  ygal  do  not  exist. 

It  only  means  that  of  the  14  different  gravity  meters  used  in  this 
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study  for  which  thsrs  appeared  to  bo  sufficient  observations  to  permit 
solving  for  tho  periodic  seraw  orror  torn,  none  exhibited  any  largo 
periodic  screw  error  amplitudes. 

This  study  also  Indicated  that  the  behavior  of  the  gravity  meters 
appears  to  be  quite  stable  with  respect  to  time,  changing  only  when 
there  was  an  actual  change  In  the  components  of  the  meter,  such  as  when 
the  long  lover  Is  replaced.  It  cannot  be  concluded  that  the  linear 
sealo  factor  changes  from  year  to  year  as  has  boon  hypothesized  since 
there  Is  no  evidence  of  this  occurring.  However,  the  selection  of 
which  absolute  sites  are  used  as  the  control  for  tho  network  has  a 
groat  Influence  on  the  value  of  the  linear  sealo  factor  determined. 

This  is  because  an  error  In  the  value  of  any  absolute  station 
propagates  almost  directly  Into  the  determination  of  the  linear  scale 
factor  term.  There  Is  no  substitute  for  good  control  In  a  gravity  base 
station  network. 

This  study  indicates  that  the  estimated  accuracy  of  observations 
made  with  a  LaCoste  6  Romberg  *8*  gravity  meter  Is  around  0.02  counter 
units.  The  accuracy  of  the  abselute  station  determinations  used  In 
this  study  Is  probably  closer  to  20  ygal  than  to  their  reported 
accuracy  of  around  10  Ug«l. 

For  owl  sting  gravity  base  station  networks  for  which  a 
variance-covariance  matrix  for  the  station  values  exist,  an  algorithm 


loped  which  could  be  used  to  determine  the  gravity  tie  that 
e  made  which  would  Improve  the  network  the  most  In  the  sense  of 


minimising  the  trace  of  the  resulting  var1anem»cover1ance  matrix  of  the 


gravity  stations  used  In  the  adjustment.  This  algorithm  was  not 
applied  to  tho  Unltod  Statos  Gravity  Base  Station  Network  because  the 
accuracy  of  tho  absolute  sites  which  were  used  as  control  for  the 
adjustment  did  net  appear  to  be  consistent  with  each  other. 

The  problem  of  matching  station  observations  with  their  proper 
station  site  description  is  real  and  at  times  very  confusing.  The 
system  adopted  for  identifying  gravity  stations,  if  not  uniformly 
followed  by  the  organizations  collecting  gravity  meter  data,  will 
result  in  needless  errors  and  confusion.  The  actual  system  use  is  not 
that  critical  as  long  as  each  station  is  assigned  only  one  unique 
identification  code.  But,  it  must  be  used  by  everyone. 

Another  problem  which  should  be  eliminated  is  that  whenever  a  change 
is  made  in  a  station*s  site  description  form,  the  change  should  be 
noted  with  at  least  the  date  that  the  change  was  made.  Otherwise,  old 
and  possibly  erroneous  information  could  be  used. 

This  study  leaves  many  areas  for  future  studies.  The  understanding 
of  how  the  behavior  of  the  LaCoste  B  Romberg  ’G*  gravity  meter  could  be 
modelled  as  the  range  of  the  observations  gets  larger  needs  to  be 
investigated.  If  the  Calibration  Table  1  information  cannot  be 
adequately  modelled  with  a  single  sinusoidal  term  for  a  large  range, 
then  tho  use  of  seme  type  of  piece  wise  continuous  curve  could  be 
investigated.  Methods  are  needed  that  will  permit  during  the 
calibration  of  the  gravity  meter  detection  of  periodic  screw  error 
effects.  This  would  involve  investigating  periods  other  than  1206/17 
counter  units.  However,  only  periods  which  have  some  physical  reason 
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for  existing  should  be  considered. 

The  modelling  of  the  other  possible  systematic  effects  eaused  by 
such  things  as  changes  in  the  distribution  of  the  mass  of  the 
atmosphere,  changes  in  the  groundwater  level,  and  possibility  of 
instrumental  drift  needs  to  be  studied  further.  However,  the  backbone 
of  any  future  studies  rests  on  the  ability  to  obtain  adequate  control 
for  the  gravity  base  station  network.  A  set  of  consistent  absolute 
station  gravity  values  is  essential;  otherwise,  wrong  conclusions  will 


almost  surely  result 
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APPENDIX  A 


Format  i  or>  Of  Hei  oht 


H»*rix  For  £ 


In  the  ‘formation  of  the  weight  matrix  for  the  gravity  observation 

equations,  time  can  be  saved  In  the  formation  of  this  matrix  by 

recognising  the  pattern  on  the  matrices  used  to  form  It.  For  each 

trip,  the  weight  matrix,  M  ^ »  needs  to  be  formed  from  the  matrix  of 

partials  with  respect  to  the  observation,  B,  and  the  covariance  matrix 

of  the  observations,  £  ,  where 

Lb 

M  -  Bj.  BC  (A.l> 

Lb 

M-1  -  (  B l  B*  )-1  (A. 2) 

Lb 

and  the  structure  of  B  has  the  form 


-  B  • 

n-1  n 


b2  -bj 


0 


0 


bn-2  “bn-l 


b  ,  -b 

n-1  n 


(A. 3) 


where  n  is  the  number  of  observations  and  bj  is  the  magnitude  of  the 

partial  with  respect  to  the  observation,  i. 

Assuming  that  the  observations  made  with  an  instrument  are 

independent  and  have  the  same  accuracy,  then  £  can  be  written  as 

Lb 
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L  *  k  x 


(A. 4) 


Hhsrs 


k  —  is  a  scalar  and  has  the  value  of  the  variance  of  each 
observation, 

I  —  is  an  nxn  identity  matrix. 

Rewriting  equation  (A.l)  as 


M  -  k(  BBC  ) 


( A.S) 


where  the  structure  of  BB1  has  the  form 

h2 ..  2  .2 

Vb2  b2 

.2  h2  .2  2 

b2  b2  b3  b3 


-h2  h2  +h2 

bn-2  bn-2  bn-l 


"bn  1  bn  l+b2 
n- 1  n-1  n 


when  both  sides  of  equation  ( A.B)  is  inverted,  M  becomes 


(A. 6) 


M’1  -  i  (BBt)”1 


(A. 7) 


If  the  parli als  with  respect  to  the  observations  all  have  the  same 


magnitude  L,  then  H  ca 


M  -  kl2Q 


(A.B) 


with  the  structure  of  4  being 
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An  analytical  Inverse  exists  -for  equation  (A. 9)  [Gregory  and  Karney, 
1968.  pp  45—46]  of  the  form 

q"1  -  - - -  C  (A.  10) 

n  +  1 

where  n  is  the  order  of  Q  and  the  elements  of  C  are  given  by  the 
relationship 

(:„  -  Cj,  -  i(n  -  j  +  1)  \$j  ‘A. 11) 

Thus  has  the  form 

m"1  -  - -  C  (A. 12) 

(n  +  1 )kL2 


APPENDIX  B 


Formation  SI  Xilft  Normal  Eq.uft.Uff.D3, 

For  off 1 cl oncy,  the  eontrl butl on  to  the  normal  equations  is  based  on 
performing  only  non-zero  multiplication.  The  non— zero  partlals  with 
respect  to  the  parameters  are  stored  on  a  tile  that  can  be  sorted  with 
a  IE  byte  structure  like 

BYTES  VALUE  TYPE 

.  1—  S  V  —  non— zero  partial  Mlth  respect  -floating  point 

to  the  parameters  number 

9—10  R  —  equation  number  within  a  trip  Integer 

o-f  the  partial 

11—12  C  parameter  number  associated  Integer 

with  the  partial 

1^-14  T  —  trip  number  Integer 

15—  IE  S  —  sequence  number  o-f  partial  Integer 

within  a  trip 

The  mlsclosure  vector,  M,  Is  stored  sequentially  In  an  array  In  trip 
number,  equation  number  order. 

The  file  of  non— zero  partlals  are  then  sorted  on  T,  C  and  R  In 
ascending  order.  Mlth  the  non— zero  partlals  In  this  order,  the 
non— zero  elements  of  can  be  formed  efficiently  and  each  element  formed 
Is  saved  on  another  sortable  file  with  the  same  structure  desribed 
above.  Since  the  structure  of  M~ '  Is  a  block  diagonal  matrix  with  a 
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block  for  each  trip,  the  non-zero  contribution  to  ACM  1  needs  to  be 
computed  only  'for  each  trip  where  i  is  the  trip  number. 

To  compute  A  Si  'a  and  A  Si  ' W  efficiently,  all  non-zero  elements  of 
A  Si  ^  are  sorted  on  R ,  T  and  C  in  ascending  order  and  all  non— zero 
elements  of  A  are  sorted  on  C,  T  and  R  in  ascending  order.  Then  if  row 
by  row,  a  row  of  non— zero  elements  A  Si  ^  were  brought  into  core,  the 
contribution  to  the  correspondi ng  row  of  the  normal  matrix,  N,  could  be 
computed  from  the  non— zero  elements  of  A.  N  is  the  relationship 
A  Si  ’a.  As  each  row’s  contribution  was  completed,  it  would  be  saved. 

At  the  same  time,  the  contribution  to  U,  which  is  the  relationship 
aSi”^M,  would  be  computed. 

Since  the  normal  matrix,  N,  is  stored  in  a  sub— block  form,  a  block 
of  rows  of  the  non— zero  elements  of  A  Si  ^  are  brought  into  core  and  a 
block  of  rows  of  N  and  U  are  formed  and  saved.  Thus  the  largest  array 
required  would  have  (number  of  rows  in  a  block)  times  (number  of 
parameters) . 
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APPENDIX  C 

Computation  £1  I&S  VBrlfnSB  J2£  Peal  duals 

Tho  computation  of  tho  covariance  matrix  of  tho  residuals  requires  a 
huge  amount  of  computational  effort.  However,  If  only  certain 
covariances  of  the  residuals  are  desired,  the  computational  effort 
required  can  be  considerably  reduced  by  taking  advantage  of  the  pattern 
and  sparse  matrices  Involved. 

For  the  combined  mathematical  model  F(X  ,L  >  =  0,  It  can  be  shown 

a  a 

that  from  Its  linearised  form 


AX  +  BV  +  W  *  0  (C. 1 ) 

that  the  least  square  solution  gives  the  residual  vector  as 

V  -  -  P_1BtM"1(AX  ♦  W)  <C.2> 

where 

H_1  -  (  B£  BC  )-1  ( C. 3) 

Lb 

X  -  -  (Atlf1A)"1AtH",W  <C.4> 


By  substituting  equations  (C.3>  and  (C.A)  Into  equation  (C.2>  gives 


IBB 


The  variance  and  covariance  propagation  for  a  linear  function  given 


in  equation  (C.5)  expressed  as 

V  -  G(V)  (C.6) 

results  in  the  covariance  matrix  for  the  residuals  of 

QV  "  *  W  *  ^  *  W  ^  <C*7> 

where 

|g  -  P_1 BtM“1 (  I  -  A(AtM“,A)“1AtM_1)  (C.«) 

-  BP'V  -  M  (C.9> 

Substituting  equations  (C.S)  and  (C.9>  into  equation  (C.7>  yields 

Qv  »  P_1Bt(H"1  -  M_1A(AV1A)"1AtM"1)BP'1  <C.10> 

letting  M  s  M  ^  —  M  ^A<A*M  ^A>  ^ A^M  ^  then  equation  (C.10)  can  be 
written  as 

Qy  -  P^B*  M  B  P_1  (C.ll) 

Exploiting  the  fact  that  P  ^  is  a  diagonal  matrix  and  the  ft  matrix  is 
sparse  and  has  a  definite  pattern(  it  can  be  shown  that  the  covarirnce 


for  a  residual  is  given  by  one  of  the  following  relationships: 


1<j<n 


(C. 12) 


“ijj  ’  k'bj(Sjj  *  2 

°vn  '  kib?  ”n 

°V  ‘  klbn  Vl,n-1 
nn 

°5lj  "  '  kibibj*  mij  " 


(C.13) 


(C. 14) 


(C.15) 


k1b?bj(  mij  +  mi -1 , j-1  "  mii 


i  *  J 

1<i<n 
1  <j<n 


(C. 16) 


klb1bo  "l.n-1 


<C.  17) 


k.b.b  (  m.  ,  , 

1  J  n 


m. 

J 


1<j<n 


where 

k  —Is  ths  variance  of  tbs  observations  -for  a  trip, 
n  —  is  the  number  of  observations  for  the  trip, 

m  —  is  the  sunblock  of  77  associated  with  the  trip. 


(C.18) 


The  M  matrix  can  be  partitioned  as  follows 


t - r  — Lj7  ~i - 

;  !  !  !  m 


(C. 19) 


where  1  is  the  number  of  equations  -formed  -for  a  trip. 

If  the  variance  o-f  the  residuals  is  desired,  then  only  the 
tri diagonal  elements  o-f  M  are  needed  which  cuts  down  on  the  amount  of 
computation  required. 

How  to  obtain  the  tri diagonal  elements  of  N  will  now  be  discussed. 
Remembering  that  M  is  a  block  diagonal  matrix  with  one  block  for  each 
trip  and  the  size  of  each  block  is  equal  to  the  number  of  equations 
formed  for  the  trip  and  that  M  s  m’1  -  M_1  ACAtrT,A>  a’VT1 ,  the  major 


effort  is  in  computing  the  tridiagonal  contribution  to  M  of 


M^ACaST’a)  ■'aV1 


Since  both  (A^M  ^ A )  ^  and  the  non— zero  elements 


of  Ail  are  available  from  the  adjustment,  it  becomes  a  matter  of  what 
is  the  best  way  of  obtaining  the  tri diagonal  contribution. 

By  sorting  the  non— zero  elements  of  A^M  ^  on  their  columns  and  rows 
in  ascending  order,  the  multiplication  required  to  form  the  elements  of 
(A^M-^ A) ”^A*M~^  can  most  efficiently  be  done.  Each  element's  value 
along  with  its  row  and  column  in  < A^M” ^ A)”^ A^M” '  are  saved  on  a 
sortable  file  and  then  sorted  in  column,  row  ascending  order.  To 


compute  the  required  tridiagonal  elements,  the  appropriate  columns  of 
-1  t  -1  -It  -1 

M  A  and  (AM  A)  AM  are  brought  into  core,  the  multiplication 


I 


performed  and  the  appropriate  tridiagonal  elements  are  subtracted  from 
-1 

M  and  the  variance  of  the  residual  is  computed. 


